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Values for the surface pressure and area of molecules 
adsorbed on solid surfaces can be calculated from the 
adsorption equation of Brunauer-Emmett-Teller. Over a 
wide pressure range these calculated values are in approxi- 
mate agreement with the two-dimensional equation of 
state for vapor-condensed films. Since the adsorption 
equation of Jura and Harkins was derived from this same 
equation of state, it is concluded that for a great many 
experimental isotherms, the mathematical agreement 
between the two adsorption equations is so good that 
extremely accurate experimental work would be necessary 
to determine which equation was more nearly correct. 

The Brunauer-Emmett-Teller equation gives calculated 
curves on Gregg-type plots that agree qualitatively with 
those reported by Gregg for multi-layer adsorptions. 
Surface-area values obtained from these curves by Gregg’s 
X-point method are generally considerably higher than 


those that would be obtained by the method of Brunauer- 
Emmett-Teller. 

The limiting area occupied by a molecule (i.e., in a 
mono-molecular layer) as calculated by the method of 
Harkins and Jura appears to decrease with decreasing 
energies of adsorption (i.e., decreasing values of ¢ in the 
Brunauer-Emmett-Teller isotherm). It is suggested that 
this decrease is artificial, rather than real, and results 
from the assumption that the constant a in the two- 
dimensional equation of state does not vary with changes 
in the nature of the solid adsorbent. 

It is not possible at the present time to suggest that one 
of the isotherms is to be preferred to the other on theo- 
retical grounds. In general, preference should be given to 
the one which fits experimental data over the broader 
pressure range. 





SYMBOLS 


a4=constant of the two-dimensional equation of state (2). 
b=another constant of Eq. (2). 
b’=intercept on the ki axis of the straight-line equation 
(18). 
¢=constant of the Brunauer-Emmett-Teller adsorption 
isotherm (1). 
k=constant used by Jura and Harkins to convert A 
values to = values. 
ki=constant of straight-line equation (18). 
k2=another constant of Eq. (18). 
m=slope of straight-line equation (18). 
b= pressure. 
Po= saturation pressure. 
v¥=volume of gas adsorbed per gram of adsorbent. 
%m=Vvolume of gas which would correspond to a mono- 
molecular layer on one gram of adsorbent. 
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vo= volume of gas adsorbed at saturation per gram of 
adsorbent. 

A =constant of Jura and Harkins adsorption isotherm (9). 

B=another constant of Eq. (9). 

M=molecular weight of gas. 

N= Avogadro's number. 

R=gas constant. 

T =absolute temperature. 

V=volume of one mole of gas. 

a=area occupied by each adsorbed molecule in a mono- 
molecular layer. 

«=surface pressure (change in free energy of the solid 
surface attributable to adsorption). 

o=area per molecule. 

2 =specific surface of adsorbent. 
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Fic. 1. Pressure-area plot calculated from Brunauer- 
Emmett-Teller isotherm with c=2; the numbers beside 
the points are relative pressures. 


INTRODUCTION 


HE multi-molecular adsorption theory of 

Brunauer, Emmett, and Teller! has been 
found to be in agreement with a great number of 
experimentally-determined adsorption isotherms 
over a fairly extensive range of pressures. In a 
critical evaluation, Brunauer? has concluded that 
the theory has limited accuracy at very low 
pressures where, adsorbents tend to be strongly 
heterogeneous, and at high pressures where 
capillary condensation can occur, but in the 
pressure range where p/po is 0.05 to 0.35 the 
agreement of theory and experiment is excellent. 
The basic equation 


Vv =CUnppo/{po— Pp) (po— Pp +ep) (1) 


has been found by many investigators to fit the 
observed data within experimental accuracy over 
the pressure range given by Brunauer. Extensive 
use has been made’ of this equation to calculate 
the specific surface of solids from 2,,. 

Jura and Harkins! later developed an alternate 

'S. Brunauer, P. H. Emmett, and E. Teller, J. Am. 
Chem. Soc. 60, 309 (1938). 

*S. Brunauer, The Adsorption of Gases and Vapors 
(Princeton University Press, Princeton, 1943), Vol I. 
Physical Adsorption. 

’P.H. Emmett, Recent Advances in Colloid Science (Inter- 
science Publishers, Inc., New York, 1942), Vol. I, p. 1. 

*G. Jura and W. D. Harkins, J. Chem. Phys. 11, 430 
(1943). 


equation for expressing adsorption isotherms. 
They considered their equation to have a thermo- 
dynamic background, since it was derived by 
combining the two-dimensional equation of state 
for condensed mono-layers 


7=b—a0e (2) 


with the thermodynamic relationship derived by 
Gibbs for adsorption on solid surfaces [Eq. (678) 
of the ‘Equilibrium of Heterogeneous Surfaces’”’ }* 
which, for a two-component system at constant 
temperature, is 

dx =(RT/VZ)(v/p)dp. 
Since 

o=ZV/N2, 

differentiation of Eq. (2) gives 

dm= —(aXV/N)d(1/v). 
By combination of Eqs. (3) and (5), followed by 
integration, the adsorption isotherm_of Jura and 
Harkins is obtained: 
Inp =Inpot+ (ad? V?2/2NRT)[ (1/20?) — (1/v*) J. (6) 
The constant terms in this equation may be 


grouped as: 


A=a>*V?/2NRT, (7) 
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Fic. 2. Pressure-area plot calculated from Brunauer- 
Emmett-Teller isotherm with c=5. 


5 The Collected Works of J. Willard Gibbs (Longmans, 
Green and Company, New York, 1928), Vol. I, p. 329. 





B.E.T. AND 
B=I|npot+(A/v0*), (8) 
=rms. so that 
rmo- Inp=B-—A/v*. (9) 
d by 


Harkins and Jura® found this equation, as well 
as Eq. (1), to be valid for almost all of the very 
large number of isotherms that they studied. 
For nitrogen adsorptions, Eq. (9) generally fits 
the experimental data over a wider range of 
pressures than Eq. (1). 

Since the specific surface is a factor in the 
constant A (see Eq. (7)), it can be evaluated 
from Eq. (9). Harkins and Jura® use the relation 


S=kA! (10) 


state 


(2) 


ad by 
(678) 
es” | 
stant 


for this evaluation. They suggest that k has a 
constant value for any adsorbate at any temper- 
ature and recommend a value for nitrogen at 
—195.8° of 4.06. Specific surface values calcu- 
lated by this method were almost identical with 
those calculated by the method of Brunauer- 
Emmett-Teller. Harkins and Jura® state that 
“this is very remarkable when the great differ- 
ence between the two theories, one kinetic and 
the other thermodynamic, is taken into con- 
sideration.”’ It is the purpose of this paper to 
discuss the agreement between these theories 
and the reason for its existence. 
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Fic. 3. Pressure-area plot calculated from Brunauer- 


unauer- Emmett-Teller isosherm with c= 10. 


*W. D. Harkins and G. Jura, J. Am. Chem. Soc. 66, 
1366 (1944), 


igmans, 
329. 
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Pressure-area plot calculated from Brunauer- 
Emmett-Teller isotherm with c=25. 


EQUATIONS RELATING THE THEORIES 


(Note: A brief summary of this section was 
published, under the same title as the present 
paper, in this journal in November, 1944.)? 

The mathematical relationship between pres- 
sure and adsorption developed by Brunauer, ' 
Emmett, and Teller, as expressed by Eq. (1), may 
be substituted into Eq. (3). The result is: 


dx =(RT/VZ)(cimpo/(po— P) 


X(po—ptep) |dp. (11) 


The integration of this equation has been shown? 
to give as the result: 


w=(RTv,,/VZ) In(po— p+cp)/(po—p). (12) 


Data calculated from Eq. (12) for c=50 were in 
agreement with Eq. (2) over a wide pressure 
range. 

It is possible to determine the relationship 
between wz and o, for isotherms obeying the 
Brunauer-Emmett-Teller equation with widely 
varying values of c. A substitution of Eq. (1) in 
Eq. (4) gives: 


¢=ZTV(po—p)(po—ptcp)/Nincppo. (13) 


The terms in Eqs. (12) and (13) that do not 
vary with p or c may be combined into two 


7H. K. Livingston, J. Chem. Phys. 12, 466 (1944). 
°H. K. Livingston, J. Am, Chem. Soc, 66, 569 (1944), 
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Fic. 5. Pressure-area plot calculated from Brunauer- 
Emmett-Teller isotherm with c= 50. 


constant terms: 
k, = V=/(RT vp, 1n10), 
ke= Nvm/ V2. 


(14) 
(15) 


Then, from Eqs. (12) and (13) the following 
equations can be formulated: 


log(po— p+cp)/(po— p) = kim, 
(po— p) (bo— pt+cp)/cppo= kee. 


(16) 
(17) 














Fic. 6. Pressure-area plot calculated from Brunauer- 
Emmett-Teller isotherm with c= 100, 


LIVINGSTON 


For any given value of c, the solution of Eqs. 
(16) and (17) for a series of relative pressures 
will give ki and koe values that can be examined 
graphically for the existence of the straight-line 
relationship: 


kin =b' —mkoo. (18) 


This, of course, is equivalent to Eq. (2). That 
is, if values of w and o calculated from the 
Brunauer-Emmett-Teller equation are in a 
straight-line relationship, then the basic assump- 
tion made in deriving the equation of Jura and 
Harkins is satisfied, and the pressure-adsorption 
isotherm .can be expressed by either of the two 
equations. 


CALCULATIONS FOR SELECTED VALUES OF c 


Values for kir and koe over a wide range of 
pressures have been calculated for the following 
values of c: 2, 5, 10, 25, 50, 100, and 250. The 
results have been plotted in Figs. 1 to 7. In 
plotting these data, considération has been given 
to the accuracy of experimental adsorption data. 
In practice, the comparison of data determined 


experimentally with a theoretical equation is 
made by applying some test of agreement such 


as a least-square calculation and, using an 
assumed value for the accuracy of the experi- 
mental data, determining the range of agreement. 

Since an assumed experimental accuracy would 
be used in determining the fit of either the 
Brunauer-Emmett-Teller of Jura and Harkins 
isotherms, the region where both appeared to be 
applicable would increase in extent as_ the 
assumed error increased. In plotting Figs. 1 to 7, 
it has been assumed that 7 and o are known 
with an accuracy of +3 percent. The position 
of the calculated points is indicated, then, by a 
cross whose dimensions in the horizontal and 
vertical directions represent 6 percent of kee and 
kim, respectively. On each figure an unbroken 
straight line is drawn which at each point passes 
within the ellipse whose axes are the indicated 
cross. The line actually drawn is the longest 
possible line fulfilling these conditions. In the 
region where the conditions no longer apply, a 
broken line is drawn. 

As an example, it will be seen from Fig. 7 
that within the assumed limits of +3 percent, 
data that rigorously obey the Brunauer-Emmett- 
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TABLE I. 








Relative pressures at 
Constants of which the equation 
straight-line fits the data within 3% 


equation b/Po, b/ Po, 
(18) low high 
~ value value 


0.28 . Lj 0.5 
63 oka : 52 
.99 : ‘ -60 
1.29 , ‘ 62 
1.61 . d 75 
1.95 3. . 72 
2.33 3.66 07 .50 
(These values are based on a recalculation of those given 
in reference 7 and should be considered to supercede the 
earlier results.) 











TABLE II. 








Relative pressures at which the 
equation fits the data within 0.5% 


b/po, low value /po, high value 


0.24 0.36 
14 65 
.14 55 
.07 19 
38 -65 











Teller equation (1) from p/po=0.07 to p/po 
=(0.50, with c=250, would be found, insofar as 
the accuracy would permit judgment, to obey 
Eq. (2) (and therefore the Jura and Harkins 
equation (7))¢ Similarly, data that rigorously 
obeyed the Jura and Harkins equation (7) over 
this pressure range, would be found to obey Eq. 
(1) (if c=250) over the same range, and no 
distinction could be made between the two 
equations. 

A summary of the range of pressures over 
which the two equations agree within 3 percent, 
for the seven c values studied, is given in Table I. 
For those ¢ values for which the range is largest, 
the range of pressures over which the two equa- 
tions agree within 0.5 percent has also been 
determined (Table II). Because of the complex 
form of the line on which the points of Fig. 7 lie, 
two sets of limits are given for c= 250 in Table II. 
In the Harkins-Jura terminology, this corre- 
sponds to the existence of two different con- 
densed phases for the adsorbed vapor. 

If accuracy were better than 3 percent, of 
course, the pressure ranges would be narrowed 
somewhat. The great majority of experimentally 
observed c values are from 25 to 250. For these 
values, the pressure range of agreement is very 
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wide, covering almost all of the range wherein 
the Brunauer-Emmett-Teller equation is gener- 
ally found to be valid. Very accurate experi- 
mental work would be necessary to determine 
which equation was more nearly correct in these 
cases. 


COMPARISON OF THESE RESULTS WITH OTHERS 
RECENTLY PUBLISHED 


Gregg*®!° has considered adsorption isotherms 
from the standpoint of two-dimensional equa- 
tions of state, in the same general manner as 
Jura and Harkins. Gregg’s procedure involves 
the integration of adsorption data and construc- 
tion of zo vs. ry plots, for a number of systems 
for which > was unknown. He has constructed 
such plots for adsorbed films on charcoal, mica, 
sodium chloride, ferric oxide, copper oxide, and 
silica gel. For gaseous films, such plots give lines 
extrapolating to ro = RT at r=0. Gregg assumes, 
after Schofield and Rideal," that for condensed 
films the slope of a r7—7= curve would be a/2. 

Gregg has developed two methods for deter- 
mining surface areas. The first method, applica- 
ble to linear plots, depends on the slope being 
a/Z, as predicted by Schofield and Rideal. The 
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Fic. 7. Pressure-area plot calculated from Brunauer- 
Emmett-Teller isotherm with c= 250. 


*S. Gregg, J. Chem. Soc. 696 (1942). 

10S. Gregg, Proceedings of the Conference on Ultra-fine 
Structure of Coals and Cokes (British Coal Utilisation 
Research Assoc., London, 1944), p. 110. 
(1988) Schofield and E. Rideal, Proc. Roy. Soc. A109, 57 
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Fic. 8. Gregg-type plot calculated from Brunauer-Emmett- 
Teller isotherm with c= 10. 


.second method is used in those cases where 
mo vs. rd plots are not linear. Where multi-layer 
formations occur, parabolic curves are ordinarily 
observed. Gregg!® takes the area corresponding 
to point X (which is essentially the vertex of the 
curve) to be a, and suggests that this value is 
approximately the same as would be calculated 
from the slope below this point. Maggs” has 
found specific surfaces calculated in this way to 
agree fairly well with those calculated from 
Eq. (1) (Emmett’s method’). Gregg!® suggests 
that the integration and plotting of Eq. (1) 
according to his method gives curves that coin- 
cide fairly well with the type obtained for his 
assumed equations of state. 

We have used the data plotted in Figs. 3 and 
6 to form a zo vs. rz plot of the Gregg type. 
The curves (Figs. 8 and 9) closely resemble 
several of those that Maggs” has calculated from 
experimental data. Gregg’s first method, which 
assumes that o=a, requires that for a plot of 
kikowo vs. kix, the slope be keo(=kex=1) (see 
Eq. 19). Since the curves drawn in Figs. 8 and 9 
are not straight lines, it is not proper to apply 
this method, but it can be stated that in neither 
case is a slope of unity approached by an 


appreciable linear portion of the section of the 


2 F,. Maggs, Proceedings of the Conference on Ultra-fine 
Structure of Coals and Cokes (British Coal Utilisation 
Research Assoc., London, 1944), p. 95. 
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Fic. 9. Gregg-type plot calculated from Brunauer-Emmett- 
Teller isotherm with c= 100. 


curves that is shown. If the second (X-point) 
method is used, the areas corresponding to the 
vertex of the curves (indicated by X in Figs. 8 
and 9), as obtained from the original kyr—keo 
data, should be o=a. Actually, the areas are 
o=a/1.16 for c=10 and ¢=a/1.29 for c=100. 
This shows that an estimation of the section area 
(a) by this method would be appreciably in 
error for adsorptions that follow the Brunauer- 
Emmett-Teller equation. 

Emmett" has recently compared the Brunauer- 
Emmett-Teller and Jura and Harkins equations 
directly by solving Eq. (1) for v at various 
pressures and plotting the results according to 
Eq. (9) [i.e., logp/po vs. (1/v?) ]. The calculations 
were made for the following c values: 2, 5, 10, 
25, 50, 100, 250, and 1000. This is equivalent to 
the operation carried out in the present papcr, 
but Emmett’s calculations were made with p—? 
data rather than r—o data. That is, in the 
present paper, Eq. (1) was integrated (by means 
of the Gibbs equation) and substituted into Eq. 
(2), from which the Jura and Harkins equation 
(9) was derived, while Emmett used Eggs. (1) 
and (9) in their original forms. 

It is interesting to compare the two results, 
especially since the same c values were used, 
except for c=1000. In both cases curved lines 


3 P,H. Emmett, J. Am. Chem. Soc. 68, 1784 (1946). 
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were obtained. Emmett has concluded that, 
inasmuch as curved lines are obtained, the state- 
ment in reference 7 that straight lines are ob- 
tained within 3 percent is misleading. This is no 
doubt true, if the significance of the 3 percent 
limitation (which unfortunately was not stressed 
in reference 7) is not appreciated. It should be 
emphasized that the 3 percent error allowed ap- 
plies to r—o data, not p—v data. 

Emmett has stated that his data, as plotted 
in Figs. 1 and 2 of his article, approximately 
coincide with a straight line over an appreciable 
pressure range. Within the pressure range p/po 
=().11 to 0.40, and for the common c values of 
25 to 250, his data show a fairly reasonable 
agreement to exist between the two isotherms. 
This leads to the same conclusion arrived at 
from a consideration of Figs. 1 to 7, namely, 
that very careful experimental work would be 
necessary in many cases to determine which 
isotherm is the more accurate. In comparing 
Emmett’s plots with those given in this paper, 
consideration must be given to the fact that his 
method of plotting is much more sensitive to 
variations in v, since he uses v~* as an abscissa, 
rather than o which is a function of v7. 


SPECIFIC SURFACES CALCULATED FROM THE 
TWO ISOTHERMS 

Emmett" has made the interesting observation 
that cross-sectional areas (section areas) of ad- 
sorbed molecules calculated from the Jura and 
Harkins equation vary with c, if they are calcu- 
lated from the synthetic data obtained by sub- 
stituting Eq. (1) in Eq. (7). His plot shows the 
calculated section area of nitrogen to vary from 
13 to 20 A?/molecule as c varies from 30 to 300. 
This same type of calculation can be made from 
the data of Table I. 

The section area, a, is defined by the equation 


a= VE/Nom=1/ke. (19) 


It bears a definite relationship to the constants 
of Eq. (18). This equation may be put in the form 


ax =b’/k, —(mk2/ki). (20) 


Comparison with Eq. (2) shows that the term 
(mk2/k1) must be equivalent to a. Therefore 


Naa 
RT», 1n10 RT 1n10 


aaV> 





m=aak,= 
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If, as was assumed by Harkins and Jura,‘ a is a 
constant for a given adsorbent at constant 
temperature, the m values in Table | may be 
used to calculate a as a function of c. It will be 
assumed that a=16 A? at c=100 (which is the 
value read from Emmett’s Fig. 3). The results of 
this calculation are plotted in Fig. 10. While 
this plot covers a different range of c values than 
Emmett’s Fig. 3, the two seem to be in close 
agreement. 

Emmett reports a trend in experimental a@ 
values roughly coinciding with the derived curve. 
While this may be a generally valid result, it 
should be pointed out that a may not necessarily 
be a constant for different c values. Harkins and 
Jura made this assumption because it is known 
that insoluble mono-layers on water tend to show 
no change in a when the pH or electrolyte 
content of the water sub-phase is changed. Such 
changes might have much less effect on surface 
forces than would result from changes in the 
solid adsorbent that cause large changes in c. 
Certainly a comparison between adsorptions of 
the same vapor on two adsorbents differing in c 
values by a factor of 10 or more might be 
expected to show appreciably different a values. 
It would seem likely then, that the values of k 
given by Harkins and Jura® for different ad- 














Fic. 10. The cross-sectional area (a) of molecules 


calculated by the Harkins and Jura method from the 
slope of Figs. 1-7 (assuming a= 16 at c= 100). 








624 G. HALSEY AND H. S. TAYLOR 
sorbents might not always have the unique _ if the kinetic theory of these authors, from which of 
values assigned. their equation was derived, is an exact expression th 
of the performance of vapors and adsorbents at th 
POSSIBILITY OF A CHOICE BETWEEN THE moderate relative pressures, for a great many - 
TWO ISOTHERMS : ; 
systems the surface pressure and area will corre- af 
It is conceivable that a mathematical analysis spond to the two-dimensional equation of state 
of the two isotherms would suggest a basis for expressed by Eq. (2), and therefore will be in a 
preferring one to the other. No such choice can agreement with the Jura and Harkins isotherm. ch 
be made as a result of the present analysis. The The nature of the relationship between the ad 
fact that the two equations are both equally in two isotherms has been most recently summar- he 
agreement with many of the experimental iso- ized by Jura and Harkins,'* who concluded that prt 
therms determined during the past ten years “it is apparent that there is no conflict.” It was the 
promises to make such a choice particularly predicted that new developments in theoretical = 
difficult. If it is a fundamental property of and experimental studies of phase transitions be 
adsorbed films under these experimental condi- would improve the agreement. In investigations ign 
tions to obey the two-dimensional equation of of this sort, consideration should be given to the mc 
state of condensed films, then it might be con- mathematical relationship between the equations to 
sidered fortuitous that the equation derived by that has been demonstrated by Emmett™ and wh 
Brunauer, Emmett, and Teller, with appropriate in the present paper. — 
egeaiaterian expresses the adsorption isotherm over 4G, Jura and W. D. Harkins, J. Am. Chem. Soc. 68, 1941 ” 
a considerable pressure range. On the other hand, (1946): 
act 
the 
ple 
THE JOURNAL OF CHEMICAL PHYSICS VOLUME 15, NUMBER 9 SEPTEMBER, 1941 y 
The Adsorption of Hydrogen on Tungsten Powders a 
GEORGE HALSEY AND HucGuH S. TayLor wha 
Frick Chemical Laboratory, Princeton, New Jersey on 1 
(Received May 12, 1947) fore 
inte 
An analysis of the extensive experimental data of Frankenburg on the adsorption of hydrogen of t] 
by metallic tungsten powder between —194 and 750°C has been made, using the Fowler- 
Guggenheim treatment of adsorption phenomena as the mode of approach. It has not been of t 
found possible to interpret the data in terms of interaction between the adsorbed species on a intel 
uniform adsorbent surface. It is shown that the experimental data can be interpreted on the able 


basis of a non-uniform surface without interaction, the heterogeneity of the surface being 

specified in terms of an exponential distribution of sites given by the distribution function I 
N=ce-*!Xm, where x represents the energy difference between the lowest energy state of the 

gas and the lowest energy of the adsorption complex. Such a distribution function yields the 

experimentally observed relation dInp/dln@=constant at constant temperature. Such a F¢ 
heterogeneity is especially significant in the case of tungsten metal which hitherto:has been 

treated as an essentially uniform surface. The quantity xm has been found to vary with the 

temperature. The significance of this variation has been briefly explored. 





COMPREHENSIVE study of the adsorp-_ the cleanliness of the metal surface, and_ the wher 

tion of hydrogen on tungsten powders has hydrogen was adsorbed rapidly at the tempera- the 

been recorded by W. G. Frankenburg,'! and he tures employed. These data form an excellent pend 
has presented a series of isotherms between example of an adsorption equilibrium in a local- F —— 
—194° and 750°. Special attention was given to jzed mono-layer. Frankenburg’s analysis of his hes 
tJ. Am. Chem. Soc. 66, 1827 (1944); 66, 1838 (1944). data led to an empirical relation, the logarithm J Engla 
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ADSORPTION 


of the adsorbed hydrogen varying linearly with 
the logarithm of the gas pressure. The plot of 
these isothermal lines, when extrapolated, ap- 
pears to intersect at a point corresponding to 
approximately 3150-mm hydrogen pressure. 

In the present communication these data are 
analyzed from the standpoint of statistical me- 
chanics, proceeding from the treatment of such 
adsorption phenomena by Fowler and Guggen- 
heim.? Their treatment of localized mono-layers 
presents initially a thermodynamic derivation of 
the Langmuir adsorption isotherm. They empha- 
size the semiquantitative nature of the equation, 
becoming in some cases ‘‘useless owing to the 
ignoring of the interaction between adsorbed 
molecules.’ Their treatment is therefore modified 
to include interaction between the molecules 
which has for its effect an apparently increasing 
resistance to adsorption as the surface sites are 
filled, if the interactions are of a repulsive nature. 

As a first approximation in solving the inter- 
active case, Fowler and Guggenheim consider 
the case in which the molecules have a com- 
pletely random distribution among the sites. 
This treatment is modified to include the effect 
of the interactive forces on the equilibrium 
distribution. This interaction on a _ uniform 
surface is rejected by the present authors, in 
what follows, as an explanation of the isotherms 
on tungsten powders. We shall consider, there- 
fore, the simple random case, because for the 
interactive case the refinement lessens the effect 
of the interaction in the sparsely covered regions 
of the isotherms. As will presently appear, the 
interaction is too small even in the more favor- 
able simple case. 


INTERACTION ON A UNIFORM ABSORBENT 
SURFACE 


Fowler and Guggenheim derive the equation, 


ila” ers (1) 


where p is the gas pressure, fo is a function of 


the condition of the materials employed, inde- 
pendent of 6, the extent of surface covered. The 


*R. H. Fowler, and E. A. Guggenheim, Statistical 
Thermodynamics emeoms ag University Press, Teddington, 
England, 1939), Chapter X, p. 421. 
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term w is a measure of the interaction. The 
differential heat of adsorption (which Franken- 
burg derives from his data) can be obtained by 
inserting Eq. (1) in the Clapeyron-Clausius 
equation (2), which we express in terms of the 
heat of adsorption, g, per molecule. 





. d mS Inpo 


= = (2) 
ee a/T di/T 


The heat of adsorption gq is, in this case, a linear 
function of 6. Such a result is not in agreement 
with the data of Frankenburg who finds a linear 
relation between log@ and g over most of the plot. 

If Fowler and Guggenheim’s treatment is 
formulated to include the effect of higher inter- 
actions than between pairs of adsorbed species, 
instead of writing for the energy of interaction, e, 


«= N,(@w), 


where JN, is the total number of sites, one may 
write 


e= N,(@Pw+6%v+---), 


where the higher terms allow for groups of three, 
etc. That is, the energy of a cluster is not the 
sum of the energies of all the pairs in the cluster. 
For the partial potential, u, Fowler and Guggen- 
heim write 


m 1 OF 6 1 de 
= -_—( )=m —Inao+ (—). (3) 
kT kT\AN4« 1-0 N.kT\ 00 
Here, N4 is the number of occupied sites and dy 
contains the partition function of the materials. 


Inserting the more general energy of interaction, 
the absolute activity, X, is 





Be (28w+367 0+ - -+)/kT 
A=erlkT = ; (4) 
(1—@)ao 





The absolute activity of the gas phase is 
pb h’ 

Rigi nee wermenensenim, (S) 
kT (2xmkT)3j4 


where ja is the partition function for vibration 
and rotation in the gas. Thus, 


@(20w+30" n+ - . “kT (6) 


1-06 


p=pPo 
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where Po is of the form 
(2xm)*?(kT)°"j6(T) 
WAT) 


wr, (7) 





where xo is the energy difference between the 
lowest energy states of the gas and the adsorption 
complex. A discussion of its variability will 
follow. From the expression for p, Eq. (6), 


d |Inp d Inpo 
=k 





—q=k = +20w+3@v+---, (8) 
d(1/T) d(i/T) 
and the slope of —q versus @ is given by 
—dq 
—— =2w+66v+--: (9) 
dé 


which shows that, when the coverage is low, 
—dq/d@=2w. Frankenburg’s curve for q versus 0 
has its maximum slope when @ is small and it 
has no inflexion. It follows, therefore, that w 
must have a positive value, and the coefficients 
of the higher terms, 2, must be negative. 
Up to a rather large value of @ (0.2-0.3) the 
slope of the curve —dq/d@ is very large. If higher 
interactions were suppressed this would be equiv- 
alent to having qg versus @ given by the projected 
tangent to this early slope. With Frankenburg’s 
data the slope is so large that the heat of adsorp- 
tion would become negative long before @ reached 
unity. Actually the tangent would reach zero 
heat with @ in the neighborhood of 0.1 to 0.2. 
Because interactions at low coverage are 
almost exclusively isolated-pair interactions and 
because, when all the sites are covered, there are 
two pair-interactions per atom, it follows that 
the lowest value of @ at which the tangent to the 
initial portion of the curve could cross g=0 
would be, with the help of all the higher order 
interactions, when @=}. Frankenburg’s data 
could conform to this condition did they involve 
low energy repulsions. On the contrary, in order 
_ to account for the extremely steep initial slope 
of the g—@ curve a value of w/kT of ~100 or 
more would be required. The early interactions 
would involve a term e! in the exact treatment, 
a condition which renders a poor argument for 
interaction actually impossible. An alternative 
statement of this conclusion is that pair-inter- 
action is more than sufficient to prohibit adsorp- 
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tion, which is clearly impossible since vigorous 
adsorption occurs on the surface at small 6. 

This argument is deemed sufficient to show 
that we cannot interpret the Frankenburg data 
on the basis of a theory of interaction on a 
uniform surface. It can further be pointed out 
that the isotherms yielded by such a theory are 
clearly not straight lines on a plot of logp versus 
log@. Nor do they tend to meet in a point as 
Frankenburg found. Rather, in common with 
all isotherms with a saturation value, they tend 
to become asymptotic to a line of constant 
coverage. 

Thus, no interaction theory on a uniform 
surface can explain a differential heat of adsorp- 
tion with increased coverage which has an initial 
slope dq/d@ greater than 0.25qinitia. Accordingly, 
the case of a heterogeneous surface must neces- 
sarily be examined.* While, in this case, inter- 
actions between neighboring molecules un- 
doubtedly exist, it is simpler, at the outset, to 
begin by considering heterogeneity alone without 
interactions. It emerges that this ideal case, 
omitting interaction, is adequate to provide a 
reasonable basis for the analysis of Frankenburg’'s 
isotherms. 






















THE NON-UNIFORM ADSORBENT SURFACE 





Consider a surface with variable energy differ- 
ences between the lowest energy state of the gas 
and the lowest energy state of the adsorption 
complex. For.any particular value of x the 
Langmuir isotherm is applied: 


0, Pp p 
1-0, po 




















ae-x!kT 






By setting p)=ae-*/*?, the effect of variable x 
has been separated from the remainder of the 
partition-function ratio. The total coverage is 


given by 
6= J mrdx, 


where NV, is the distribution function of sites, 
usually normalized to §N,d,=1. Solving the 








(11) 








3 There are many such efforts previously made. Set 
Fowler, Statistical Mechanics (Cambridge University Press, 
Teddington, England, 1936), second edition, p. 833. Also, 
Brunauer, The Adsorption of Gases and Vapors (Princeto® 
University Press, Princeton, 1943), p. 74. 
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Langmuir expression for 6,, 


1 


é VIET 4. 


P 
the expression for 6 becomes 


N,dx 


(ef, 
a 
re kT 4 1 
p 


Case 1: N is Constant 


To evaluate the integral it is necessary to 
assume a form for N,. One form that permits 
integration with indefinite limits obtains when 
N, is constant. Then if (—N,dx=1, N=1/Ax 
where Ax is the range of integration, Xupper— Xo- 
Inserting this value in the expression (13) for 6 
vields 


— T4y 
p 


p 1 4g xetay ner 
kT p 
«t+ — be ——__—__. 
Ax a 
1 +—e-xolk gl 


P 


To compute the differential heat, —g, by means 
of the Clapeyron-Clausius equation (2), the 
expression for @ must be solved for p. Two 
approximations, corresponding to the low and 
high pressure approximations in the simple 
Langmuir theory can be made. For low pressure 
the poorest sites are totally vacant, in which 
case (a/p)e-*"*?>>1. The isotherm then becomes 


a 
—e~xo/kT 
kT p 
1—9=— In— 
Ax a 
1 +-—e~(x0t Ax) /kT 


p 





kT pb 
St ats coe In( “ewer pe-aenr), (15) 


Ax a 
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Hence, 


kd |Inp 
—g=——= const. 


. d1/T 


(1 — A) Aye —O-MAx/kT — Ay e-AxlkT 


— —___—_—_—_—_— (16) 
e (1) AX/kT _ Ax/kT 

At high pressure, (a/p)e~%°t4”)/*? becomes small 

and the best sites are covered. Then, 


kT a 
1—9=— n( +-e-witr) 
Ax p 


ae~xolkT 


p= (18) 


el Ax/kT _ 4 


(17) 





and 


kd Inp (1—8)Axet-Oax/e7 
—g=——— =const. —-———__ rm , 
d1/T eB Ax/kT _ 4 


(19) 

These equations, complex as they are, do not 
account for Frankenburg’s data, unless a number 
of segments are arbitrarily patched together. It 
is usual to measure g and assume that this 
function represents the distribution pattern of 
the surface. These equations indicate that this 
procedure is not sound. There is always a com- 
plex relationship even in the simple case when J, 
is constant. An exception may be noted in the 
simplification possible when the two approxima- 
tions which have been used apply simultane- 
ously. In this case, adsorption has no beginning 
or end in either direction and the normalization 
becomes meaningless. It is replaceable by a 
factor m adjusted to make @ emerge in the 
proper units. 


a 
1—0=mkT |n-e-*°/*7, 


p 


(20) 
Replacing a e~**/*? by po, 


p 
1—6=mkT |In—. 
Po 


(21) 


From this logarithmic isotherm (21), it follows 
that 
1-6 
Inp =——+Inpo 


mkT 
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kdInp 1-6 kd Inpo 
 @/T . m~ d(1/T) 





“' 


This is a linear relationship between @ and g. In 
this case, therefore, the distribution and the 
heat have the same functional form. Whenever 
the heat function is not linear, as in Franken- 
burg’s results, analysis of the problem with 
finite limits is difficult. Applicability of the two 
approximations already employed, permits a 
more complex function to be used, which, at the 
same time makes the analysis possible. 


Case 2: Exponential Distribution for N 
With a distribution of the form 


N,, = ce7x!xm (22) 


over the whole range of x, positive and negative, 


+o e~Xx!xm 
o=c f ioremensinilig 
a 


—w 
—e7xskT 4 1 


p 


Inserting Z for (a/p)e-*/*? the expression can be 
written 


ay —kT/xm © F (kT /xm—1) 
0-(-) aT f ——dZ. (24) 
p 0 1+Z 


This integral has the value‘ 7 cosecrkT/xXm. This 
expression has the approximate value x»,/k7 for 
large values of this quantity. This expression 
will not be carried in further equations because 
other terms involving 7 are indistinguishable. 

In the expression for @, there are no constants 
involving p under the integral. Now, if the 
exponential retarding adsorption in the region 
of negative x(e+*/*”) increases more rapidly than 
the available sites (e*/*"), the infinite integral 
will converge and become independent of the 
variable. The condition for convergence is 
xm/kT>1. Hence, 


(23) 


(25) 


kT lxm 
oer () f(RT/xm); 


a 


where f(k7'/xm) is written for the infinite integral. 


4D. Bierens de Haan, Nouvelles Tables 
Definies (Amsterdam, 1867), Table 16, No. 1. 
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From this expression it follows that 


kT 
Iné=— Inp+const. 
Xm 


dinp xm 
d Ing kT" 


(26) 


(27) 


which is the functional relation between p and @ 
found by Frankenburg. Similarly, as found by 
Frankenburg for the main portion of his data, 


kd \Inp 


—q= 28 
Y DT (28) 


=const.+xm logé. 


For some value of 6, g is zero. If the value for 
the constant in the equation for g is independent 
of temperature, then the value of @ that makes 
g equal to zero is the same for all temperatures. 
Extrapolation of Frankenburg’s data supports 
this constancy with respect to temperature. If 
we select, as a normalization condition, this 
value as 6=1, the constant becomes zero and 


— {= Xm logé, (29) 


and 


Xm 
Ind = — Iné+Inpo. (30) 
p tr Po 


When p=» all of the isotherms intersect. 
Frankenburg finds this value by extrapolation 
to be 3150 mm. At the point of intersection, for 
powder 9799, @ has the value S=1000X10™ 
mole He per g tungsten powder. At this point 
all of the sites with positive energies would not 
be filled (unless T=O0°K) but enough of the 
sites of negative energy are filled to make the 
net heat effect zero. Frankenburg did not reach 
this point experimentally. 

Below a certain coverage, 6=0.015, Franken- 
burg found that the differential heat of adsorp- 
tion was constant. This indicates a breakdown, 
in the sparsely covered region, of the hypothesis 
in the present treatment that the best site is 
always covered and that, therefore, —q should 
increase indefinitely. This low pressure region 
doubtless has an effective x and can be treated 
on the basis of a simple Langmuir isotherm. 

Granting the validity of the treatment here 
proposed, it should be pointed out that the slope 
(Inp versus In@), 2.0-2.06, indicates that the 
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adsorption is atomic rather than molecular. This 
is the accepted view. The total effect of this 
modification is that Eqs. (29) and (30) become, 
respectively, 


(p/po)* = Oxm/k? (31) 
and 


—g=2x», Iné. (32) 


The constant xm, which is the energy at which 
the sites are a fraction 1/e of those when x=0, 
is a function of temperature which is approxi- 
mately linear. It is somewhat difficult to recon- 
cile this continuous change in xm over the 
temperature range with the calculations of 
Frankenburg in which he found g to be constant 
below 529°C. He himself indicates the reason 
for the discrepancy when he details the calcula- 
tion of g, in which the slopes are adjusted on the 
basis of his modification of the Langmuir 
treatment. 

In Tabie I are given Frankenburg’s average 
slopes and the variation of x» with temperature. 
The values of x» were calculated by use of the 
formula, based on atomic adsorption, 


kT dnp 
ta —— . (33) 
2 diné 


Since xm was considered constant with respect 
to 1/T throughout the preceding derivations, 
the equations will be re-examined in respect to 
—q. The derivation of the isotherm is not 
impaired ; hence, from Eq. (26) it follows that 


kd Inp 
en a(/T) 


OXm 
=const. +1n0( xs — Tr), (34) 
oT 


The slope of the differential-heat curve plotting 
—q versus |né, 


is still independent of @ but now is a function of 
I’. The constancy of g with respect to tempera- 
ture is so often encountered or implied that it is 
of interest to deduce the most general law that 
will yield it. If xm° is a constant, the differential 
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TABLE I. Variation of d Inp/d In@ and xm 
with temperature. 








Temperature d Inp 
°K d |Iné 





99 (67.4) 
200 (24.1) 
273 16.63 
373 10.81 
473 7.58 
573 5.60 
673 4.17 
773 3.06 
873 2.31 
973 1.85 








equation 


OXm 
pink [— + Xm = Xm" 
oT 


has the solution 


Xm=Xm°+cT or —q=x»m" Iné, (37) 


where c is the constant of integration. If c=0, 
Xm =Xm°, which is the special case of constant xm 
treated earlier. The data of Frankenburg are not 
specially suited to this case but are fairly satis- 
factory over the middle range. The expression 


Xm = 16.6 —0.005(T — 273°) 


where the data are given in kcal. per mole and 
the corresponding expression for the heat of 
adsorption, Q, per mole, 


—Q=2.3X 16.6 logé, 


conform more closely but not entirely satis- 
factorily to the Frankenburg data. The form 
that appears to fit the data best is 


Xm = Xm9+ce In7T, (38) 
whence 


— 9 =xm'—c(InT+1) In. (39) 


Numerically, in kcal. per mole and °K, the 
expression becomes 


T 
Xm = 4.45 —5.0 log ° 
300 


Hence, the numerical expression for Franken- 
burg’s data over the main portion of his iso- 
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therms, would be 


Pmm mMue 

1350 &= 10 J 
In this expression, mH» is moles H» per gram of 
tungsten, powder No. 9799, which was used in 
the high temperature runs. For powder No. 9798 
the term 1.0010-5 must be replaced by 0.360 
x 10~-* because of the smaller surface per gram 
of this powder. 

The variation of x, with temperature indicates 
that it is not alone x which varies from site to 
site on the surface of the tungsten. A term in x 
proportional to T is to be associated with a 
varying entropy of adsorption over the surface. 
The exact meaning of an entropy term may not 
be immediately apparent. If the term in the 
adsorption equation that varies over the surface 
is not e~x/*? but an expression of the form 


e*xe-x (ET ae e~x(l—a T)/k Tr (40) 


[9.90—10 log(7T /300)] 
RT 
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then the expression for In@ becomes 


kT 
In@d=— Inp. (41) 
m(1—cT) 


lf x is independent of the temperature, as just 
assumed, a variation in x, such as was found 


necessary in the preceding must result. Taking 
the values found for the constants into account, 
it means that when the energy of adsorption 
changes by 16.6 kcal., the change in the entropy 
of adsorption is 5 entropy units or 0.3 entropy 
unit per kcal. change in energy of adsorption. 
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The electronic absorption spectra of the SO2, SeQ., and TeO2 molecules have been studied 
under widely varied experimental conditions. The vibrational analysis of the ultraviolet 
absorption region of SO: with center of gravity near 2000A reveals that it is formed by at 
least two electronic systems with the following vibrational constants: 


, 


System a+ ++v9=42,170 cm™, = »)'=963 cm, v2’ =379 cm, x’ ve’ =2 cm". 
System ag: ++ vo =44,236 cm™, 9»! =775 cm, = ve’ = 375 cm™. 


The presence in the same region of a third electronic system, with v2’ = 360 cm™ and »;’ = 845 


cm7!, seems highly probable. 


The character of the vibrational structure of the band systems and a preliminary analysis 
of the rotational structure of several bands have enabled us to obtain some information about 


the molecular structure of the excited states. 


The vibrational analysis of the c-system of SeO» in the visible, and the TeO» bands in the 
visible and ultraviolet, have enabled us to correlate the electronic states of the three similar 


molecules. 


INTRODUCTION 


CONSIDERABLE amount of work has 

been done in the last few years, chiefly by 
Mulliken and his collaborators, in the field of 
electronic spectra of triatomic molecules. During 
the war we were engaged on similar problems 
and started inyestigations on the electronic 
spectra of the SOQ2, SeOs, and TeO: molecules. 
The wartime conditions prevented us from pub- 
lishing our results, with the exception of one 
paper on the b-system of SeO..! The present 
paper deals with the analysis of the a-system of 
SO», in the far ultraviolet (center of gravity near 
2000A), the c-system of SeQ, in the visible, and 
with TeO, bands in the visible and near ultra- 
violet. 


HISTORICAL CONSIDERATIONS 


The a-system of SO. was first studied in 
absorption by V. Henri and by Wieland.? They 
correlated the observed bands with three funda- 
mental frequencies of the ground level and 
assigned the frequencies 960 and 379.5 cm to 
the upper level. In 1933, Chow* performed a new 


* Associate fellow of the Belgian National Foundation 
of Scientific Researches. 

' J. Duchesne and B. Rosen, Physica 8, 540 (1941). 

*V. Henri, The Structure of Molecules, edited by P. 
Debye (Blackie and Son Limited, London, 1932), p. 121; 
K. Wieland, Nature 130, 847 (1932). 

* Tung-Ching Chow, Phys. Rev. 44, 638 (1933). 


analysis and suggested that the bands belong to 
three electronic systems, each of them being 
characterized by a frequency of about 750-770 
cm in the upper evel. 

Chow’s analysis was criticized in 1938 by 
Price and Simpson‘ who considered that all the 
most prominent bands in the region 2300—2000A 
form a single progression with v2’ =380 cm. 

The a-system of SO, was observed by Chow? 
in emission as well as in absorption. Later, 
Lotmar® studied the monochromatically-excited 
resonance spectra in SO, vapor and Kornfeld® 
extended the emission spectra, using electrodeless 
excitation. 

While Chow and Lotmar interpreted the sys- 
tem using an asymmetrical vibration, v;’’, Korn- 
feld succeeded in arranging the observed bands 
using only symmetrical frequencies in both 
states. In view of the application of the selection 
rules and the discussion of molecular configura- 
tions, these divergencies must be cleared up. On 
the other hand, the values for the band heads 
given by different authors are in some cases 
widely divergent, so that new measures of this 
system seemed to be of interest. For the present 
study, we were able to use some spectra taken 
by V. Henri in 1922 and 1930. The dispersion 
~ 4W.C. Price and D. M. Simpson, Proc. Roy. Soc. A165, 
272 (1938). 


5 W. Lotmar, Zeits. f. Physik 83, 765 (1933). 
6G. Kornfeld, Trans. Faraday Soc 32, 1487 (1936). 
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was either 2A/mm or 5A/mm at 2200A, the 
pressure varied from 0.5 to 100 mm, and the 
temperature from — 80° to 500°. The absorption 
lengths varied from 5 cm to one meter. 


VIBRATIONAL ANALYSIS 


In the first column of Table I we give the 
mean value for the band heads obtained on a 
series of plates taken at — 80°C and correspond- 
ing to widely varied conditions of pressure and 
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absorption paths. The dispersion was 2A/mm 
at 2200A, and the difference between individual 
measures did not exceed 3 cm~. Some additional 
band heads measured only on low dispersion 
plates, and marked with an asterisk, are also 
included. For these bands the precision is not so 
good. This is also the case for the last bands 
situated towards the ultraviolet, where accurate 
measures of our high dispersion plates are no 
longer possible. 











TABLE I. 
D. and R. Chow D. and R. Chow D. and R. Chow D. and R. Chow 
(42,170(?) (44,663 44,663 (46,218 46,217 47,707 
\42,183(2) 44,677 44,674 46,227 47,742 
_ 44,885 a 46,353 a: 46,348 47,889 
(42,576 42,875 won 
42,589 42,589 (44,966" 44,978 46,444 47:65 
42,759 \44,980 “44,987 46,464 48,000 
(42,955 42,954 (45,032 45,030 (46,471* 47,472 48,015* a2 48,005 
42,968 —~ 42,966 145,046 = 46,484 46,482 48,083 
43,018 45,160(?) 46,574 a3 46,570 48,121 
(43,138 43,137 45,247* a 46,697 48,355 
43,148 43,145 45,337 45,333 46,704* as 46,704 48,390* a2 48,375 
43,332 43,329 Sas “45,351 46,820* az 46,818 a, 48,724 
i, “43,340 45,399 45,393 46,829 a, 48,781 
43,407 ae - 46,866 46,865 as 49,133 
(43,539 43,538 45,499 as 146,878 «46,874 a: 49,477 
143,552 “43,550 45,601* a, 45,572 46,925* as a: 49,550 
(43,701 7 45,672 47,057* as 47,052 ew: 49,910 
\43,713 “43,714 45,686 47,074 a, $0,248 
(43,919 43,918 (45,710 45,706 47,198 a: $0,302 
143,930 “43,929 45,726 45,725 47,223 
44,086 a 44,085 45,759 os 45,757" 47,232 
— (45,858 45,859 47,257* ay 47,246 
44,236 a2 44,236 45,869 47,391 
44,293 44,293 45,923 47,548 
a 44,305 46,077 47,576 
44,491(?) on 46,085* a2 46,084 47,623 
44,543(?) 44,570 ~ 46,104 
47,627* a2 47,636 
44,603 a2 44,603 46,122 46,118 
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The third column of Table I contains the band 
heads observed by Chow® at room temperature. 
The two sets of values are in general in quite 
satisfactory agreement. 

In 1931, V. Henri? pointed out that he had 
observed 120 band heads between 2445 and 
2100A. A systematic study of our plates reveals, 
however, that many of these band heads are, in 
fact, sub-heads in the fine structure. This point 
will be discussed in detail below. In any case, 
the bands given in the first column of Table I are 
by far the most characteristic of the whole 
system. The analysis of the spectrum shows that 
the occurrence of the 380 cm~ »,’ frequency is 
beyond any doubt, but we cannot agree with the 
statement of Price and Simpson‘ that all the 
prominent bands are members of a single v»’ 
progression. 

The microphotometer record of Fig. 1, corre- 
sponding to the long wave-length end of the 
spectrum, shows clearly that there are at least 
three superposed series with characteristic in- 
tensity distribution around a single maximum. 
Moreover, intensity-temperature variations indi- 
cate definitely that the observed series cannot 
arise from excited vibrational levels in the normal 
state. Thus, the distance between the corre- 
sponding members of two series must be corre- 
lated either to a frequency (fundamental or 
combination) of the upper level or to the 
distance between two neighboring electronic 
States. 


SYSTEM {%- 0 





The essential results of our new interpretation 
are summarized in Tables II to V. In Table II 
the columns represent v2’ series with v2’~380 
cm. Taking into account the constant difference 
of about 963 cm! between bands of the first 


TABLE II. System ai. 








1 


42,170 [43,138 
42,183 \43,148 


406 402 


[42,576 [43,539 
\42,589 \43,552 


379 379 394 


(42,955 (3 (43,919 44,885 
\42'968 {43/930 


376 375 


Poet 44,293 
43,343 44,306 


370 370 


{43,701 (44,663 
143,713 \44,677 


373 369 


(45,032 
\45,046 


367 


{45,399 
145,413 


360 
{ 45,759 





44,491 


f enna 
\44,086 
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TABLE III. System az. 














44,236 
367 


44,603 
377 
on _ 


44,980 (779) 45,759 
373 363 


(45,337 (769) a 


45,354 (768) (46,122 
372 363 


(45,710 (761) Lieve 


(758) 
(45,726 (758) 46,484 
345 


75 395 


46,085 (781) {46,866 


4 47,627 (763) 
146,878 
359 
{(46,444)* 


| (46,404)° 
356 


46,820 


48,015 (766) 


48,390 


$8,781* (769) 49,550* (752) 50,302* 


RY 352 360 


19, 910* 


334 344 


18,724* (753) 19,477* (771) 








* Measured by Chow only. 


two columns of Table II, we obtain the second 
characteristic frequency mentioned by V. Henri 
and correlated with »;’. All the values 963+ 380 
cm! are of course alternative values for »;’. All 
the bands in the first two columns seem to have 
the same fine structure. For instance, they are all 
doublets with the same mean separation of 
about 13 cm™. 

The first band in the second column of Table 
Il shows irregular behavior inasmuch as its 
distance from the next one has an abnormally 
high value of 402 cm~'. On the other hand, the 
doublet separation is definitely smaller than for 
the other bands (it is only 8 cm~ according to 
Chow’s measure). Both kinds of irregularities 
result probably from a perturbation affecting 
the first level of the second series (the first 
doublet of the first series is too faint and accurate 
measures are impossible; thus it is not possible 
to state definitely whether this level is also 


perturbed). The other levels seem regularly 
spaced, indicating an anharmonity x’v2’ of about 
2cm". 

The distribution of intensity in the first two 
columns of Table II is very regular around a 
single maximum occurring for v2’=2 or 3. The 
intensities indicated (between brackets) are 
rough visual estimations, valid only when com- 
paring bands in the same column. It must be 
noted that bands of the first column appear only 
at pressures about 10 times larger than those 
necessary to observe the bands of the second one. 

A third group of bands can be fitted into the 
scheme of Table II as a v;/=2 series. They are 
blended by strong bands of other series (see 
Fig. 1) and it is rather difficult to estimate their 
intensity. Nevertheless, it seems that the maxi- 
mum intensity for the bands of Table II is 
reached for the (1’v2’0) series, and this seems 
normal in view of the assumed ratio »;//1" 
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=963/1151. The following series could not be 
observed because their position coincides with 
strong bands belonging to another electronic 
transition. 

The first column of Table II! contains a very 
strong series, with nearly the same characteristic 
frequency difference of about 375 cm~. It also 
consists of doublets but, contrary to the first 
two series of Table II, the strongest head is 
situated towards the ultraviolet, and the doublet 
separation is larger (17 cm~! on the average) 

The series given in the second column of 
Table III is not so regular, and at first sight it 
seems that it could form a continuation of the 
first one, the bands 46106-46122 following the 
bands 45710-45726. This arrangement, which 
would suppose some perturbations, cannot be 
definitely ruled out. But several arguments can 
be advanced against it. First of all, the characters 
of the bands in the first and second columns are 
somewhat different, the latter being less regular, 
for instance in regard to the relative intensity in 
the doublets. The fine structure seems also to be 
different, and the structures of the bands in the 
second column are on the whole more complex, 
probably as a consequence of superpositions of 
at least two series. A possible confirmation of 
this point of view is given by the observations 
made by Chow (at room temperature) of 4 bands 
at 46444, 46464, 46472, and 46482 cm™', instead 
of two bands 46471-46484 observed by us. These 
4 bands can be interpreted as nearly superposed 
doublets belonging to two series. It is to be 
noted that the first doublet cannot be interpreted 
as arising from excited vibrational levels of the 
lower state. It is included (between brackets) in 
the first column of Table III. 

Assuming the existence of two separate series, 
it is possible to interpret two additional bands 
which otherwise cannot be easily fitted in the 
scheme (bands 46085 and 46820 cm7'). 

The constant Ay difference of about 770+370 
cm~' between the first two columns of Table II] 
must correspond to the frequency v;’. We admit 
with Chow that the most probable value is 
v1'=770 cm-. The next columns of Table III 
contain further bands measured towards the 
ultraviolet, as well as the most characteristic 
bands measured by Chow. In this region the 
characteristic frequency difference of about 770 
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TABLE IV. System a3. 








0 1 2 


45,499 (854) 46,353 (845) 


47,198* 
(846) 46,704 (844) 
(832) 47,889* 
e 
46,574 


351 
46,925 





* Measured by Chow. 


cm is clearly visible on the records given in 
Chow’s paper. It is predominant further to the 
ultraviolet, according to Herreng,’ and to not 
yet published results obtained in our laboratory 
by P. Migeotte.** 

For the system of Table III, the v;’ series are 
much more developed than for the system of 
Table II corresponding to the decreasing v;'/»;"’ 
ratio. 


TABLE V. Absorption bands of SOz arising from excited 
vibrational levels of the ground state. 








Corresponding band 
arising from (0’0’0’’) 


43,919 
43,930 


Arising from: 





d A v? 7 
1411633 (2"°0"0") 


41,807 (1'0’0"") 42,955 
{41,987 
\41,998 


(1 0'0'’) 
(1 erg") 


(43,138 
43,149 


{42,059 
\42,073 


(0’'1 "'0'") 
(0’1''0’) 


{42,576 
\42,589 
( 1 7 1 — 
( 1 orn) 


(42,249 
42,387 


43,919 
43,539 
( 1 "00'’) 
( 1 “—-— ) 


{42,766 
\42,776 


(43,919 
43,930 
43,020 


(O71"0") 43,539 519 








Mean values: #1’ =1,151.5 cm~! from infra-red® v1" =1,151.38 
ve’ = §17.3 v?’= 517.84 


7P. Herreng, Rev. d’optique 15, 413 (1936). 

** According to preliminary results of P. Migeotte, it 
seems probable that a new electronic transition is re- 
sponsible for bands in the 1900-1800A region. 
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The two schemes of Tables II and III contain 
nearly all the strong bands observed by us and 
by Chow, with the exception of the bands 
denoted as in the Table I. On account of their 
great relative intensity, it is difficult to correlate 
the latter with an asymmetrical frequency. For 
this reason, the possibility that they belong to 
a third electronic system must be taken into 
consideration. For this hypothetical third system, 
the scheme of Tablé IV with »,/=845 and 
vo’ = 360 cm— seems to be the most probable one. 

The third column of Table IV corresponds to 
bands observed by Chow only. They cannot be 
explained as arising from excited vibrational 
levels of the ground state. A study further to 
the ultraviolet is of course necessary before 
reaching definite conclusions concerning the 
reality of this system. 

Under the conditions of our experiments, 
bands corresponding to excited vibrational levels 
of the ground state are practically lacking at 
— 80°, but they are easily observed at 20°. An 
accurate measure of these bands is possible only 
on the long wave-length side of the (0/0’0’) 
— (000) band of the a; system, and the mean 
values obtained for several bands at tempera- 
tures ranging up to 335° are given in Table V. 
For v;”’ and y.’’ the mean values obtained are in 
good agreement with the values given by Barker’ 
on the basis of infra-red analysis. Av values for 
the first two bands of Table V indicate a small 
but definite anharmonicity for »,’’. 


ROTATIONAL STRUCTURE 


The dispersion used did not allow a complete 
resolution of the fine structure, but some very 
characteristic features of this structure were 
observed. For instance, the coarse K-structure 
was analyzed for several bands situated at the 
long wave-length end of the spectrum, where 
blendings are less frequent than in other parts of 
the spectrum. 

It can be seen easily that the structure of all 
bands in the first two columns of Table II is 
similar, thus reinforcing our previous assumption 
that these bands belong to the same electronic 
transition. On the other hand, we have observed 
two series of sub-bands (p andr) in the K- 


8 E. F. Barker, Rev. Mod. Phys. 14, 198 (1942). 
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structure of several a;.bands; this indicates that 
they belong to the perpendicular type. Using the 
second differences, we obtain C’—C”=-—0,7 
+0.10 cm. Taking into account the value 
C’”’ =1.8 cm~'*** and the fact that the K-struc- 
ture is degraded towards the red, it can be 
stated that C’=1.1+0.10 cm™. 

With the dispersion used, the J-structure re- 
mains completely unresolved, but there are 
indications that the convergence is slow. 

For this reason, it can be reasonably admitted 
that the change of B correlated with the elec- 
tronic transition is small; even the sense of the 
degradation in the J-structure is not clear, but 
it seems to be rather degraded in the same 
direction as the K-structure, i.e., toward the red. 
This statement is enhanced by the lack of 
marked absorption on the short wave-length side 
of the K-heads, even at high pressure. 

It should be of great interest to resolve the 
fine structure of the a2 bands given in the first 
two columns of Table III. All that can be said 
for the time being is that it seems to be very 
different from the structure of the a, bands 
(Table I1). Our previous arguments in favor of 
separating the observed bands into at least two 
electronic systems are thus reinforced.f 





























INTERPRETATION OF THE RESULTS 





The molecular orbitals for Cz, symmetry are 
distributed among four classes, following the 
symmetry of the wave functions with respect to 
the planes of symmetry of the molecule. These 
classes have been described by Mulliken as 4, 
a2, b;, and be. On this basis and by the application 
of the molecular orbital theory, Mulliken® has 
proposed a scheme of the electronic levels o 
angular molecules and particularly of SO». He 
has given a preliminary interpretation of the 
observed electronic spectra of this molecule and 












*** In our previous letter (Nature 157, 692 (1946)) we 
had erroneously given the value C’” =1.6 cm™. This erro 
does not affect our results. ; 

tIt must be stated, however, that the alternative 
arrangement, according to which the bands of Table ! 
and Table III form a single electronic system, thoug! 
rather improbable cannot be definitely ruled out. In fact 
from a strictly numerical point of view it would be possibk 
to consider the first two columns of Table II and the 
columns of Table III as consecutive v’ series with 71’ it- 
creasing from 965+380 to 1062+375 and reaching 119! 
+370 cm for 71’ =2. 
°R. S. Mulliken, Rev. Mod. Phys. 14, 204 (1942). 
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has considered the a-system as resulting from an 
allowed transition 3b.— 5a. 

In the first part of this paper we have ex- 
plained that the ‘‘a-system’”’ must be divided into 
at least 2 electronic transitions corresponding to 
Tables II and III and described, respectively, 
as a, and as. 

The molecular constants of the upper state of 
the two systems are the following: 


ay: ++ vy’ = 963 cm—, vo’ =379 cm, 
x’ ve! =2 cm™, vp = 42,170 cm. 

ag? vy’ =775 cm, ve’ =375 cm, 
vo = 44,236 cm. 


The high transition probability for both sys- 
tems and the lack of antisymmetrical vibrations 
indicate that they are presumably due to allowed 
transitions; it must be stated however that the 
a, bands are much more intense than the a, 
bands. As already remarked, the a; bands appear 
to be of perpendicular type; according to group 
theoretical rules they must thus correspond to 
one of the three transitions: a;—4a;, a;—6,, and 
a,—bs. The first assignment is partially in agree- 
ment with the theoretical scheme proposed by 
Mulliken.tf In fact this scheme predicts two 
transitions in the ‘2000A region, one of them 
corresponding to 4a,—5a; (1 type), the other 
to 3b2—5a, (|| type). It seems thus possible to 
correlate a; with the first and a. with the second 
one. The relative intensities of both systems are 
also in agreement with Mulliken’s prediction. 
This interpretation, however, cannot be definitely 
accepted because Mulliken predicts Aa>0O for 
the 4a,—5a, transition, while, as we shall see 
below, we have deduced Aa <0 from our analysis 
of rotational structure. 

It is also to be noticed that, while the a 
system is weak relatively to a2, it is much 
stronger than the b-system at 3000A, corre- 
sponding to 1a,—2,; this seems also contrary 
to Mulliken’s predictions. 

For this reason, the other possibilities men- 
tioned above cannot be disregarded. 

As the transition a,—}6, has already been 
correlated with the c-system (A max~3700A),° 
we have to emphasize only the third one, i.e., 

tt We are very indebted to Professor Mulliken for his 
valuable suggestion concerning this possibility, erroneously 


neglected by us in our previous letter (Nature 157, 692 
(1946)). 
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la2—4b2. However this would require, according 
to Mulliken’s scheme, 10 ev, instead of about 
5 ev, according to our analysis. But it must be 
noted that even the energetic position of 3b. has 
not been determined with certainty experi- 
mentally and that according to Mulliken all the 
calculations have not much more than a quali- 
tative meaning. If a, corresponds really to 
la2—4b2, it should be necessary to reconsider 
theoretically the variation of 4b. with respect to 
the apex angle a of the molecule. 

On the other hand, we have said that B’ must 
be nearly equal to B”’ while C’ is very different 
from C’’. If so, we have a case already discussed 
by Metropolis'® and related to Aa<0 and Ar~0. 
Then from the variation of C and B, we can 
conclude that the angle @ has decreased by 
about 20°, while 7, has slightly increased by 
excitation. . 

These conclusions are in complete agreement 
with the general considerations of Metropolis!® 
and would explain the fact that the a; system is 
not well developed toward high »,’ values. 

For the a2 system, the lack of rotational data 
prevents a more extensive analysis. A transition 
such as 3b.—5a,, as proposed by Mulliken for 
the whole a-system, seems not to be excluded 
and would give rise to parallel bands. In any 
event, the angle a cannot be determined by the 
analysis of the vibrational potential function of 
the excited state as it has been done by Metrop- 
olis'' for the b-system. The reason lies in the fact 
that the value of a is very sensitive to the cross 
terms for angular molecules and, as it has been 
proved by one of us,” the cross terms cannot be 
neglected in such molecules. Furthermore, as 
the parallel character of the bands in the b- 
system must be considered with caution, we 
shall understand that the assigned transition 
(1a2.—2b,) cannot be considered as being estab- 
lished with certainty. The full discussion of the 
scheme proposed by Mulliken must thus be 
postponed until new experimental data farther 
to the ultraviolet will be available. On the other 
hand, it seems that the az system is much more 
developed in the sense of high v;’ values than 
the a; system, and this indicates that in the 

10N. Metropolis, Phys. Rev. 60, 283 (1941). 

11 N, Metropolis, Phys. Rev. 60, 295 (1941). 


2 J. Duchesne, Mémoires Soc. Roy. Sci. (Liége) 1, fasc. 
2, 429 (1943). 
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former r, must be modified by excitation as 
well as a. 

Before new steps in the interpretation of the 
SO, spectrum can be undertaken, it would be 
important to check the reality of the third a; 
system of Table IV, to get more information 
concerning the distribution of intensity in the 
bands, and to have more reliable data concerning 
bands in the 2000—1800A region. 

P. Migeotte in our laboratory is now engaged 
on such experiments. 


SeO, 
Introduction 


The electronic diffraction patterns of SeO, 
vapor, studied by Palmer and Elliott,'* indicate 
an angular symmetrical structure of the SeO, 
molecule, with a Se—O distance of 1.61A. The 
diffraction patterns being relatively insensitive 
to the apex angle, its value has not yet been 
determined. 

The SOz, SeQ2, and TeOe molecules belong to 
the same class of 18-electron molecules, and their 
comparative spectroscopic study presents, there- 
fore, considerable interest. Up to now, three 
electronic systems of SeO2, have been reported 
between 2000 and 5000A. Preliminary descrip- 
tions of these systems were given by several 
authors,'* and in a recent paper! we have ana- 
lyzed in some detail the vibrational structure of 
the b-system situated between 2300 and 3300A. 
In the present paper we deal with the vibrational 
study of the c-system situated in the visible, and 
we will try to develop some general conclusions 
regarding the electronic structure of the excited 
states. The analysis of the a-system in the 
farther ultraviolet is now in progress and will be 
published later. 


Vibrational Analysis 


For the present study we have used a two- 
prism Bourguel spectrograph for the visible, 
giving a dispersion of about 30A/mm at 4500A, 
and a Jobin quartz spectrograph giving a dis- 
persion of about 15A/mm at 3000A. The source 


13K, J. Parlmer and N. Elliott, J. Am. Chem. Soc. 60, 
1389 (1938). 

“4S. F. Evans, Nature 125, 528 (1930); Choong-Shin 
Piaw, C. R. Acad. Sci. Paris 202, 127 (1936); zhid. 203, 
239 (1936); Ann. de physique 10, 191 (1938). 
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tungsten-strip filament bulb. The absorption of 
the SeO. vapor was investigated in 4-meter, 
25-cm, and 2-cm absorption cells at temperatures 
varying between 300 and 500°C. 

The SeQz crystals were carefully purified by 
repeated distillations in vacuum and introduced 
in a side tube sealed to the main absorption cell. 
This side tube could be heated independently 
from the furnace heating the main cell, and thus 
the vapor pressure could be varied in the usual 
manner independently of the temperature. We 
used pressures between 50 and 450 mm. 

By varying the conditions of pressure and 
temperature, the SeO, absorption spectrum was 
investigated between 5000 and 3400A. In this 
region about 100 bands were measured, and the 
measurements are in general in agreement with 
those of Evans and of Choong-Shin Piaw." 
These authors have reported (especially in the 
ultraviolet part of the spectrum) more bands 
than we. Some of them must probably be corre- 
lated with rotational structure, as we have 
already noted in the case of the b-system.! 

The absorption in the bands of the c-system is 
considerably weaker than in the bands of the a- 
and b-systems. The number of molecules needed 
to reveal absorption in the visible is at least 10 
times greater than the number needed to provoke 
absorption in the ultraviolet. 

A characteristic feature of the c-system is the 
diffuseness of the bands in the central region and 
the increased sharpness toward both ends. On 
both extremes, the bands are relatively well 
defined and the majority of the heads can be 
measured with good precision, although errors 
up to 20 cm~ cannot be excluded for fainter 
bands. In the central region, heads cannot be 
measured at all, and we have determined the 
position of the maxima of intensity on micro- 
photometric records. Another striking feature of 
the c-system is provided by long series of bands 
extending through all the spectrum with the 
characteristic periodicity of about 200 cm™. 
These series are given in the columns of Table V1. 
The 200-cm™ progression must surely be corre- 
lated with the bending frequency v2’ of the 
excited state. The presence of long v2’ progres 
sions indicates a considerable change of the apex 
angle of the molecule. The apparent lack o 
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TABLE VI. 











23,462 
192 
23,650 
184 
23,834 
190 
24,030 
200 
24,230 
192 
24,422* 
200 
24,622 
198 
24,820 
182 
25,002 
203 
25,205 
185 
25,390 
245 
25,635 
181 
25,816 





22,930 
206 
23,136 
209 
23,345 
165 
23,510 
194 
23,704 
185 
23,889 
201 
24,090 
204 
24,294* 
186 
24,480 
210 
24,690 
212 
24,902 


885 


890 


880 


896 


931 


880 


911 


20,201 


20,615 
205 
20,820 
215 
21,035 
203 
21,238 
210 
21,448 
207 
21,655 
188 
21,843 
207 
22,050* 
~~ 190 
22,240 
174 
22,414 
216 
22,630 
163 
22,793 
198 
22,991 
202 
23,193 
197 
23,390 
182 
23,572 
190 
23,672 
197 
23,955 
208 
24,163 
205 
24,368 
196 
24,564* 
211 
24,775 
188 
24,963 


21,680 
192 
21,872 
200 
22,072 
211 
22,283 
199 
22,482 
214 
22,696* 
198 
22,894 
208 
23,102 
208 
23,310 


919 


21,176 
215 
21,391 
194 
21,585 
198 
21,783* 








* The underlined bands correspond to intensity maxima in the series. 
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TABLE VII. 








_¢-system 
weighted values 


b-system 
mean values 


894 911 
893 898 
907 889 
902 903 





3,596 3,601 








stretching frequencies means that the valency 
bond remains practically unchanged. It seems, 
however, not excluded that some fainter bands 
in the ultraviolet part of the spectrum are due to 
a weak excitation of the valence frequency. 

The first column of Table VI contains the 
most prominent bands of the whole system, and 
must, therefore, undoubtedly be correlated to 
the zero level of the normal state. The intensity 
distribution in this series is regular, and the 
maximum occurs at about 24,500 cm. The 
assignment of the other series to definite vibra- 
tional levels cannot be done without some 
ambiguity. This is due to lack of regularity in 
the intensity distribution and to frequent irregu- 
larities in the distances between neighboring 
bands. From a strictly numerical point of view, 
several arrangements are possible; however, the 
scheme given in Table VI seems to us to be by 
far the most reliable; in particular, it gives for 
the ground level a frequency of about 900 cm—", 
already assigned to the stretching frequency of 
normal SeOz on the basis of our analysis of the 
b-system. 

It is important to check this last conclusion 
by comparing the mean values of AG” resulting 
from the study of the b- and the c-systems. The 
agreement of the mean AG” values, correspond- 
ing to both systems and given in Table VII, can 
be considered as satisfactory and is even im- 
proved if we consider only those bands which 
can be measured with good accuracy. In doing 
so we obtain for the symmetrical stretching 
frequency in the ground level of SeO, the value 
of »;/’=910 cm. 

In adopting the vibrational scheme of Table 


TABLE VIII. Isotopic constitution of selenium. 








Atomic weights 74 76 77 78 80 82 
Abundance 0.9 9.5 8.3 24 48 9.3 
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VI for the c-system, we must admit that the 
irregularities in the progressions are due to some 
kind of perturbation. On the other hand, the 
occurrence of two intensity maxima observed in 
the bands of the second and third columns of 
Table VI (underlined bands) must be correlated 
either with the special form of potential energy 
surfaces or with the superposition of at least 
two series. 

We have already noted that the band edges 
become increasingly sharper toward both ends 
of the spectrum ; whereas on the long wave-length 
side the bands are shaded toward the red, they 
are shaded to the opposite sense at the other 
extremity of the spectrum. In the central part, 
the sense of the shading is not well defined; 
however the change of shading can be approxi- 
mately located at 22,000 cm“. 

A similar anomaly in the shading was observed 
by us in the TeO, bands, as will be reported in 
detail in the last part of this paper. On the 
other hand, the bands of all known SO» band 
systems are uniformly degraded toward the red 
(if we neglect the J-structure which has not been 
observed in the spectra of SeOQ, and TeQz). It 
seems to us that the above-mentioned behavior 
of the SeO, and TeO: spectra must be explained 
as due to an isotopic shift and can be correlated 
with the fact that Se and Te have several 
abundant isotopes, whereas S has practically 
only one. We give in Table VIII the isotopic 
constitution of selenium: 

Table VIII shows that the center of gravity of 
the isotopic group lies on the side of greater 
atomic weights. Thus the groups of isotopic 
bands are apparently shaded toward the red on 
the long wave-length side with respect to the 
origin, and toward the violet at the opposite 
extremity. In conformity with this point of view 
the change in the sense of degradation occurs at 
about 22,000 cm-', and coincides with the origin 
of the system according to the proposed vibration 
scheme. 


Interpretation of Results 


The present data do not permit stating 
whether the c-bands of SeO2 correspond to an 
allowed or to a forbidden transition. The in- 
tensity distribution in the vibrational series indi- 
cates strongly that the electronic excitation does 
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not change appreciably the Se—QO distance, 
whereas the apex angle is considerably modified. 

We have already noted that most of the bands 
are rather diffuse, and this could presumably be 
due to the fact that the J- and K-structures are 
shaded in opposite directions. Not much more 
can be said for the moment about the structure 
of the molecule in the upper level. 

It is worth while pointing out that no analogy 
can be established between this system and the 
c-system of the SO, molecule, studied by Me- 
tropolis and Beutler.’5 For the former, Ar=0 
and Aa#0, while for the latter Ar#+0 and 
Aa=0. On the contrary, the c-system of SO: is 
much more similar in regard to these properties 
to the b-system of SeO: for which we have 
previously proved! that 7, is the only parameter 
to change appreciably. 

As Metropolis remarked,!® the ratio of the 
moments of inertia J./J, is practically constant 
(~1.17) in the normal state of a series of angular 
molecules such as SOs, NOs, and ClOs. If we 
extend the validity of this rule to SeQs, it 
provides an apex angle of about 130°. On the 
other hand, we observe that the ratio of the 
antisymmetrical and symmetrical frequencies of 
the same molecules also follows such a sort 
of regularity and equals approximately 1.17. 
Though the identity of the constants seems 
purely fortuitous, the extrapolation of this rule 
to SeO2 gives us the possibility to locate v;’’ at 
about 1060 cm. The bending frequency »,” 
can also be approximately estimated ; in fact the 
apparent disappearance of this vibration in the 
spectrum may be explained by assuming that 
the fundamental frequency »2”’ is a multiple of 
v', the most probable value being v2’’ =400 cm-. 
lf this estimation is correct, the columns of 
Table VI could result from a superposition of 
several types of progressions, and the apparent 
lack of the 400-cm- progression, as well as the 
anomalies of intensity referred to above, could 
be satisfactorily explained. 

With the three fundamental frequencies thus 
obtained (910, 400, and 1060 cm-) and by using 
a valence-deformation force field, we compute 
an apex angle of about 160°. The reason for the 
divergency with respect to 130°, which seems 


tone “ Metropolis and H. Beutler, Phys. Rev. 57, 1078 


much more plausible, could be explained by the 
neglect of the cross terms in the potential func- 
tion. The introduction of a small stretch-stretch 
cross term (f12) of about —0.4 10° dyne/cm is 
sufficient to remove the divergency. On the 
other hand, by taking into account a stretch- 
deformation cross term (g,) only, the best value 
equals about +0.410°5 dyne/cm, and it is 
impossible to obtain quite satisfactory results. 
However, in order to obtain such results, it is 
sufficient to add in the potential function a 
slightly negative value of fis. Thus, as we have 
no reason to believe that fiz has a considerable 
value, while g.=0 (or the reverse), it seems 
established that fiz and g. have opposite signs, 
the former being negative and the latter positive, 
with numerical values smaller than 0.4105 
dyne/cm. Furthermore, the presence of the cross 
terms proves that the apex angle of molecules 
like SeO. and SOz cannot be determined by 
assuming a valency force field, contrary to the 
assumption by Metropolis for the case of the 
c-system of SO». 


TeO, 


The TeO, absorption bands were first studied 
by Choong-Shin Piaw.'® He arranged a part of 
the observed bands in a vibrational scheme 
according to the formula: 


v= 256364+545v,'’ —v;'2 + 206v,’ 
—_ 3u,"2 +67 5v,’ — 2. 703". 


Several objections can be raised against this 
analysis: 

1—It seems quite abnormal that at the high 
temperature used, the frequencies of the ground 
level do not appear in the spectrum. 

2—only about 60 percent of the observed 
bands fit into the scheme. 

3—the selection rules are not respected. It 
thus seemed necessary to re-investigate the TeO. 
spectrum before any attempt to compare it with 
the SO. and SeO. band systems could be under- 
taken. 


Experimental Procedure 


The TeO. vapor reacts with quartz at high 
temperatures necessary to cbtain sufficient vapor 


16 Choong-Shin Piaw, C. R. Paris 201, 1181 (1935); 
Thesis, Paris (1937). 
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pressure. The absorption was therefore studied 
in thin open porcelain tubes filled with carefully 
purified TeO, crystals. The tube with TeO, was 
introduced into the electric furnace after the 
desired temperature was attained. We have 
thus operated in conditions far from equilibrium 
and contrary to the experiments concerning the 
SOz and SeO, absorption we were not able to 
vary the pressure independently of the temper- 
ature. In some cases, it was necessary to avoid 
excess oxygen and then the absorption tube was 
provided with water-cooled quartz windows and 


TABLE LX. Intensity maxima in the bands of 
the b-system of TeQOs. 











21,950 24,680 (29,460 
22,130 24,810 29,490 
22,240 25,000 \29,515 
(22,380 {25,150 (29.650 
22,390 25,170 29,675 
\22,400 {25,190 {9,700 
(22,560 25,340 29,870 
| 22,570 25,470 29,950 
22,580 25,525 30,080 
30,165 a 
22,650 (25,620 30,320 
22,665 { 25,640 30,400 a 
22,680 {25,660 30,515 
25,825 30,610 a 
26,000 30,750 
22,750 26,080 30,790 
22,765 26,160 30,830 a 
22,780 26,280 30,960 b 
26,460 31,050 a 
26,640 31,160 b 
22,920 26,820 31,250 a 
22,930 26,920 31,385 b 
22,940 27,070 31,450 a 
27,095 31,580 b 
23,060 127,120 31,660 a 
27,260 31,800 b 
23,200 27,460 31,865 a 
27,550 32,000 b 
23,385 27,730 32,090 a 
(27,880 32,205 b 
(23,550 { 27,910 32,215 a 
| 23,575 {27,940 32,410 b 
23,600 28,125 32,500 a 
28,210 32,620 b 
(23,710 (28,570 32,700 a 
{ 23,730 { 28,385 32,830 b 
{23,750 (28,400 33,040 b 
28,540 
23,870 28,860 
(24,000 (29,020 
24,020 29,030 
{24,040 29,045 
24,200 29,190 
(24,360 29,345 
{ 24,375 
{24,400 
24,540 
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filled with nitrogen to prevent rapid distillation 
of TeO:, outside the high temperature zone. 
The dispersion used was the same as in the 
case of SeOs. For the bands situated near 4000A, 
it was possible to use a multiple-prism Jobin 
spectrograph in Littrow mounting, either with 
5 quartz prisms giving a dispersion of about 
9A/mm, or with 5 glass prisms with a dispersion 
of about 2A/mm at 4000A. 


Vibrational Analysis 


The most characteristic band system of TeO, 
is situated in the blue and near ultraviolet, 
between 3000 and 4550A, with a maximum near 
3550A. In the conditions of our experiments 
(absorbing layer of a few centimeters, rapid 
distillation from the middle of the furnace 
towards the end) it appears at about 900°. At 
about the same temperature, we have observed 
very strong absorption bands in the region of 
2000A, but owing to the presence of Os, it has 
been up to now impossible to investigate these 
in detail. It is not even excluded that this 
absorption is due to impurities. Choong Shin- 
Piaw observed in the ultraviolet above 3000A 
only continuous absorption. At somewhat higher 
temperatures (about 1000°C) very weak and 
diffuse absorption bands appear in the region 
between 2450-2650A77T). It may belong also to 
some impurities. At a temperature of about 
1100°, some weak and very diffuse bands appear 
between 4650 and 5100A. We shall refer to them 
at the end of this paper. At still higher tempera- 
ture, the TeO. bands gradually disappear, and 
they are replaced by TeO absorption bands which 
will be discussed in detail elsewhere. In this 
paper we will deal chiefly with the main system 
between 3000-4550A, which will be designed as 
the b-system. 

In the b-system, the bands situated in the 
blue are sharp and shaded towards the red. 
Towards the ultraviolet they gradually become 
more diffuse, and even the sense of degradation 


is reversed in the extreme ultraviolet end of the 


spectrum. This change of shading is probably 
due to an isotopic effect in the same way as 


tit The wave numbers of the band heads, which are 
shaded toward ultraviolet, are the following : 40740, 40260, 
39760, 39400, 39260, 39060, 38920, 38740, 38580, 38420, 
38240, 38080 cm™. 
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described for the corresponding case occurring 
in SeOe. In fact, the proportion of the different 
isotopes of Te is such that the general character 
of the isotopic shift is the same in both cases. 
The general character of the b-system of TeO: 
has many common features with the b-system of 
SeO». As in the case of SeOs, it is rather difficult 
to locate the heads accurately, owing to the 
diffuse character of most bands. The spectra 
taken with great dispersion show a very compli- 
cated sub-structure which, even at the dispersion 
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OA, 
»bin 
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640 





30,320 
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of 2A/mm, is not resolved into lines. The present 2 

eOr data are not sufficient to decide whether this . 

let, observed structure is due to isotopic effect or 3 

and represents sub-bands in the coarse structure of = 

ents the system. Probably both kinds of substructure Ri 

pid interfere and explain the complicated character 3 

navel of the whole system. ee 

A The positions of the most prominent band = 2 

ved heads are in good agreement with those given by | A 

n ol Choong Shin-Piaw, but for the reason given a < . 

— above these figures must be considered with 2 E z 
_ caution. For the vibrational analysis, we pre- ‘Ss ” ” 

“ ferred to use chiefly the measures of the band = 3 

od maxima given in Table IX. A le 

het As in the case of the b-system of SeQz, the < [5 

pale" TeQ, b-system is characterized mainly by long * ig 

ston series of bands with a periodicity of about 650 2 i 

ads cm~'. Several irregularities are observed in the a a 

— distances between neighboring bands; these ir- | 

pear regularities are less marked if we use band 2 

leoon maxima instead of band heads. With small dis- EF: 

— persion the complicated structure of the bands |8 

and resolves into very characteristic doublets, with $ 

hich irregular doublet separation. As this doubling is i" 

this presumably due to some kind of perturbations, lz 

stem we have used the mean values for the vibrational 


660 


d as analysis. In the Table IX these doublets are 
joined by brackets. The vibrational analysis 
the finally adopted is reproduced in Table X. Most 
red. of the bands of Table IX are fitted into this 
ome scheme with only one characteristic frequency in 
tion each state; both are to be considered as sym- 
the, metrical valency frequencies v;/=670 and py,” 
ably =810 cm“. 
y as lt is worth while to note that the value 


25,340 


660 


660 24,680. 


24,020 


v1 =810 cm~! for the valency frequency of the 
ground state of TeOs is nearly equal to the 
vibrational frequency w”’ =796 of TeO given by 
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Choong-Shin-Piaw.'*t The same relation is cor- 
respondingly valid for the valency frequencies 
of SeO. and SO, (910 and 1152) as compared 
with those of SeO and SO (909 and 1124). 
Although at present, no theoretical considera- 
tions supply a basis for such a relation, it seems 
nevertheless that the parallelism in the behavior 
of the six molecules constitutes a further argu- 
ment in favor of our analysis. The scheme of 
Table X contains nearly all bands observed 
between 22,000 and 30,000 cm-. In the ultra- 
violet end of the spectrum, however, there are 
two rather long series of faint bands which 
appear only when the absorption in the central 
part is complete and which cannot be fitted in 
the scheme. These bands are marked a and 0 in 
Table IX and form two progressions with a 
characteristic frequency of about 200 cm~. On 
the other hand, some of the weaker bands 
farther to the visible, which cannot be fitted in 
the scheme of Table X, are situated at about 
200 cm= on the violet side of strong bands 
belonging to the main series. It seems reasonable 
to admit that 200 cm represents the deforma- 


tion frequency of the upper state and that each 


t According to the new analysis of TeO bands (not yet 
published) by one of us, w’’=790 cm™. 
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v;' level is combined with a series of v2’ levels, 
This statement cannot, however, be considered 
as definitive, as the possibility that the 200-cm~! 
progression belongs to a separate electronic 
transition cannot be ruled out. If all the bands 
belong to a single system, it would mean that 
the corresponding transition is accompanied not 
only by considerable change of the apex angle a 
but also of the Te—O distance. 

After having established the analogy existing 
between the b-systems of SeO2 and TeOs, we 
have searched for a system of TeQO: situated 
farther toward the red in analogy with the 
c-system of SeOQ,. At the highest pressure ob- 
tainable, corresponding to temperature above 
1100°C, some weak and diffuse bands, without 
marked heads, appear between 4650 and 5100A 
on microphotometer records; we have measured 
the following maxima: 21,520, 21,030, 20,560, 
20,150, 19,750 cm—. These bands bear no analogy 
with the c-system of SeQse. 

Our best thanks are due to Professor R. S. 
Mulliken for valuable criticism. Some of the 
spectra used in this investigation were taken and 
measured by R. Migeotte and A. Robert; we 
wish to express to them our thanks for their 
effective help. 





levels. 
idered 
)-cm~! 
‘tronic 
bands 
n that 
ed not 
ingle a 


xisting 
Je, we 
tuated 
th the 
ire ob- 
above 
vithout 
5100A 
-asured 
20,560, 
inalogy 


rR. 3 
of the 
en and 
ort; we 
wr their 


THE JOURNAL OF CHEMICAL PHYSICS 


VOLUME 15, 


NUMBER 9 SEPTEMBER, 1947 


Restricted Rotator Thermodynamic Properties from the Old Quantum Theory 
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Accurate partition functions and free energies are calculated for the hindered internal 
rotator from the old quantum theory based upon the quantization rule of Wilson and Sommer- 
feld. The description of the system is incomplete and the energy levels may even be incorrectly 
placed, but the partition functions are surprisingly close to those now accepted as correct. 
In the course of the development the lowest energy state of the harmonic oscillator turns 
up in its correct position at the half-quantum level. 





IGH precision can be achieved in calculating 

the contributions of vibration and pure 
rotation to thermodynamic properties by meth- 
ods which omit most of the detail of the energy 
states of the molecule. Vibrations are taken to 
be harmonic, interactions are assumed to be 
negligible, and the contributions are expressed 
by simple equations involving not more than 
one molecular constant for each degree of 
freedom. 

Before these simplifications could be used 
with confidence they had to be justified mathe- 
matically or by reference to experimental data. 
For example, Mayer, Brunauer, and Mayer! 
derived the energy levels of the unsymmetrical 
three-dimensional rotator and the free internal 
rotator of the ethane type from the wave me- 
chanics, and showed that the partition functions 
could be represented by summations involving 
only moments of inertia and fundamental con- 
stants. Their partition functions proved to be 
the same as those obtained by the classical 
statistics. Consequently, it was possible for 
Kassel? to use the classical method when he 
developed the general case for the rigid frame- 
work with attached symmetrical rotators in 
extension of the important generalization of 
Eidinoff and Aston.* Kassel’s result proved for 
practical purposes to be the same as had been 
predicted by the writer.‘ 

Another less sweeping simplification is applica- 
ble to the rotator in regions where the classical 


J. E. Mayer, S. Brunauer, and M.G. Mayer, J. Am. 
Chem. Soc. 55, 37 (1933). 

*L. S. Kassel, J. Chem. Phys. 4, 276 (1936). 

1933) L. Eidinoff and J. G. Aston, J. Chem. Phys. 3, 379 


‘J. O, Halford, J. Chem, Phys. 2, 694 (1934). 


statistics are no longer accurate. Here the energy 
levels based upon the old quantum theory give 
the same partition functions as those derived 
from the wave mechanics. For the harmonic 
oscillator the old quantum theory may be used 
in spite of the fact that the energy levels are 
incorrectly located. Consequently, except for 
restricted internal rotation, it appears ‘that the 
finer details of the wave mechanics contribute 
nothing of practical significance to the ordinary 
calculation of thermodynamic properties. 

In this paper it will be shown that a simplified 
method may also be used for the internal rotator. 
The old quantum theory based upon the Wilson- 
Sommerfeld quantization rule yields partition 
functions which do not differ significantly from 
those obtained by more detailed methods. The 
description of the energy states is incomplete 
and the individual levels may be in incorrect 
positions, but the partition functions and, there- 
fore, the thermodynamic properties always as- 
sume values very close to those now accepted 
as correct. 

The old quantum theory handles this problem 
without any assistance and leads to a minor 
uncertainty only in the region of close approach 
to free rotation. The uncertainty in the free 
energy, however, remains within about 0.02 
cal./mole/deg., and only a small and unimpor- 
tant part of the tables of thermodynamic prop- 
erties is affected. 

Even the half-quantum of energy in the ground 
state of the harmonic oscillator is revealed cor- 
rectly by the old quantum theory of the hindered 
internal rotator. This unexpected result is a 
consequence of the necessity for establishing a 
consistent scale of quantum numbers to be used 
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in both the vibrational and the rotational regions. 
Thus one of the major faults of the method 
disappears during the present application. 

It is doubtful if any procedure which had been 
so successful in its day has ever been more 
completely discarded than the Wilson-Sommer- 
feld quantization method. Confidence in its 
results has reached such a low point that nearly 
everyone, including the writer, would be un- 
willing to trust its application to the restricted 
rotator if tabulations derived by methods ac- 
cepted as more rigorous were not available for 
comparison. As far as thermodynamic functions 
are concerned, the present study indicates that 
it has not been wise to discard this useful tool 
so completely. 

It becomes possible to calculate the thermo- 
dynamic properties of a compound, except in 
very unusual cases, without recourse to the 
detail of the wave mechanics. This statement 
has valuable practical implications, because it is 
evident that the complications inherent in the 
solution of the Schroedinger equation can con- 
stitute a serious obstacle and may lead the user 
into making approximations which tend to 
nullify the advantages of the method. 

The effect of varying the assumed potential 
energy can now be studied in a simpler manner 
than the one outlined by Pitzer and Gwinn.° 
There appears to be some theoretical justification 
for a function of the form 


V =(Vo/2)[1+cosn0+k(1—cos2n68) ], 


which corresponds roughly to curve C in their 
Fig. 2. Calculations are- now in progress to 
determine whether a detailed study of this 
function would be justified. 


RESULTS 
Table | shows the differences between the free 
energies of Pitzer and Gwinn and those obtained 


TABLE I. Deviation of old quantum theory free energies 
(cal./mole/deg.) from the results of Pitzer and Gwinn. 








0.1 0.2 


+0.0007 —0.004 —0.003 

—0.0004 +0.002 +0.006 

+0.0006 —0.002 +0.003 $-0.005 

— 0.0002 0.000 —0.001 ' +0.001 0.00 
—0.0003 —0.001 —0.001 J 0.000 0.000 











®K.S. Pitzer and W. D. Gwinn, J. Chem. Phys. 10, 428 
(1942). 
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in this study at twenty-seven points distributed 
over the entire useful range of the variables 
V/RT and 1/Q;. A positive difference means 
that F/T or (F—F;)/T is greater than the 
Pitzer and Gwinn value. V is the height of the 
simple cosine potential barrier, R the gas con- 
stant (1.9869), ZT the absolute temperature, Q,; 
the partition function of the free rotator, 


OQ, =2.8148(10°77)3/n, (1) 


and J is the reduced moment of inertia. 

The largest difference appears at 1/Q;=0.55 
and V/RT=1.0. At this point an adjustment of 
the zero state to insure asymptotic approach to 
the free rotator, as described later, becomes 
most uncertain, but the spread of reasonable 
adjustments should not cause deviations larger 
than 0.02 in the column for 1/Q;=0.55, 0.015 
for 1/Q;=0.50, and 0.01 for 1/Q;=0.45. Within 
each column, it has been found by calculations 
not included in the table that the deviations 
decrease rapidly from the maximum as V/RT is 
changed in either direction. 

For further application of the old quantum 
theory it would be permissible to adjust the 
values in this region toward those of Pitzer and 
Gwinn. For the present, however, they are left 
as they are in order to show what the method 
produces without assistance from other sources. 
In any case the deviations are unimportant in 
relation to the errors of thermodynamic data, 
the arbitrary selection’ of the potential function 
and the probability that the region of the largest 
deviations will not be used unless methyl alcohol 
proves to have an unexpectedly low potential 
barrier.® 

It is obvious without further calculation that 
the other thermodynamic properties will be close 
to those found by Pitzer and Gwinn. 


ENERGY LEVELS 


The total energy E is expressed as the sum of 
the kinetic and potential energies by the equation 


2E=16+V(1—cosn8), (2) 


in which @ is a rotational coordinate and 1 is the 


symmetry number or the number of potential 
minima in the cycle. The momentum P, equal to 


6 J. O. Halford, J. Chem. Phys. 15, 364 (1947). 
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16, becomes 


p= [[2E— V(1—cosn8) }? 
= (IV)(2(E/V—1)+(1+cosn@) }?. (3) 


For states with energy less than V, E/V—1 is 
equal to —(1+cosn@o)/2, where 6 is the ampli- 
tude of vibration. 

The Wilson-Sommerfeld rule requires that 


=  pdd= (IV) f (2B V—-1) 
+ 


(1+cosn@) }'d0, (4) 


with r to have only integral values. The sym- 
metry number can be factored out of the integral 
by the substitution g=76. Limits are then 
assigned for a complete cycle in ¢ to give 


rh=2(1V)! nf (2(k/V—-1) 
: +(1+cosg) dy. (5) 


This equation states that the rotational quantum 
number at a given total energy is one th of the 
quantum number for the same energy in the 
single minimum case and means that only each 
nth rotational doublet is to be counted. The 
corresponding point in the Pitzer and Gwinn 
development is the assignment of one state to 
each of the regions delimited by their theory. In 
Eq. (5), for convenience, the integration is 
performed over half the cycle and the result is 
multiplied by two. 

The quantum number obtained in this manner 
gives the mean position of the rotational doublet. 
If the integration is performed in the vibrational 
region strictly according to the rule, the upper 
limit is changed to the amplitude, go, and the 
length of the cycle is doubled. Here, to keep 
rotations and vibrations on the same number 
scale, the doubling is omitted and the rotational 
scale is retained into the vibrational region, 
where each quantum number r now gives the 
mean position of two vibrational states. But in 
the harmonic oscillator limit the individual states 
are equally spaced on both the energy and 
number scales. The two number scales, vibra- 
tional and rotational, become consistent if the 
vibrational states are placed at r=1/4, 3/4, 5/4, 

‘*, corresponding to the vibrational quantum 
numbers 1/2, 3/2, 5/2, In this way the 
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harmonic oscillator acquires an energy hy/2 in 
the ground state by an argument based upon 
the Wilson-Sommerfeld quantization rule. 

The key to the use of the old quantum theory 
in the internal rotation problem is to retain 
under all conditions the rule that the rotational 
quantum number 7 is the arithmetical mean of 
the quantum numbers of the two components of 
the doublet. 

The manner of arriving at Eq. (5) suggests 
that it should give a solution for the energy 
levels in the single minimum case when there is 
no over-all rotation of the molecule around the 
axis of the internal motion. This interpretation 
is in qualitative agreement with diagrams pub- 
lished by Koehler and Dennison’ for the meth- 
anol case in illustration of their general solution 
for the energy levels with the simple cosine 
potential function. 

For states entirely in the rotational region the 
theory fails to designate the doublet separation. 
This is known, however, at the top of the barrier, 
and can be reasonably assumed to approach 
zero for high energy levels, a situation that can 
be approximately represented by assigning the 
positions r+ V/4E. Because of a tendency toward 
cancellation of the errors introduced by this 
assumption, its effect is felt only in those cases 
where the first doublet is at a short distance 
above the top of the barrier. 

The integration of Eq. (5) has been performed 
graphically at chosen values of go or E/V. At 


TABLE II. Energy of the internal rotator in terms of the 
quantum number index. 





E/V A (radians) p 


0.06724 0.008539 
0.26385 0.033505 
0.07297 


0.5746 

0.9749 0.12348 

1.4305 0.18165 

1.9005 0.24134 

2.3344 0.29644 
0.3394 


2.6729 
2.8284 0.3592 
0.4133 


3.2550 
3.7267 0.4732 
0.5548 


4.3689 
0.6861 


5.4026 
7.0071 0.8898 
1.1959 


9.4175 
12.9503 1.6445 


0.03016 
0.1170 
0.2500 
0.4132 
0.5868 
0.7500 
0.8830 
0.9700 


wipe 
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7J. S. Koehler and D. M. Dennison, Phys. Rev. 57, 
1006 (1940). 
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TABLE III. 6 (Eq. (9)) as a function of p. 








p 6 


0.008539 0.00002 
0.033505 _ 0.00024 
0.07297 ~ 0.00098 
0.12378 0.00215 
0.18165 - 0.00235 
0.24134 0.00000 
0.29644 — 0.00863 











each point from twenty to forty ordinates were 
calculated and the area was approximated for 
the separation Ag and the ordinates yo, yi, ---, 
Yn, according to the equation 


A =Agl[(yotyn)/2+yi+ saa +Yn—1 |. (6) 


Apparent areas A for several values of Ag were 
plotted against Ag and extrapolated to Ag=0. 
For this purpose, if the number of ordinates is 
chosen conveniently, only one addition is neces- 
sary. For example, with forty ordinates the first 
total is based upon yo, ys, etc. To this are added 
Y4, Vio, etc., to give the second total for half the 
spacing first used. The rest of the process is 
evident. Simpson’s one-third rule might have 
been used, but the writer is a little suspicious of 
it and preferred the method just described. In a 
few tests the error made with Simpson’s rule for 
a given Ag proved to be about equal to that 
made with Eq. (6) for Ag/2. Any advantage of 
Simpson’s rule is partly offset by the extra work 
of adding the ordinates in new combinations and 
the need for enough ordinates to provide three 
points for the extrapolation. 

As the ratio E/V increases above 1.0 the area 
rapidly approaches that of a rectangle erected 
on the ordinate at g=2/2 and can be obtained 
accurately with quite large spacing of the ordi- 
nates. At the top of the barrier the integration is 
carried out analytically to give an area of 2! 
radians and provide a check on the accuracy of 
the graphical method. The difference appears in 
the fifth digit beyond the decimal point. Inci- 
dentally the classical picture of a vibration at 
this point is rather amusing. At the maximum 
amplitude the restoring force is zero, suggesting 
that the period should be infinite and the 
frequency zero. 

Equation 5 is combined with the equation 


QO; = (8x°IkT)*/nh (7) 
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to give Eq. (8), in which p is a convenient quan- 
tity that may be called the quantum number 
index. 


r/LOsV/RT)*]=A/(22°)' =p. (8) 


Since the relation between E/V and A or p is 
the same for all systems within the range of the 
variables a single curve or its equivalent in 
tables and equations will summarize the results 
of the quantization. 

Table II shows in successive columns the 
related: values of E/V, the integral A, and the 
quantum number index p. 

The information in Table II can be expressed 
in equations to be used with tables of differences. 
For E/V between zero and 0.75 the equation is 


E/V=3.545p—1.812p?+6 (p<0.24), (9) 


for which 6 is shown in Table III. 
At higher levels the equation 


E/V=2p?+0.5+A (p>0.24) (10) 


is valid. The values of A appear in Table IV. 
For interpolation in Tables III and I[V it is 
convenient to use the graphical method. 

It is evident that the curve of E/V against p 
is asymptotic to the harmonic oscillator line 
E/V=2rn'p in the low energy limit, and ap- 
proaches the parabola E/V=-p?+0.5 at high 
energies. The latter limit is consistent with the 
statement by Koehler and Dennison that the 
high levels approach the free rotator energies 
plus V/2. The relation to the harmonic oscillator 
and the limiting rotator is illustrated in Fig. 1. 

In terms of the. variables used here the 
equations for the harmonic oscillator and the 
free rotator take particularly simple forms. For 
the oscillator the expressions E=2rhv and 


TABLE IV. A (Eq. (10)) as a function of p. 








A 


— 0.01683 
+0.06703 
0.10695 
0.10792 
0.09471 
0.06325 
0.04641 
0.03303 
0.02132 
0.01259 
0.00695 
0.00371 
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y=n(V/2I)'/2m are combined with Eqs. (7) and 
(8) to give E=2zn'pV. For the free rotator, the 
relation E=(nrh)?/8x°I together with Eqs. (7) 
and (8) yields E=zp*?V. Although p becomes 
infinite as V approaches zero, the energy remains 
definite for all finite values of the quantum 
number 7. 

To find the allowed energies for a given ro- 
tator, the quantum numbers r+} and r+V/4E, 
where r is an integer, are divided by the product 
Q;,(V/RT)? (which is independent of the temper- 
ature) to give the required values of the index p, 
after which the energy ratios E/ V are determined 
from Eqs. (9) and (10) and Tables III and IV. 
The sum of the quantities e~/"—#0/”) (V/R7) is then 
the required partition function. 

This simple procedure gives highly accurate 
partition functions and free energies for all cases 
for which three or more individual energy levels 
lie within the potential valley, and covers practi- 
cally all situations for which the thermodynamic 
properties would probably be needed. The special 
treatment used between this point and the free 
rotation limit is considered in a separate section. 

HIGH TEMPERATURE LIMIT 

When 1/Q; approaches zero the spacing of the 
energy states becomes infinitesimal and the 
partition function may be expressed by means 
of the following integral : 


g=2f 


As Pitzer* originally found, the ratios Q/Q; and 
the corresponding free energy differences are 
finite, although the partition functions themselves 
are not. 

The integral from p=0 to p=0.24 was found 
by the graphical method in the manner described 
above in connection with the quantizations. The 
remainder of the integral was evaluated closely 
enough by a rapid procedure which enlisted the 
aid of tabulated values of the probability integral 


DL 


e- HIRT dy = 20;(V/RT)! f e-EIRTdy, (11) 
0 


71 
Pr=(2/n) f exp(—x*)dx, (12) 


found in Mellor’s Higher Mathematics. 


°K. S. Pitzer, J. Chem. Phys. 5, 469 (1937). 
_ "J. W. Mellor, Higher Mathematics for Students of Chem- 
— Physics (Dover Publications, New York, 1946), p. 
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The integral between the limits x; and xz is 
P,—FP,, and for our purposes takes the form 


ap2 
P,—P\= (2a/r!) f exp —a’p" dp (13) 
ap} 


from which, if a=(rV/RT)!, the fraction of the 
partition function of a free rotator which comes 
from the region between p; and p, is given 
directly by P,—P, as found in the table for 
X,=ap; and x2=ap2. If p; and pe are not far 
apart, the correction A of Table IV may be 
regarded as constant, and the contribution of 
the region to Q/Q; will be given by 


[Q/Q; oS as (Po—Py)e—O-5t SVT, (14) 


As the energy ratio approaches the limit 
E/V=rp’+0.5, the fractions may be taken over 
increasingly wide intervals. The result is nearly 
as accurate as the direct graphical integration 
when only twelve areas are summed in this 
manner. The procedure is equivalent to replacing 
the E/V curve by a series of parabolic sections 
parallel to the free rotator curve before per- 
forming the integration. 


APPROACH TO FREE ROTATION 


In the harmonic oscillator limit the ground 
state is found at the position ryo=} or 
po=1/[40;(V/RT)']. If this relation were applied 
in the region of close approach to free rotation, 











y 


iil 




















0.2 p 0.4 


Fic. 1. Ratio of energy to barrier height as a function 
of the rotational quantum number index. A, harmonic 
oscillator. B, restricted rotator. C, limiting rotator. 





650 541 @. 


TABLE V. Energy levels of methyl alcohol for V=770 cm. 








E/V 
(this study) 


0.2598 
0.7144 
1.0329 
1.5241 


E/V (K and D) 
Min. Max. p 


0.2482 0.07594 
0.6784 0.22782 
0.9644 0.37970 
1.2152 0.56196 





0.2506 
0.7187 
1.1697 
1.6766 








with Q; constant and V/RT approaching zero, 
the zero state would appear to emerge from the 
potential valley at a finite value of V/RT, and 
the ratio Ey/ V would increase indefinitely as the 
barrier decreased toward zero. It is easily shown 
that in the limit the energy would be one- 
sixteenth that of the first rotational doublet, and 
that when the zero state was above V/2, the 
partition function would exceed the free rotation 
value. 

It is, therefore, necessary to adjust the zero 
state downward if an asymptotic approach to 
free rotation is to be achieved. While it is 
obvious that the zero-point energy must disap- 
pear in the limit, there is no evident reason why 
the partition function should not exceed the free 
rotation value through some limited range of 
very low potential barriers. 

An examination of the graphical illustrations 
of Koehler and Dennison raises the same question 
in connection with the region over which the 
zero state may roam in the triple-minimum 
problem as the potential barrier decreases toward 
zero. The answer would come from a detailed 
study of the combined partition function of the 
coupled internal and external rotations, and will 
be sought at some future date when time 
permits. 

For the present it is assumed that the zero 
state will remain low enough to insure an asymp- 
totic approach of the partition function to the 
free rotator limit from lower values. An adjust- 
ment is made according to the arbitrary relation 


apo’ +po=1[40,(V/RT)}]=p1/4. (15) 


There appears to be no way to determine a and 
8 uniquely, but satisfactory evaluations can be 
made by considering the behavior of the function 
in two limiting regions. It is required, if 7 is to 
be at the mean position of the doublet, that the 
first term of Eq. (15) shall become negligible 
when the first doublet is entirely within the 


HALFORD 


potential valley, where the separation is fixed 
by the quantization over the vibrational cycle. 
For this situation po is 0.072, and an error of 
about 0.0004 can be made without introducing 
a significant uncertainty into the partition func- 
tion. If 8 is assigned, the corresponding «@ is 
determined by this assumed error. The behavior 
of the function near free rotation for various 
values of B is next examined. It is assumed that 
at 1/Q;=0.55 and V/RT =0.10 the partition func- 
tion will be below and very close to the free 
rotator value, and the lowest integral value of 8 
is taken for which this is true. The equation 
then takes the form 


7.7 X10®p0° + po = pi /4. (16) 


Actually the deviations in Table I are smaller 
for the next higher power of po. A change in the 
permissible uncertainty at pp =0.072 will alter the 
selection of a and 6 without producing much 
change in the size of the corrective term. The 
procedure is quite arbitrary, but no reasonable 
change in the assumptions produces serious 
differences in the partition functions and free 
energies. 

For p; greater than 0.29, Eq. (16) is used to 
find the position of the ground state. If the 
energy at 7o+1/2 falls within the barrier, the 
lower component of the first doublet is placed at 
this point and the second component is moved 
upward to 3/2—ro in order to keep p; for r=1 at 
the mean position. When ro+1/2 is in the rota- 
tional region the doublet components are placed 
at r=12+V/4E. 

As might be expected, the largest deviations 
of Table I occur near the midpoint, on the po 
scale, of the region spanned by Eq. (16). 


METHYL ALCOHOL 


The energy levels given by the old quantum 
theory can be compared with those calculated by 
Koehler and Dennison for methyl alcohol with 
an assumed potential barrier of 769.43 cm. 
They have described the first four levels in 
detail, showing exactly how the components of 
each triplet change in position with the quantum 
number of the external rotation about the axis 
of the methyl top. In Table V, the first and 
second columns show the ratio E/V for the 
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minimum and maximum positions of the states 
within each triplet, and the third, fourth, and 
fifth give the values of r, p, and E/V obtained 
in this study. The reduced moment of inertia 
was calculated from other constants in their 
paper to be 1.116X10~*°. 

At the upper levels the old quantum theory 
result appears alternatively toward the upper or 
lower end of the region delimited by the wave 
mechanics. At the same end of each region, 
when the external quantum number is zero, is 
found a component of a separated doublet 
derived from a rotational state whose quantum 
number in the complete system of levels is a 
multiple of 3, or, more generally, a multiple of . 
The Wilson-Sommerfeld rule, therefore, gives an 


average position weighted toward the true level 
taken by a component of each mth state when 
there is no external rotation. The upward dis- 
placement of the odd numbered states, on the 
vibrational scale of quantum numbers, is greater 
than the downward displacement of the even 
numbered states. 

For the same zero level the old quantum theory 
would lead to a lower partition function than 
the one obtained from the mean positions of the 
regions derived from the wave mechanics. In 
Table V, however, the zero state is a little high, 
and our experience indicates that it is just high 
enough in general to compensate for the different 
description of the higher levels and lead to the 
correct partition functions. 
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A Gaussian network is defined as a network of flexible 
chain segments, linked to each other and to a system of 
fixed points, in which each unbranching chain segment can 
take on a number of configurations which is a Gaussian 
function of the distance between its ends. Real molecular 
networks, such as those of rubber, can under certain cir- 
cumstances be treated as Gaussian networks. The present 
paper carries out a systematic mathematical discussion of 
the statistical properties of Gaussian networks: the total 
number of possible configurations of the network as a 
function of the fixed-point coordinates, the probability of 
finding a given element of the network in a given position, 
or of finding two elements of the network in given relative 
positions, and so on. All probability-density functions 


UBBER-LIKE materials consist of molec- 

ular networks formed by the random 
bonding together of long flexible molecular 
chains.! In developing the theory of rubber-like 
elasticity it is necessary to study the statistical 
properties of these networks—to determine how 
the total numbers of configurations possible for 
the network depends on the external constraints, 


1H. M. James and E. Guth, J. Chem. Phys. 11, 455 
(1943). 


appear as exponentials of quadratic forms, with constants 
explicity expressible in determinant form. An explicit 
reduction to a sum of squares is given for all quadratic 
forms occurring in the theory of coherent Gaussian net- 
works, and an explicit general formula is found for in- 
tegrals of the form 


fi, exif axXe--- fax, exp -= 2 wiXiX;}, 


There is described a mechanical analog of a Gaussian 
network, by consideration of which the statistical proper- 
ties of the Gaussian network can be determined. The 
method is applied to the discussion of the statistical 
properties of a Gaussian network with the connectivity of 
a regular cubic lattice. 


to compute the probability that a given element 
of the network is in a given position, or that a 
given pair of elements are in given relative 
positions, and so on. It has been shown else- 
where ! that for such purposes it is a reasonable 
approximation (in the case of soft rubber) to 
assume that each unbranched chain segment in 
the network can take on a number of configura- 
tions which is a Gaussian function of the sepa- 
ration of its ends, and this independently of the 
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Fic. 1. Terminology. a—fixed points; O—points neces- 
sarily treated as “‘junctions”’ in the present paper; X—points 
which can be treated as junctions. 


configurations assumed by other segments. Such 
chain segments will be termed Gaussian seg- 
ments, and the networks will be called Gaussian 
networks. Attention will be confined in the 
present paper to the statistical properties of 
Gaussian networks. 

In real molecular networks the molecular 
chains cannot take on all configurations inde- 
pendently of each other; they are subject to 
mutual constraints which, for instance, deter- 
mine the volume occupied by the material. It is 
a consequence of these constraints that the most 
probable positions of all elements of the networks 
do not lie at the same point (as they would for 
an unconstrained Gaussian network), but are 
distributed with uniform average density through 
the material; some elements of the network are 
found always near the surface of the material, 
some deep in the interior. It would be exceedingly 
difficult to take detailed account of these con- 
straints. Instead, one can* represent the real 
system by a model consisting of a corresponding 
Gaussian network of which certain elements— 
those corresponding to elements of the real 
system that lie near the surfaces of the material 
—are constrained to lie at fixed points on the 
bounding surface of the model. Change in the 
macroscopic form of the material thus corre- 
sponds to displacement of these fixed points of 
the model. 

The networks to be considered in the present 
paper are Gaussian networks which extend 
between a set of fixed points. The problem at 


? This point will be discussed in more detail in a forth- 
coming” paper. 
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hand is to determine the statistical properties of 
these networks as they depend on the nature of 
the network and on the positions of the fixed 
points. The subject will here be developed sys- 
tematically from an almost purely mathematical 
point of view; physical applications will be found 
in papers quoted here and in others now in 
preparation. 


I, TERMINOLOGY AND NOTATION 


The networks to be considered are formed 
when long chains are bonded to each other and 
to fixed points. The points at which the chains 
are bonded together are usually called junctions 
of the network, and the chains are considered to 
be divided into segments by the points of bonding 
and the fixed points. 

In the present paper it will be convenient to 
extend the meaning of the terms ‘‘segment”’ and 
“junction.” A chain will be considered to be 
divided into segments at points of bonding, at 
fixed points, and at any other points which it 
may be desired to single out for consideration. 
The term ‘‘junction”’ will be applied to all points 
at which segments are terminated; it will thus 
include points of bonding, the arbitrarily chosen 
points of divisions of a chain into segments, and 
the free ends of chains (see Fig. 1). Since the 
division of an unbranching chain into segments 
can be made arbitrarily, the term junction can be 
made to apply to any chosen point of the net- 
work. In the theory to be developed, attention 
will be focused on the junctions and fixed points 
of the network, each of which is connected by 
an unbranched Gaussian chain to one or more 
other junctions or fixed points. It is to be re- 
membered that the statistical properties of 
“junctions” thus derived will be properties 
shared by all points of the network. 

The fixed points of the network will be dis- 
tinguished by Greek subscripts a, 6, ---, the 
junctions by Roman subscripts 7, j, k, ---; the 
Greek letters 7, u, v will apply to both types of 
points. Thus }-. will indicate a sum over all 
fixed points, >; a sum over all junctions, and 
>, a sum over all points of both types. The 
positions of these points will be specified by the 
Cartesian coordinates x,, y;,2,, or by the vector f;. 

The number of configurations possible for the 
chain segment linking the points 7 and » will be 
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expressed, as a function of the separation of these 
points, by the Gaussian segment-configuration 
function. 


(Nv, \r,—r,|) 


=k,, exp{ —(3/2PN,,)|r,—r1,|?}. (1-1) 


If the segment can take on finite numbers of 
configurations when the ends occupy points from 
a discrete set in space—as in some models con- 
venient for the discussion of physical problems— 
then the segment-configuration function can give 
directly the number of possible configurations. If 
the configurations of the segment are con- 
tinuously variable and the ends can take on all 
positions in space, then c(N,,, |r,—r,|) can give 
only the relative number of configurations pos- 
sible for the segment. The configuration function 
cis expressed in the form derived (as an approxi- 
mation) for a segment of N,, independent links 
of fixed length /. The parameter / will be the 





Cit... Tg, ines i r;, 


This function of the coordinates of all fixed 
points and junctions will be called the primary 
configuration function of the network. 

From the primary configuration function one 
can derive the total number of configurations 
of the network consistent with given positions of 
the fixed points by summing over all possible 
positions of the junctions. In the continuous 
case, to which attention will here be restricted, 
the corresponding integral over the coordinates 
of all junctions is 


C(r,, Tg, °° -) 
= fav. favy--Clte ost od Cl 


The relative number of configurations of the 
network consistent with given positions of the 
fixed points and of some one of the junctions, 
say the ith, can ‘be obtained by integrating C 
over the coordinates of all junctions except the 
ith. Denoting such an integration by /dV,’, we 
have for this configuration function for the junc- 
lion i 
€i(r., tg, +++; r;) 


= favictte Ts, °°°Ti, Tj, °°+). (1-4) 


-++)=K exp{ —(3/2P) 2 (1/ Nev) L (Xr —%)? + (91 — Iv)? + (2 — 2)? J}. 


7>yv 


same for all segments; each segment will be 
characterized by the parameter N,,, which may 
be regarded as the effective number of links in 
the segment. The constants k,, will be of no 
further interest. 

When the points 7 and » are not directly 
connected by a segment, one may set 1/N,,=0. 
The nature of the network is completely specified 
by giving the complete symmetric matrix 
(1/N,,); this, and the positions of the fixed 
points, fe, will determine all the quantities in 
which we shall be interested. 

Let the fixed points and junctions of the 
network have definite positions (x,, y,, 2,). The 
relative number of possible configurations of the 
entire. network, consistent with these specified 
positions, is simply the product of the number of 
possible configurations assumed independently 
by the component segments of the network. 
This is 


(1-2) 





The physical importance of this configuration 
function is the following. If a molecular network 
is effectively Gaussian and its potential energy 
is independent of its configuration, then the 
probability dP of finding the ith junction with 
coordinates in the ranges x; to x;+dx;, y; to 
yitdy;, 2; to 2,+dz;, will be proportional to the 
number of configurations of the network con- 
sistent with this position of the ith junction, 
€(r., --+; 1;)dxdydz;. We define the probability 
density function for the junction 1, 


P(e, elite ri), 


such that this probability is 
dP =P (fa, tp, 


-+-;rijdxdydz;. (1-5) 


P; differs from €; only by a normalizing factor 
determined by the condition that 


favre fg, °°°31,)=1. (1-6) 


Since 


faveicre fe, °°°3 Ti) =E(ra, Ts, °°), (1°7) 
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we have 


93 8%) 


=@,(ra, Tg, | 


Pile Ig, 
*s ri)/E(ra, Tp, °° *). 


This probability density function describes. the 
fluctuations in the coordinates of the 7th junction 
of the thermally agitated network; its use to 
determine the time-averaged mean values of the 
coordinates, their average deviation from these 
mean values, and so on, will need no description 
here. 

In the same way one can obtain a con- 
figuration function for the junctions 7 and 7 by 
integrating the primary configuration function C 
over the coordinates of all junctions except the 
ith and jth: 


(1-8) 


@5(ta, Ts, °° °5 Ti, T{) 


-fa Vij’ C(ta; Ea, ** "Ty ty ** -). (1 -9) 


This gives the relative number of configurations 
of the network consistent with the specified 
positions of the fixed points and the two junc- 
tions. From it one can obtain the probability 
density function 
Pij(Tay + °°; Ti, T5) 
=@ (ra, — 
such that 
dP =P;;(Ta, Te, °° 


°5 Ti, r;)/C(ra, Tg, coe), (1-10) 


-;1,, r,)dx dy dzdxjdyjdz; 
(1-11) 


is the probability that junctions i and j of the 
thermally-agitated flexible network will simul- 
taneously have coordinates in the indicated 
ranges. 

From P;; one can determine the probability 
that the relative coordinates of the two junctions, 


E=xX,—Xj, N=Vi-Vjy F=2i—2;, (1-12) 


lie in given ranges. Let P;, be expressed as a 
function of r; and é, 7, ¢, instead of as a function 
of r; and r;. The probability that these coor- 
dinates will lie within ranges of magnitudes 
dx;, dy;, dz;, dt, dn, df, respectively, is easily 
shown to be 


dP=P;;(ta, Tp, °° 5 Tj, 0, §)dx dy dzjd&dnd¢. 


(1-13) 
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The probability that &, , ¢ will lie in the given 
ranges, without regard for the values of the 
other variables, is 


dP’ = R;;(Ta, Ig, ° 
where 
Rij(Ta; 


“+5 & , S)dédndg, (1-14) 


“+36 @ 2) 
= [av Pite sietig’ rj, g, Up o) (1-15) 


is the probability density function for the relative 
coordinates &, 7, ¢. 

In the same way one can derive from the 
primary configuration function other configura- 
tion functions for larger sets of junctions, and 
from these can obtain other probability density 
functions, by normalization and integration over 
unwanted coordinates. Only the functions defined 
above will be of immediate interest here. 

The integrals defined above, and in particular 
the integral of Eq. (1-3), will not all exist if it is 
possible to remove any part of the network from 
the system without breaking any chains. For 
example, if the network includes a set of junc- 
tions linked by chain segments to each other, 
but not to any other junctions or to any fixed 
points, then one can increase the coordinates of 
all these junctions together indefinitely without 
increasing the separation of any linked junctions, 
or decreasing the value of the primary con- 
figuration function C. Integration over the coor- 
dinates of these junctions, as in Eq. (1-3), will 
then lead to non-convergent integrals. A network 
in which every junction is linked to a fixed point 
by at least one continuous path in the network 
will be termed a coherent network. In this paper 
attention will for the most part be restricted to 
coherent networks, for which the integrals thus 
far defined always exist. 


II. TRANSFORMATION OF THE PRIMARY 
CONFIGURATION FUNCTION 


For many purposes it is convenient to express 
the primary configuration function C in terms of 
new variables—the Cartesian coordinates Xj, Vi 
Z; of each junction relative to its most probable 
position. 

Let xo:, Yoi, 20; be the coordinates of the junc- 
tions which, taken all together, make C al 
absolute maximum: (It will be shown in Section 
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V that these are also the mean and the most 
probable positions of the several junctions of the 
network.) These quantities depend upon the 
coordinates of the fixed points in a way which 
we shall now determine. 

The quantities xo;, Yoi, 20 are the values of 
Xi, Viy Zi, respectively, which solve the equations 


dC /dx;=dC/dy;=dC/dz;=0 (2-1) 
for all 7. If we write for the fixed-point coor- 


dinates 
Xa =Xoa; (2-2) 
Eq. (2-1) becomes, by Eq. (1-2), 
Lix(1/Nzi) (or — X01) = 0,) (2+ 3a) 
= +(1/N-:) (yor — yi) =0, >for all 7. (2-3b) 
¥+(1/N 1) (Gor —20:) = 0.) (2-3c) 


Yr=d 1/N,, | 
ver 4 
Yvw=—1/Nn, rxv.J 


(2-4) 


Then Eq. (2-3a) can be written as 


bs VisX0j= — Le ViaXas (2-5) 


with similar equations in the y’s and 2's. The 
terms on the left involve the unknowns Xo; 
whereas the terms on the right involve the known 
coordinates of the fixed points.. These linear 
equations can be solved by Cramer’s rule to 
determine the xo; in terms of the xe. To give 
explicit form to this solution it will be con- 
venient to use a determinant notation which will 
later be employed extensively. 

Let I’ be the determinant of the coefficients y;;, 
with one row and column for each of the 
junctions in the network: 


Yu 


(Note that the fixed points of the network are 
not represented by rows and columns in this 
determinant.) Let I,’ be the cofactor of y,; in T’ 
—that is, (—1)**? times the determinant ob- 
tained by striking out of T the ith row and jth 
column. Let I;,° be (—1)*+*+! times the 





determinant obtained by striking out of T rows 
7 and j, columns k and /. Let I'’~* be the deter- 
minant obtained by replacing the jth column of 
I by the column of constants Ye, Yea, ***Yia; ***; 
and let I'’-* be the determinant obtained by 
replacing the jth column of [ by the column of 
quantities Dia YiaXa, Dw Y2aXay ***; Da YiaXa, *** 
It is evident that 


"=x!" (2-7) 


The Cramer’s rule solution of Eqs. (2-5) is, 
in this notation, 


= = T/T. (2-8) 


Proof that the determinant [ cannot vanish will 
be given in Section III. To make more explicit 
the dependence of the x; on the fixed-point 
coordinates, one can write 


Xoi= Dia CiaXay (2-9a) 


and similarly, 


(2.9b) 
(2-9c) 


Vou= i. CiaVa- 


201 = ; Ciara; 
where 


(2-10) 


Thus the xo;, Yo:, 20¢ are linear functions of the 
Xa, Va, Za, respectively, with coefficients cia which 
depend only on the N,,: a change in the position 
of any fixed point of a Gaussian network will 
produce a proportional change in the mean 
position of any other element of the network. 
It should be noted that the same constants 
appear in Eqs. (2-3a), (2-3b), and (2-3c) and 
hence in Eqs. (2-9a), (2-:9b), and (2-9c). It 
follows that a shift of a fixed point of a Gaussian 
network in any direction will produce a shift 
of the mean position of any other element of the 
network, in the same direction, and with mag- 
nitude the same for all directions of shift. 

The coordinates of the junctions relative to 
their mean positions are 


X 5 =Xi—Xoi, Y;=Yi—Yoi, 


for the fixed points, X2= Ya=Z_2=0. In terms of 
the variables xo;, Yor, 2,, and X,, Y,, Z,, the 
primary configuration function C becomes, by 
use of Eqs. (2-3), 


C=K exp{ —(3/2P) OY (1/Nz»)[(xor — x0»)? + (Yor — Yor)? + (Sor —Zor)* J} 


™>yv 


Xexp{ — (3/27?) 2 > (1/N.)[(X-—-X,)?+(¥,— Y,)?+(Z,—Z,)* ]}. 


t>v 


(2-12) 
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The first exponential is a constant determined by 
the positions of the fixed points; the second ex- 
ponential depends only on the coordinates of the 
junctions relative to the points f;, and is, aside 
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from a normalizing constant, just the proba- 
bility density function for these relative coor- 
dinates. A useful alternative form for C employs 
the constants y,, of Eq. (2-4): 


C=K exp | aa (3/2/?) > > 2 YrvLXorXov + YorVorv +20r20y |} 


Xexp{ — (3/2?) Vdls vil XiX7+ ViVj;4+Z,Z;]}. 


As is indicated by the subscripts, the sums in the 
first exponential are sums over all fixed points 
and junctions, whereas the sums in the second 
exponential are over the junctions only. 

In terms of the variables x, y, 2, we have cor- 
respondingly 


C=K exp{ —(3/22) © doy vax +r +212r]}. 
(2-14) 


III. INTEGRATION OF THE PRIMARY 
CONFIGURATION FUNCTION 


We must now consider the problem of inte- 
grating C over some, but not all, of the variables 
X;, Yi, Z; or Xi, yi, 2:. The integrations over the 
X’s, the Y’s, and the Z’s can be considered 


+00 
fax. 





+0 +2 n n 


dX::: dX, exp| 


—o —D 


= [(*!)/ {Dag }*] exp| 


provided only that D;®>0 for i<g. Otherwise 
the integral does not converge. 

Proof of this theorem involves partial rear- 
rangement of the quadratic form in the ex- 
ponent of the integrand into a sum of squares, 


with repeated use of the relation 

(p+1) 
Di; D pp” — Dp Dp =D,°mD), y-1- (3°3) 
It will be noted that the determinants on the 
left of Eq. (3-3) are (except for possible inter- 
change of rows and columns) first minors of 
D;;‘°+» obtained by striking out the rows p or 1 


and the columns ? or 7; oe is the second 
minor obtained by striking out rows p and 1, 
columns p and j. The truth of Eq. (3-3) then 
follows from a well-known relation of a deter- 
minant to its adjoint.’ 


3M. Bocher, Introduction to Higher Algebra (The Mac- 
millan Company, New York, 1929), p. 33. 


(2-13) 





separately and offer identical problems. They 
are easily accomplished by the use of the fol- 
lowing basic formula: 

Theorem I. Let the constants y;; form a sym- 
metrical matrix T. Let D®=1, and let D;;‘” 
denote the p-rowed minor of |I'|, obtained by 
striking out all rows but the first to the (p—1)st 
and the 7th, and all columns but the first to the 
(p—1)st and the jth: 


Yi Vij 


Y21 


Y12 * Jug 


Yp-1, 3 
Vii 


oe *Yp-1, p-l 
. Vi, p—l 


Yp—1,1 Yp—1,2 
Vil 7 i2 


Then 


—2 2 WwXiX,| 
ij 


—[1/Due IEE DsleXXj}, (3-2) 


i>q j>a 





Rearrangement of the quadratic form proceeds 
as follows: On completing the square of the 
terms in X,, one finds 
ei es VX X5= ul Xit(1/yn) DV Wj:X 5}? 

j>1 


+ DVlvis— (menu) |XX; 


i>l j>l 


= (Dy©/D) Xi¢+> (Di; /Du™)X; : 


j>l 


+X LU (D2, /Du)X :X j. 


i>1 j>l 


(3-4) 


On completion of the square of the terms in X:, 
the last sum on the right becomes 
LY U(Di;s/Du)X Xj 
i>1 j>l 
= (Dx /Di®){Xx+Z (Dos? Dn) X 3)" 
j>2 


HY LU (Dis/Do) X 5X 5, 


t>2 j>2 


(3-5) 





oroba- 
coor- 
iploys 


(3-2) 


oceeds 
of the 


(3-4) 


i n Xa, 
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when use is made of the fact (Eq. (3-3)) that 
(Di /Du™) — (Dai D2;5?/Du Da) 
=Dj;°/Dn™. (3-6) 


Continuation of this process leads to no dif- 
ficulty if all D,,“ #0; one finds 


EL wXiXy 


P . n 2 
=¥ (Du0/DE2s)| XE (Dy /Du)¥,| 


j>t 
HDL (Dis?t?/Dyp)X Xj, (3-7) 
i>p j>p 


if D;, 40 for <p. When p takes on its maxi- 


+a +a 
r= f dX dXq:-- 


i=l 





—O —oO 


mum value, , the second sum on the right 
disappears, and one obtains 


p Dd VisX Xj 
. # 


= L(Ds/DEs {XAT (Dy/D)X 5h , 


j>i 


(3-8) 


valid if all D;“° 40. The possibility of such a 
reduction of a quadratic form to a sum of squares 
and the coefficient of the last term in X, are 
well known,’ but the full form of the reduction 
is not easily found in the literature.’ Let us now 
consider the integral 


+n 


dX, exp | _ >: po VijX X 5} (3 -9) 


over the first g coordinates X;. By Eq. (3-7), with p=q, this becomes 


I=exp| —(1/D.) == DsexXX;1 


i>@ j>q 


+ +2 +00 q : 2 
x f dX,| 4X;--- dX, II exp| — (Dui DSN] XA LDu'/DWXy] |, (3-10) 


—RQ -—® —o i=l 


if D;“° 40 for «<q. The integration over X; can 
be carried out at once, since this variable appears 
in only one exponential factor. If Du™>0, the 
integral of this factor over X; converges: 


+0 2 
f dX, exp = Dx(%| K+ Z(Dy/Du”)X, 


j>l 
= {¢/Dy™ } } 


After this integration the integrand involves X2 
only through a single factor; integration over X2 
replaces this factor by the constant factor 


(3-11) 





I=[rt"/{ Doo} TIT exp} - (Dsio/ DN] Xe LDul0/Du)X] ; 


j>i 


valid if D;“°>0 when i<q and D,;,;“ +0 when 
i>q. 

The possibility of applying Theorem I in all 
cases in which we are interested is assured by 


Theorem II. The constants y;; of a coherent 
network form a symmetric matrix such that |T| 
and all its principal minors have positive values. 

The truth of this theorem follows easily from 


j>i 





{rDy™/Doo}3, if Do >0. Integration over X3 
then replaces an exponential factor by the con- 
stant {D2 /D3;}?, and so on. Integration 
over all g variables yields 


I= [Lwt!?/{ Dag }*] 

xexp| —(1/Dae) LX Dse0x.X;, (3-12) 
i>q j>q 

if D;;® >0 when i <q; otherwise the integral does 

not exist. This completes the proof of Theorem I. 

Similar use of Eq. (3-8) yields 


(3-13) 





the theorem:* The principal minors of a deter- 
minant |A| of a generalized quadratic form that 


4M. Bocher, Introduction to Higher Algebra (The Mac- 
millan Company, New York, 1929), p. 140. 

5 The author is indebted to Professor A. E. Ross for 
calling to his attention that this is a known result: See, for 
instance, W. Tartakowsky, Bull. Akad. Sci. U.R.S.S. 7, 
111 (1929). He wishes to acknowledge other helpful dis- 
cussions with Professor Ross on the material of this section. 

6 See, for instance, H. W. Turnbull and A. C. Aitken, 
An Introduction to the Theory of Canonical Matrices 
(Blackie and Sons, Ltd., London, 1945), p. 91. 
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is positive definite are all positive. The quadratic 
form 
i» » VewXeXyt KX" 

na + ee (1/Ny) (x; —X»)?-+« La Xa" 


jy 


(3-14) 


is positive definite, if the y’s are the constants 
defined for a coherent network, x. is the coor- 
dinate of fixed point a, and x is a positive con- 
stant. It is obvious that this quantity can vanish 
only if x.=0. For every junction or fixed point 
in the network directly connected to some fixed 
point a, there will be a term (1/Nya)(x;—Xa)? 
on the right of Eq. (3-14); this can vanish only 
if x,=xXa=0. Similarly, there will be a term 
(1/N,,)(x,—x,)? for every junction or fixed point 
v directly connected to one of the junctions or 
fixed points 7; vanishing of the quadratic form 
will require that for all of these x,=x,=0. By 
extension of this argument it follows that the 
quadratic form can vanish only if x,=0 for every 
junction or fixed point that can be reached from 





fan far... fav.ctr. fg, °°°3%1, To, «°° Tn) 
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some fixed point a by a path through the net- 
work. Since the network is coherent this includes 
every junction and fixed point; the quadratic 
form can vanish only if every variable in it 
vanishes, and it is thus a positive definite form. 
The corresponding determinant |A| is a deter- 
minant with a row and column for every junction 
and fixed point in the network; its elements are 
the y,, except that the diagonal element cor- 
responding to fixed point @ is Yaa+x. By the 
theorem, every principal minor of |A| is posi- 
tive. Since the determinant |I°|, with which we 
are concerned, is a principal minor of |A|, ob- 
tained by striking out the rows and columns 
corresponding to all fixed points, the truth of 
Theorem II follows immediately. 

These theorems can now be applied in the 
integration of the configuration function C over 
the coordinates of any chosen set of junctions of 
the network. If these junctions are numbered 1 
to qg, the integration of C, as given by Eq. (2-13), 
yields 


= K (2xl?/3)34[ Dag } exp { poi (3/21?) ; a V rv. XorXov + YorVou + Z0r20v_| } 


Xexp{ —(3/2?)[Da J" XY LX Dis LX X54 ViVj4+Z2Z;]}. 


To obtain the configuration function of any 
given junction of a Gaussian network one has 
only to assign this the highest number, 2, in 
numbering the junctions, and then to apply Eq. 
(3-15) with g=n—1. Similarly, to obtain a con- 
figuration function for an arbitrary pair of junc- 
tions one can number these as junctions »—1 
and n, and apply Eq. (3-15) with g=n—2; the 
other junctions can be numbered arbitrarily. 





C(r., ts, «+ +) =Ki exp{ —(3/2P)(1/|T |) Le Yig Vas XaXp+YaVetZaze |}, 


where 


(4-2) 


is the determinant |T'| bordered by the column 





C(r., Tg, °° ) =K, exp { i (3/2/?) +d (1/N2») [(x0r — Xor)? + (Yor — Mor)? + (Zor — Zor)? J} ° 


jv 


(3-15) 


i>q j>4@ 





IV. TOTAL NUMBER OF CONFIGURATIONS 
POSSIBLE FOR A GAUSSIAN NETWORK 

To determine @, the total number of con- 
figurations possible for a Gaussian network, as a 
function of the positions of the fixed points, it 
is necessary to integrate the primary con- 
figuration function C over all the coordinates 
X;, Vi, Z; or xi, Yi, 2. The most compact result is 
found by applying Theorem | in integrating C, as 
given by Eq. (2-14), over all x;, yi, 2;. This yields 


(4-1) 





of constants yi2 and the row of constants ‘g; 

A useful alternative form expresses this result 
in terms of the mean coordinates xo;, Yor, Zor of 
the junctions and fixed points. It is evident from 
Eq. (2-12) that integration over the coordinates 
Xi, Y;, Z; yields 


(4-3) 
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If the network is coherent the constant K, has 
the form shown in Eq. (3-15) with g=n. If the 
network is not coherent it can be divided into a 
coherent network and unattached material. The 
presence of the unattached material will increase 
K, by a factor equal to the number of con- 
figurations which it can take on, but will not 
change the exponential in Eq. (4-3). Further, 
for any pair of junctions 7 and v in the unattached 
material, one has xo; =Xov; Yor = You, Z0r = 20. Thus 
Eq. (4-3) has the same form for non-coherent 
as for coherent networks; the same is true for 
all other results given in this section. 

A case of special interest is that in which the 
system of fixed points can be subjected to a 
homogeneous deformation. For example, in the 
simple model for a cube of rubber subjected to 
unilateral stretch! the fixed points of the network 
are constrained to lie on the surface of a rec- 
tangular parallelopiped with edges of length 
L,, Ly, Lz in such a way that 


Xea=XeQLe, Ya=VaLy, Za=ZaL,. (4-4) 


By Eq. (2-9), the mean coordinates are then also 
proportional to L,, L,, L.: 


Moi=xOL,, 2,=2,L,, (4-5) 


Voi = y OL,, 
where 

x = as CiaXa, yi = Te as CiaVa, 

2; = Zs Cate™, 

are the mean coordinates of the ith junction 
when L,=L,=L,=1. The number of con- 
higurations possible for the network then depends 
only on L,, Ly, Lz; Eq. (4-3) becomes 
€(L,, L,, L.) 

= K, exp} —(K,L/7+K,L7+K.L,)}, 
where 


K, - (3 /21*) > ( 1 / Nyy) (x, = xF, 


ty 
K, - (3/2P)5 5 (1 / New) (vy, — y,)?, 
t>v 


K.= (3/2P) 0 (1/Nw) (2-5 —2,)?. 


t>v 


(4-6) 


(4-7) 


(4-8) 


The number of configurations possible for any 
coherent Gaussian network is thus a Gaussian 
function of the parameters L,, L,, L, which deter- 
mine the external dimensions of the system. 
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If the network is isotropic on the average, 
then K,=K,=K,; one can write 
€(L., L,, Lz) 
= K, exp| —K.(L+L,7+L)}, (4-9) 
with 


Ke = (1/2?) (1/Nw) C(x 


tT>v 


“+ (y, — y,)2+ (z, —z,)?], 


—x,)? 


(4-10) 


A variant form of Eq. (4-10) has been noted 
elsewhere.’ Let the vector extension of the 


rv-segment be denoted by 


E,, =i(x,—x,)+j(v;—y) +k(s,—2,). (4-11) 


It will be shown in Section V that the average of 
this vector over all configurations of the network 
is 
_ Nov) +j(vor = Yor) 

+k(2, —Z,). 


\ E.,)m = i(Xo, 


(4-12) 


The “‘vector-mean fractional extension”’ of this 
segment is the magnitude of this vector divided 
by the ‘‘total length’ N,,J of the segment. Let 
the value of this vector-mean fractional exten- 
sion, when L,=L,=L,=1, be 
hey _ [ (x, —x,)24 (vy, — y,)? 

+ (2, —s,)2 P/N,1. (4-13) 
Then one has simply 


K2=4 FD Norn P. 


T>v 


(4-14) 


The significance of this form of the result is dis- 
cussed elsewhere.* 

V. PROBABILITY DENSITY FUNCTION FOR A 

SINGLE JUNCTION 

It is easy to show, without entering into de- 
tailed calculations, that the probability density 
function for any element of a Gaussian network 
is a Gaussian function about a mean or most 
probable position : 


P(r) = {3/2rP?M;}? exp! — (3/2PR,) |r; —1Tox|?}. 
(5-1) 


7H. M. James and E. Guth, J. App. Phys. 15, 294 
(1944). 

5H. M. James and E. Guth, J. Chem. Phys. 15, 669 
(1947). 
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This function involves two parameters: f;, the 
mean position of the junction, and 9t;, which 
determines the breadth of the probability dis- 
tribution about this point. The constant Jt; is the 
number of links in a chain which has a con- 
figuration function [Eq. (1-1)] of the same 
form; in other words, the function P,(r;) is the 





Pa(tn) = { (3/22) (Dan /Dy1,»-1)}4 exp} 


Thus 


(n—1) 


Mn = Dp, a—1/D,,‘” = r.*/T. (5 ' 3) 


Since any junction can be taken as the nth, it 
follows that 


N, =T,/T. (5-4) 


(It will be observed that only the notation in- 
troduced after Eq. (2-6) is suited for this 
generalization. ) 

Equation (5-2) shows that the mean and most 
probable values of X,, Y,, and Z, are zero. It 
follows that the mean and most probable values 
of x:, yi, and 2; are the quantities %o;, Yoi, 20: 
defined by Eq. (2-1) and discussed in previous 
sections. Thus in Eq. (5-1), fo;=(xo:, Yo, Z0:)- 
From this the truth of Eq. (4-12) follows without 
difficulty. 

The parameter Jt; depends on the V,,—on the 
structure of the network—but not at all on the 
positions of the fixed points. When a fixed point 
of a Gaussian network is shifted, the mean posi- 
tion of each junction undergoes a proportionate 
shift, but there is no change in the fluctuation 





Pei n—l (Xn —1, Xn) 





— (3/2P) (Dan /Da—r, na) (Xn? + Vn? +Zn*)}. 
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same as the probability density function for the 
positions of one end of a chain of Yt, links, when 
the other end of the chain is fixed at the point r»,. 
We can now determine these quantities in terms 
of the quantities previously discussed. 

It follows at once from Eqs. (1-8) and (3-15) 
that 


(n—1) ~ 
(5-2) 


of the junction coordinates about their mean 
values. 

It should be noted that Eqs. (2-1) will deter- 
mine the same values of x;, y;, 2; for 77, whether 
one treats x;, y;, 2; as variables of the problem or 
assigns to them the fixed values x0;, Yoi, 20:. In 
other words, one can fix any chosen junction of a 
Gaussian network at its mean position without 
thereby altering the mean position of any other 
junction. On the other hand, this conversion of 
a junction into a fixed point will change the 
values of some of the Yt;; it will reduce the fluc- 
tuation of the coordinates of at least some of the 
other elements of the network. 





VI. PROBABILITY DENSITY FUNCTION FOR 
A PAIR OF JUNCTIONS 


The probability density function for the 
(n—1)st and nth junctions of a Gaussian network 
is easily found by application of Eq. (3-15). It 
factors into three parts of identical form, de- 
pending on the X’s, the Y’s, and the Z’s, re- 
spectively. Each of these factors is itself a 
probability density function for the coordinates 
involved. The X dependent function is 


(n—2) (n—1) 


= (3/2ml?){ Dan /Dyo oo} # exp { —(3/22)(1/Dy-2 ns)[Dona naXn— 







(n—1) 
It should be noted that Dn n—1 corresponds to 


—T,’, not to +T,’, and that 
r/>0 (6-3) 


for all 7 and 7. If the two junctions under con- 


In the ['-notation this becomes, for an arbitrary pair of junctions, 


Pi (Xi, Xj) = (3/2)eP {T/T ;7}) exp{ —(3/2P)(1/T i) (T PX? — 20 PX Xj + THX 7}. 





(n—1) 


+ 2D, naXnXaatDan X97}. 


(6-1) 


(6-2) 


sideration are numbered as the (n—1)st and 





nth, one obtains rae? by striking out the mth 
row and (n—1)st column of €. If these junctions 
are numbered as the ith and jth, with a cor- 
responding increase in the numbering of other 
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junctions, one obtains (—1)*’T' on striking out 
the ith row and jth row column of I’. These two 
determinants differ, however, by n—j trans- 
positions of rows and »—1—1 transpositions of 
columns, if 1 <j, and by »—j—1 transpositions 
of rows and n—i transpositions of columns, if 
i>j. In any case they differ by a factor of 


(—1)i+41; ited and I';/ accordingly differ by a 
factor of (—1). Now, all cross terms of the 
primary configuration function C have coeffi- 
cients which are zero or negative, and it is easily 
shown that this remains true when C is partially 
rearranged into a sum of squares, as in Eq. (3-7). 


: (n—2) acs 
Since D,~2,n-2 is positive, for any coherent net- 


work, it follows from Eq. (3-7) that Do <0, 
no matter how the junctions are ordered; the 
truth of Eq. (6-3) follows at once. 

If the 2th and jth junctions are not mechani- 
cally connected—if there are no paths between 
them in the network which do not pass through 
a fixed point—the coordinates of the two junc- 
tions will be statistically independent. There will 
then be no term in X;X; in Eq. (6-2), and 
P;(X;, Xj) will factor into independent prob- 
ability density functions for X; and Xj, in which 
the coefficients in the. exponential will be 
—3/2PM; and —3/2PN;, respectively, by Eq. 
(5-1). To exhibit the effect of mechanical 
coupling between the ith and jth junctions it is 
sometimes convenient to write Eq. (6-2) as 





Pi(w, Xj) =(3/2nP){ T/T ;;*}3 


Pi(Xi, X 3) =K’ exp{ —(3/2P)[(a/M)X? 
+ (b/MN)X PF + (c/(MMNj))X Xj J}. 


The constant ¢ is necessarily zero or positive 
[see Eq. (6-3) ]. When the junctions are not 
coupled, c=0, a=b=1. It is therefore possible 
to take c as a measure of the strength of the 
coupling between junctions of given J; and N;. 
Comparing Eq. (6-4) with Eqs. (5-3) and (6-2), 
one finds 


(6-4) 


a=b=P/T/TT,;', 
c=2Pe{ryery}i/rr,;. 
By the use of 
rerg—(P)2=CT yj", 


(6-5) 


(6-6) 


a relation which can be proved in the same way 
as Eq. (3-3), it is easily shown that 


c?=4a(a—1). (6:7) 


VII. FLUCTUATION OF RELATIVE COOR- 
DINATES OF TWO JUNCTIONS 


To obtain the probability density function for 
the relative coordinate §=x;—x; of the ith and 
jth junctions, one can integrate P;;(X;, Xj) over 
X;, say, while keeping & constant [see Eq. 
(1-15) ]. Let 

£o = X0i — X05, (7-1) 
and 


w=t—f=X;—X;j. (7-2) 


Equation (6-2) can be rewritten as 


Xexp{ — (3/2P)[(Vi— 20 e+ 0 )/Vy WX s+ LPP —PA)/PA—2VI+T) OP) 


with the aid of Eq. (6-6). Integration over X; 
vields 


Ri(&) = {(3/2e2)(V/[Pé—2re+rj)})}3 


Xexp { — (3/27) (T/(T#—2ry+Ty) ][é—& }*}. 
(7.4) 


The probability density function for the relative 
coordinate £ is thus a Gaussian function about a 
mean value £. 

The probability density function R;;(£, 7, ¢) is 
the product of three formally identical factors 
for the three variables. It can be written, in 


Xexp{ — (3/22)(P/Pi—2P +0 y)w}, (7-3) 





analogy to Eq. (5-1) as 
R;;(é, n, ¢) = {3/2eP?M,;|3 exp} —_ (3/2PN:;) 


x [(E— £0)? + (n— m0)? + (6 — fo)" J}, (7-5) 


where 


Nij= (1 i—2PZ+T,/)/T. (7-6) 


This important result implies that im a thermally 
agitated Gaussian network the relative coordinates 
of any pair of junctions will fluctuate according to 
a Gaussian distribution about mean values &, no, 
¢o. The mean values in question will depend on 
the coordinates of the fixed points, but the 
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breadth of the distribution, specified by 9t;;, will 
depend only on the structure of the network. 
It will be noted that 


NM 5;=N4AMN;— 27/7, (7-7) 
or, by Eqs. (6-5) and (6-7), 
Ni p=N AM, — [2c/ (14+ (A +e*)4} JOM)? 
>(MNZA—NAP. (7-8) 


For two statistically independent junctions 


the breadth of the distribution of values for 
t, n, ¢ is greater than that for the coordinates of 
either junction. Any coupling between the junc- 
tions reduces Yt;; below this maximum value, 
thus reducing the fluctuations in. the values of 
t, n, ¢. There is a lower limit of Jt,;; consistent 
with given values of Jt; and Yt; Very strong 
coupling of the junctions will, however, make Yt; 
approach ¥i;, and reduce this lower limit to zero; 
the fluctuations of values of &, 7, ¢ in a Gaussian 
network can be arbitrarily small. 


(7-9) 


VIII. ELASTIC THREAD MODEL EQIVALENT 
TO A GAUSSIAN NETWORK 


The preceding sections have given a formally 
complete treatment of the statistical properties 
of a Gaussian network; the probability density 
functions have all been expressed (or can very 
easily be) as functions of the fixed-point coor- 
dinates and certain determinants involving the 
N,,. The explicit evaluation of these deter- 
minants may, however, offer serious difficulties. 
It is, therefore, very useful to have an alternative 
method of evaluating the statistical parameters, 
one that is better adapted to the detailed treat- 
ment of special problems. This method is based 
on the consideration of a mechanical system 
equivalent to the given Gaussian network—a 
network of elastic threads. Each of the quan- 
tities X0;, Cja, Ji, Mi; which describe the statistical 
properties of the Gaussian network has its 
analog in the equivalent thread network: the 
Xo; and Cj are the same, Yt; is related to the dis- 
tortion produced by unit force applied to the ith 
junction of the thread network, and Yv;; to the 
distortion produced by equal and opposite forces 
applied to the 7th and jth junctions. This reduc- 
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tion of the statistical problem to an equivalent 
mechanical problem makes applicable to its solu- 
tion a quite different set of mathematical tech- 
niques, as will be illustrated in Section IX. 

A Gaussian network can be defined by giving 
the equivalent number of links N,, in the segment 
connecting the rth and vth junctions or fixed 
points, for all pairs of indices. The equivalent 
thread network’ has the same connectivity, with 
the zv-segment of flexible chain replaced by a 
classical elastic thread for which the stress- 
strain relation is given by 


F=(3kT/PN,,)d, (8-1) 


F being the tension in the thread and d the total 
distance between its ends. In the configuration 
of static equilibrium of the thread network the 
rv-segment should run straight from the rth to 
the vth junction; that is, the threads are not to 
be looped over each other in a way incom- 
patible with this condition. (This corresponds to 
neglect of steric hindrances in the Gaussian net- 
work. ) 

The potential energy of the thread network, 
as a function of the position of the fixed points 
and junctions, is 


-+)=(3kT/2P) VD (1/N) 


t>yv 
x L(x, —x,)?+ (y; —4¥,)*+ (2,— Zy)* ]. 


There is thus a relation between the number of 
configurations C accessible to the Gaussian 
network and the potential energy V; of the 
equivalent thread network: 


(X¢.. +9, vee) 
=K exp{|—Vii(ta, °° 


Vi(Ta, Tp, °° Ti, Tj, ° 


(8-2) 


‘Ti? -)/kT}. (8-3) 


The positions of the junctions which make C a 
maximum for the Gaussian network are those 
which make the potential energy of the thread 
network a minimum. It follows that when the 
thread network is in equilibrium its junctions 
will lie at the points (xo:, ¥Yo:, 20;)—the mean 
positions of the junctions of the corresponding 
network of Gaussian flexible chains. All the 
previously given relations for the xo, thus apply 
to the thread network; conversely, any method 
by which one can determine the equilibrium 


* The usefulness of this model has already been noted 
in reference 1. 
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configurations of the thread network will supply 
information about the Gaussian network. 

A constant external force applied to a junction 
of a Gaussian network will change the mean 
positions of this and other junctions of the 
network. It is easy to show that an equal force 
applied to the corresponding junction of the 
thread network will produce the same shift in 
the equilibrium positions of each junction. In 
either case an external force F; in the x direction, 
acting on junction 2, introduces into the problem 
a potential energy which may be written as 


V.(x;) = — Fix;. (8-4) 


The potential energy of the thread network 
becomes V;+ V,. In the presence of this external 
force the probability of occurrence of a given 
configuration of the junctions of the Gaussian 
network is not simply proportional to the con- 
figuration function C; it is proportional to 


C(fa, -**Ti, +++) exp} — V.(x,)/RT} 


=K exp{|—(V.t+Vi)/RT}. (8-5) 


The most probable positions of the junctions of 
the Gaussian network (still the mean positions) 
maximize the quantity on the right; the equi- 
librium positions of the junctions of the thread 
network are those which minimize V;+ V,. Ob- 
viously, then, these remain the same. 

When the force F; acts on the ith junction of 
the Gaussian network, the probability density 
function for r; is the probability density already 
derived [Eq. (5-1) ], multiplied by the weight 
factor exp[ — V./RkT] and renormalized: 


Pi'(r,) 
=K exp{ — (3/2PM,) | r;—1oi |? + Fix:/kT}. (8-6) 
On rearrangement, this becomes 
Pi'(r;) = K’ exp{ — (3/27?M,) 
« C(x; — x01 — (PRP /3RT))? 


+(vi:—yoi)?+(s:—20:)? ]}. (8-7) 


The force F; thus changes the mean value of x; 
by an amount equal to PR; F;/3kT; by the result 
of the preceding paragraph this is just equal to 
the displacement X;; of the ith junction of the 
thread network due to a force F; applied to that 
junction. That is, 


X ;:=PRF;{/3kT, (8-8) 
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whence 


Mi= (3kT/P)(X;i/ Fi). (8-9) 


Thus 9; for the Gaussian network is proportional 
to the displacement X;; produced by unit force 
applied to the ith junction of the corresponding 
thread network. 

The %;; can be determined similarly. Let equal 
and opposite forces in the x direction, F; and 
F;=—F;, be applied to junctions 7 and j, re- 
spectively, of corresponding thread and Gaussian 
networks. Since, with both types of network, 
displacements due to different forces are additive, 
these paired forces will produce equal changes 
Aé in the relative coordinates of the thread- 
network junctions and in the mean-relative 
coordinates of the Gaussian network junctions. 
In the case of the Gaussian network the dis- 
tribution of relative coordinates £, n, ¢ is that 
already derived [Eq. (7-5) ], multiplied by the 
weight factor 


e VelkT — pF i(xi—zj) [kT — pF ik/kT (8 : 10) 
and renormalized. A derivation entirely parallel 
to that of Eq. (8-9) leads to the result that 


MNis= (3kT/P)(Abo/ Fi); (8-11) 


N,; for the Gaussian network is proportional to 
the change in the relative coordinates of junc- 
tions 7 and j of the thread network produced by 
oppositely directed unit forces applied to these 
junctions. If we denote by X;; the displacement 
of the 7th junction produced by a force F; applied 
to the jth junction, and so on, we can write also 


Finally, it should be noted that the Gaussian 
and thread networks exert equal forces on each 
fixed point. If the junctions of the thread network 
are allowed to take on their equilibrium posi- 
tions, the potential energy becomes 


+++) =(3kT/2P) > (1/N2) 


t>y 


x [ (Xor = Xo)” + (Yor a Yor)” + (20; 


Vo(Ta, Ta, 


aes Zor)? |. 


This is a function of the fixed-point coordinates 
only, since these determine all the xo,, Yor, Zor} 
it is related to the total number of configurations 


(8-13) 
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possible for the Gaussian network [Eq. (4-3) ]: 
C(ra, ts, «+ +) = Kye-V/*?. (8-14) 


Now the mean x-component of force exerted on 
fixed point a by the thermally agitated Gaussian 
network is, by the usual methods of statistical 
mechanics, 


kT(0/0%q) INn@ =kT(8/AxXa)(— Vo/kT) 


=—9Vo/dxe. (8-15) 


This is, of course, also the x-component of force 
exerted on fixed point a by the thread network 
in its equilibrium configuration. 


IX. REGULAR CUBIC NETWORK 


The regular cubic network is the simplest 
typical three-dimensional network. It will be 
discussed here both for its intrinsic interest!’ and 
as an illustration of the power of the methods 
described in the preceding section. 

A regular cubic network will be defined to 
consist of segments all having the same number 
of links, connected together like the edges of a 
set of equal cubes filling a cubical region on the 
surface of which the fixed points of the network 
are regularly spaced. Let the network segments 
consist of M links each, and let there be N of 
these segments in each chain extending between 
opposite faces of the bounding cube. There will 
be in all 3(N—1)? of these chains of MN links, 


ater \ 
Cescee nam am ene «2 an amen an am as aman Ge nee am on 





Fic. 2. Regular cubic network, N =3. Solid lines—Gaussian 
chains; X—fixed points; O—junctions. 


10 The results are employed by James and Guth, refer- 
ence 8. 
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extending between fixed points 

[0, (m/N)L, (n/N)L and (L, (m/N)L, (n/N)L], 
or 

[(m/N)L, 0, (n/N)L and ((m/N)L, L, (n/N)L), 
or 

[(m/N)L, (n/N)L,0)and [(m/N)L, (n/N)L, L), 


for all integral values of m and n from 1 to N—1. 
Such a network is illustrated in Fig. 2, for the 
case VN =3. 

The junctions of any cubic network are con- 
veniently specified by three indices p, g, r, which 
give the number of segments lying between this 
junction and the planes x=0, y=0, and z=0, 
respectively. Let N(r, s, t; u, v, w) be the number 
of links in the segment between the junctions 
(r, s, t) and (u,v, w). Then in the regular cubic 
network 

N(r, s, t; u,v, w)=M, 
if w=r+1, v=st1, w=t+1, (9-1) 
N(r, s, t; u, v, w) = ©, otherwise. 


Equation (9-1) defines the structure of the net- 
work; its state of extension is specified by giving 
the coordinates of the fixed points, 


x(r, s,t)=(r/N)L, 
y(r, s, t) =(s/N)L, (9-2) 
2(r, s, t)=(t/N)L, | 


where one of the three indices 7, s, ¢, must be 0 
or N. 

By Eq. (8-1), the equivalent elastic-thread 
network consists of a cubic array of threads for 
each of which the stress-strain relation is given by 


F=(3kT/MP)d =kxd. (9-3) 


The nature of the equilibrium configuration of 
this model is obvious: 


xo(r, s, t) =(r/N)L, 


yolr, 5, t) = (s/N)L, 
ao(r, s, t)=(t/N)L, 


(9-4) 


for all values of the indices. These are then the 
mean positions of the junctions of the Gaussian 
network. 

The main problem to be considered is the cal- 
culation of N(r, s, t) and N(r, s, t; u, v, w), which 
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determine respectively the fluctuations of the 
coordinates of the junction (r, s, t) of the Gaus- 
sian network and the fluctuations of the relative 
coordinates of junctions (r,s,¢) and (u,v, w) 
about their mean values. The determinants of 
the general theory are not easily calculated; we 
shall therefore apply the methods described in 
Section VIII. Specifically, we shall consider the 
effect of applying a unit force in the x direction 
to the junction (u,v, w) of the elastic thread 
network. This will produce a displacement 
X(u,v,w;u,v,w) of this junction from its 
previous equilibrium position, and a displacement 
X(r, s,#; u,v, w) of the junction (r, s, t). These 
displacements can be calculated by a method 
closely related to the normal-mode calculations 
for cubic crystals; the Jt’s can then be obtained 
by application of Eqs. (8-9) and (8-12). 

Let any system of external forces F(r, s, t), all 
having the x direction, be applied to the junc- 
tions of the thread network. They will produce a 
system of junction displacements X(r, s, ¢), all 
having the x direction. The x component of the 
force exerted on junction (7, s,¢) by the thread 
running to junction (r+1, s, f) will be increased 
by 


AF=x{X(r+1, s,t)—X(r,5,t)}, (9-5) 


and similarly for the other threads entering this 
junction. The conditions for the static equi- 
librium of all junctions under the applied forces 
are thus 


F(r, s, t) =x{6X(r, s, £) 
—X(r+1, s, t)-—X(r—-1, s, t) 
—X(r, s+1, t)—X(r, s—1, t) 
—X(r,s,t+1)—X(r, s,t-—1)}, 
r,s, t=1,2, ---N—1. 


(9-6) 


For the fixed points, 
X(r, s, t)=0, if r or s or t=0 or N. 


(9-7) 


The quantities X(r, s, t; u,v, w) in which we are 
particularly interested solve these equations for 





X(r, s, t3 u, v, w) 


2 sin(rau/N) sin(xar/N) sin(xbv/N) sin(xbs/N) sin(rcw/N) sin(rct/N) 
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F(u, v, w) =1, 
(9-8) 
F(r, s, t)=0 unless r=u, s=v, t=w. 
It is easy to solve Eqs. (9-6) and (9-7) for a 
normal-mode distribution of forces given by 


Fase(r, Ss, t) =sin(rar/N) sin(rbs/N) 


Xsin(rct/N). (9-9) 


The indices a, 6, c which distinguish the normal 
modes can each run from 1 to N—1. It is natural 
to expect that the normal-mode forces will 
produce proportional displacements of each 
junction: 


X abel’, S, t) =have sin(xar/N) 


<sin(rbs/N) sin(xct/N). (9-10) 


The truth of this conjecture is easily proved by 
introduction of Eqs. (9-9) and (9-10) into Eq. 
(9-6); one finds that this is then satisfied for all 
r, s, t, provided only that 


have = {4x[sin?(xa/2N)+sin2(rb/2N) 


+sin*(rc/2N) ]}-'. (9-11) 


Since Eqs. (9-6), (9-7) form a linear system, 
any sum of solutions will be a solution. In par- 
ticular, any sum of normal-mode distributions, 

F(r, s, t) = > gaveFare(r, 5, t), (9-12) 
a, be 
will give rise to the displacements 
X(r, 5, t) — 2. LabeX abel’, 5, t). 


a, b,c 


(9-13) 


To solve the problem at hand it is only necessary 
to find constant gave such that the F(r, s, t) are 
those specified by Eq. (9-8); the corresponding 
X (r,s, t) are then given by Eq. (9-13). 

The constants gare can be determined from 
Eqs. (9-8), (9-9), (9-12) by the usual method of 
Fourier analysis: 


Lar. = (2/N)* sin(xau/N) 


Xsin(rbv/N) sin(rcw/N). (9-14) 


Thus 





xN3 a, b,c 


By Eq. (8-9), one has then 


(9-15) 


sin?(xa/2N)+sin?(rb/2N) +sin?(rc/2N) 


2M sin*(zar/N)+sin*(rbs/N)+sin?(xct/N) 





Nir, s, H= 


N3 z sin?(ra/2N)+sin*(rb/2N) +sin?(rc/2N) 
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Similarly, Eq. (8-12) yields 
Nir, s, t; u,v, w) 


2M (sin(rar/N) sin(xbs/N) sin(rcb/N) —sin(xau/N) sin(rbv/N) sin(xcw/ N) } 
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(9-17) 














NP abc sin?(za/2N)+sin?(rb/2N) +sin?(rc/2N) 


To make clear the physical nature of these results it is necessary to evaluate the sums in Eqs. 
(9-16) and (9-17). In Eq. (9-16) it is convenient to replace sin?’rar/N by 
similarly for the other terms in the numerator. Let 


cos(2rra/N) cos(2rsb/N) cos(2rtc/N) 


+(1—cos2mar/N), and 





1 
I(r,s,)}=— > 























- 


r— sto it~ 
er eas 
2 2 2 





’ ’ 


2 2 2 


It remains to evaluate the sums [(r, s, f). 


The sum /(0, 0, 0) requires separate attention. 
When JN is large—as will be assumed hence- 
forward—the summand in /(0,0,0) changes 
slowly with the indices a, b, c. The sum can then 
be approximated very closely by an integral: 


r/2 a/2 r/2 
7(0, 0, 0)=(2 ay f ax f dy dz 


X {sin’x+sin*?y+sin’z}— 


(2 =(2/x)*f vas dy{sin?x+sin*y}— 


X {1+sin’x+sin?y}— ° 


=(2/m fo dx} 1+sin?x}~! 


x K({1+sin’x}-!), (9-21) 


where K(_ ) is the complete elliptic integral. This 
integral has resisted evaluation in closed form; 
numerical integration shows that 

I(0, 0, 0)=1, (9-22) 
to within 1 percent. This may very well be the 
exact value. 


, (9-18) 
N3 ,.5,< sin?(ra/2N)+sin?(rb/2N) +sin?(rce/2N) 
One then has 
Mir, s, t) =(17/4) {1(0, 0, 0) — I(r, 0, 0) —I(0, s, 0) —I(0, 0, £) 
+TI(r,s,0)+J(r, 0, 1) +](0, s,))—T(r, s, t)}. (9-19) 
In a similar way one finds 
M(_(r-u s—v t- r+u s—v (—w 
Nir, s, tru, v, w) =N(r, s, O+N(a, v, w) -—-— (? —, — —*)- i(° : ‘) 
2 2 2" 2 2 





iy (> $S—v Bets rt+u stv (—w 
on Ca 9, ies paren 


rt+tu s—v tt+w r—u s+v t+w ci hea stv t+w 
AEE AEE SALE) 0m 


’ 


2 2 2 


= 


All the quantities J(r, s,?t) except J(0, 0, 0) 
are sums of terms with oscillating signs. They 
will all be smaller than /(0, 0, 0), and the smaller 
the more rapid the changes in sign—that is, the 
larger (7?+5?+7*)?, so long as 7, s, t<N/2. Thus, 
except for junctions very near one of the bound- 
ing surfaces, one can neglect all but the first 
term in Eq. (9-19), and write 


Nir, s, )=M/4. (9-23) 


To compute the corrections to Eq. (9-23), 
wherever they are appreciable, it will suffice to 
compute I(r, s,¢) for r/N, s/N, t/N<&1. The 
summand in I(r, s, t) then changes slowly enough 
to permit replacement of the sum by an integral: 


=f dof af dc 
c atens er ma cos24—c 
N N 
xX<- . (9-24) 
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The integrand oscillates sinusoidally with a 
variable amplitude 


(sin*(ra/2N)+sin?(rb/2N) +sin*(rc/2N) |. 


The integration over the regions in which a/N, 
b/N, or c/N are large, makes a very small con- 
tribution to the integral: the amplitude is there 
small and, more importantly, it changes so 
slowly that contributions of opposite sign to the 
integral 
nearly cancel each other. The value of the 
integral is thus determined mainly by the 
integration over the regions of smaller a/N, b/N, 


c/N, where one can write with good accuracy 
sin?(ra/2N)=(ra/2N)?, (9-25) 


and so on; even for the larger a/N this is not a 
very bad approximation. It is thus a good ap- 
oaeeus to write 


llr, s =f wf 7 
rN 


cos2 mi cos2 what cos2 —s 
J N 


» i l 
—. (9-26) 





a’+b?+c 


Extending the range of integration, we can also 
write 


1 +N +N +N 
I(r, s, n= f da f av f dc 
2n’?N J_n -N —N 


xp | (221/N)(ra+sb+tc)} 
Pt Aat/ NY rat b+t} 9.94 
e+P+e 





the sine parts of the complex exponentials con- 
tribute nothing to the integral when the range of 
integration is thus extended. 

To evaluate this integral it is convenient to 
use polar coordinates. The radial coordinate will 
of course be 


= (a+b +¢)'. (9-28) 


As the polar axis we choose the direction of the 
vector 
P=ri+sj+ik, (9-29) 


where i, j, k, are unit vectors in the a, b, and c 
directions, respectively. Then 


ra+sb+tc=(P-o) =(r?+5s?+£)'p cosé. 


(9-30) 


from neighboring regions will very 


We shall also replace. the integration over a 
cubical region by integration over a sphere of 
equal volume. This change in the range of 
integration is of little importance, since it 
affects only a region that contributes little to the 
integral—how little will be made evident im- 
mediately. We have then 


1 r Qn (6/n)4N 
I(r, s, D\=—_ f sinado f def d 
2r°N 0 0 0 . 


Xexp{(22i/N)(r?+s°+F)'p cos0}. (9-31) 


Integration over the angle coordinates is simple; 
on a change of variable in the remaining integral, 
the result becomes 


1 1 
I(r, s, )0=— 
7 (P?+s?+ +s°+P 2) 


(6/x)42x(r? +5? +-4)3 sinx 
x f dx——. 
0 x 


The asymptotic series for the integral on the 
right, for large values of the upper limit, is 





(9-32) 


f | dx(sinx/x) =2/2—cosx/x+---, 
0 


x>1. (9-33) 


Thus we have, to terms of the first order in the 
small quantity 1/22(7?+s?+?), 


I(r, s, )221/2e(r?+s?+F)3, r,s, tKN. (9-34) 


It will be noted that the approximation of 
neglecting higher order terms corresponds to 
extension of the integration of Eq. (9-31) over 
the whole of space; the smallness of the resultant 
change justified the change in limits which led to 
Eq. (9-31). 

Since Jt(r, s, f) is unchanged on replacing 7 by 
N-r, s by N-s, t by N-—t, it will suffice to 
consider 3t(r, s, 2) when all three indices are less 
than N/2. By use of Eq. (9-34) one can then 
write Eq. (9-19) as 


M ina ce. 
Nir, s, )=T\1- |-+-+-- 
4 2rLr s ¢t (r?+57)! 





(9-35) 





1 1 1 
_ (e+e) rarer d 


This approximation appears to be good for all 
junctions, for which r, s, and ¢ are all necessarily 
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1 or more. It does not apply to the fixed points 
(r, s, or t=0), for which the correct value is zero. 

There is an important difference between Eq. 
(9-35) for a three-dimensional network and the 
corresponding equation for a single chain. Let us 
consider a single chain consisting of NV segments, 
each containing M links. The value of Jt for the 
rth junction is easily computed as follows. The 
chain can be treated as a network with one 
junction, the rth, linked to two fixed points by 
chains with Nig=Mr and Nig=M(N-—r) links, 
respectively. The network determinant T has a 
single element: 


TP =yu=N/Mr(N—-r). (9-36) 

By Eq. (5-2), 
Nr) =1/yn=Mr(N—r)/N, (9-37) 
since T;/=D®=1. In the three-dimensional 


network all junctions except those very near the 
fixed points have essentially the same value of Yt, 
and fluctuate about their mean positions with 
_ the same freedom. In the single chain, however, 
9 varies much more strongly along the chain, 
being proportional to 7 over a large range; the 
junction coordinates have a mean-square vari- 
ation proportional to r(N—r). It thus appears 
that there is danger in extending conclusions 
reached by considering the properties of a single 
Gaussian chain to three-dimensional networks. 
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An expression for Jt(r, s, t; u, v, w) can now be 
written down by use of Eqs. (9-20), (9-34), and 
(9-35). We shall restrict our attention to pairs of 
junctions neither of which is close to a bounding 
surface—that is, to cases in which 7, s, t, u, v, w 
are all large but not close to N. Then 


Nir, s, HLMN(u, v, w)\=M/4. (9-38) 


Furthermore, each of the J’s in Eq. (9-20), except 
possibly J[(r—u)/2, (s—v)/2, (t—w)/2], has at 
least one large argument and is negligibly small. 
Thus we have 


Nir, s, t; u, v, w) 


M 1 1 


2| = C(r—u)?+(s—v)?+ (t—w)? FJ’ lias 





valid for junctions far from the bounding surfaces 
but possibly near to each other. The correct value 
of I(r, s, t; u,v, w) is of course 0 if r=u, s=z, 
t=w. Equation (9-39) seems to be very accurate 
in all other cases. The worst case is that of 
|r—u| =1, s=v, f=w, which requires the use of 
the approximate Eq. (9-34) for Z(3, 0, 0). It is, 
however, easy to show that the exact value of the 
integral J(4, 0, 0) as given by Eq. (9-24) is 4, in 
excellent agreement with the approximate value 
1/m given by Eq. (9-34). 
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The rigidity of rubber, considered as a network of 
flexible molecules with Gaussian configuration functions, 
can be calculated for a particular sample if one is given a 
complete description of the molecular network or certain 
types of statistical description. In particular, it is sufficient 
to know the distribution of lengths and vector-mean 
extensions of the segments of the network, or only the 
latter distribution if it is of Gaussian form. The assumption 
that the vector-mean extensions have a Gaussian distribu- 
tion corresponds to a similar postulate in the theory of 
Wall, and has been applied to a simplified network theory 
by Flory. In the complete theory it leads to calculation 
of the same proportionality between the rigidity of the 
material and the number G. of segments, per unit volume, 
in the “‘active’’ part of the molecular network. However, 
consideration of the process of cure shows that this 


assumption cannot be expected to be correct, though it 
does lead to results of the right order of magnitude. An 
alternative approach to the problem is based on a study 
of the increase in rigidity of the material as cure proceeds. 
It uses a more realistic picture of the process of cure than 
those hitherto employed, and uses Gaussian distribution 
functions only where they can be logically justified. The 
final result is of the same form as that given by the previous 
theory, except that Ga is replaced by Ba, the number of 
bonds formed within the molecular network during cure. 
Since B, is smaller than Gg by a factor of about two or 
more, the improved theory leads to the prediction of 
somewhat lower rigidities per bond formed. Available 
experimental results of Flory on Butyl rubber vulcanizates 
do not appear to give an adequate check on the quantita- 
tive features of the theory. 





N another paper on the network theory of 

rubber! we have shown that one can derive 
the form of the stress-strain relation for a unit 
cube of ideal soft rubber, 


on the following assumptions: 

(1) Soft rubber consists of a liquid-like mass 
of long, flexible molecules united into a network 
by a relatively few intermolecular bonds. 

(2) The molecular network is isotropic on the 
average. 

(3) Each segment of the network can take on 
a number of configurations which is a Gaussian 
function of the separation of its ends. 

(4) The steric hindrances between the mole- 
cules (as opposed to the constraints implied by 
their union into a network) can be neglected. 

(5) Compressibility and thermal expansion of 
the material are neglected. (More complicated 
forms arise when these factors are taken into 
account.) 

Equation (1) is valid whatever the structure 


‘H. M. James and E. Guth, J. Chem. Phys. 11, 455 
(1943), 


of the network. In this equation & is the Boltz- 
mann constant. The constant K, which deter- 
mines the rigidity of the material, will vary from 
sample to sample; to determine it theoretically 
one must know more about the network— 
whether there are few or many molecules in- 
volved in it, and so on. In particular, K for a 
piece of rubber will change as cure proceeds. It 
is the purpose of this paper to discuss the relation 
of K to the factors in terms of which one might 
specify the structure of a piece of rubber. It 
includes the calculation of K for a completely 
specified network, a discussion of the problem of 
determining K for a molecular network described 
only in a statistical sense, and the development 
of a third approach to the problem-.through a 
consideration of the process of cure in a more 
complete and more realistic way than has 
heretofore been attempted. 


I, CALCULATION OF K FOR A COMPLETELY 
SPECIFIED NETWORK 


A complete specification of the structure of a 
Gaussian molecular network will make possible 
calculation of the corresponding constant K. 
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670 H. M.. JAMES 
For the purposes of the network theory of 
rubber, a complete description of the molecular 
network will include specification of (a) the 
“configuration function’’ for each segment of 
the network, (b) the connectivity of the network 
into which these segments are united, and (c) the 
positions of any points of the network which are 
subject to external constraints—the ‘“‘fixed 
points” of the network. The configuration func- 
tion for a network segment gives the relative 


number of configurations possible for the seg- - 


ment, as a function of the separation of its ends. 
If this is of Gaussian form—an approximation 
which has been thoroughly discussed in reference 
1 and will be used throughout the present paper 
—it can be written as 






c(N,,, |t,—1,|) 


=k,, exp{ —(3/2PN,,)|r,—1,|7}. (2) 
Here r, and r, are vectors specifying the positions 
of the end points of the segments, and |r,—r,| 
is the distance between the end points. This is 
the form valid in the limit of small extensions 
for a chain of N,, independent links each of fixed 
length /; for other types of flexible chain it serves 
to define an effective number N,, of links in the 
segment. In the description of a network, / can 
be assigned the same value for every segment. 
If the network is considered as a system of 
junctions and fixed points joined by unbranching 
chain segments, its structure can then be speci- 
fied completely by giving the number of links 
N,, in the segment connecting each such pair of 
points. Specification of the positions of the fixed 
points then completes the description of the 
network. 

It has been shown elsewhere? how the number 
of configurations possible for a network of 
Gaussian segments can be expressed as a function 
of the N,, and the coordinates of the fixed points. 
Let the indices 7, j, refer to junctions, the 
indices a, 8, --- to fixed points, and the indices 
T, My V; to points of both types; >; then 
indicates a sum over all junctions, }}4 a sum 
over all fixed points, and >>, a sum over all 
points of both types. The number of configura- 
‘tions possible for the network can then be 


2H. M. James, J. Chem. Phys. 15, 651 (1947). 
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written as 


3 1 
€=K, exp} ——LE—L (wo, - 0») 
‘ ’ 2/? Ty N,v 


° 


a (Yo = Yor)” + (Zo; _ Zo»)? | : . (09) 


Here the double sum is over ail pairs of junctions 
and fixed points, or over all segments of the 
network; if there is no segment of the network 
joining the points 7 and », then 1/N,, is to be 
taken as zero. The quantities xo;, Yo:, Zo: are the 
mean values of the coordinates of the ith junc- 
tion, whereas Xoa, You, 20a are identical with the 
fixed coordinates Xa, Ya, Za Of the fixed point a. 
The xoi, Yo:, 0: can be expressed in terms of the 
Xa, Ya; Za as follows. Let 


Tre —1/N,,, THY 
Ce die > 1/N,, (4) 
vet 


and let T be the determinant of the coefficients 
ij, With one row and column for each of the 
junctions of the network: 


man Y12 rds Y in|} 
soe Ceeeeereoeece ieee | (5) 

enn 

1Yni Yn2 os Van | 


Let ['-* be the determinant obtained by re- 
placing the ith column of [ by the column of 
constants oa YieXa, doa Y2aXa, ***; let I and 
I‘ have the same form with y. and 2a, respec- 
tively, replacing xa. Then 


Pe /T. 


X0i = f% 
Yo = —T*9/P, (6) 
201i = a_i! je 


Equations (3) to (6) together give an explicit 
expression for @ in terms of the N,, and Xa, Ya; 2a.’ 
The constant K;, is of no interest. 

In deriving Eq. (1) we may consider a molec- 
ular network which fills a unit cube when the 
material is unstretched; the fixed points lie 
on the surface of this cube with coordinates 
Xa”, Ya, Za. To subject the material to a 





’ A more compact expression is derived in reference 2, 
but it is of less interest for, the developments of the present 


paper. 
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homogeneous deformation, the fixed points may 
be constrained to lie of the surface of a cuboid 
with edges L,, Ly, L, at points with coordinates 


Ka=XeML,, 

= 
Va = Va™ La, ( /) 
Sa=2_Ls. 


Then the Xo:, Yo:, 20: are also proportional to 
L,, Ly, L., respectively : 


Xo =xX,OL,, 
Voi = pS OL, (8) 


20: =2,L,., 


where aes 
x; = —| ~» (9) | # (9) 


and similarly for y; and z;. Then the number 
of configurations possible for the network be- 
comes a Gaussian function of the external 
coordinates only : 


€(L., L,, L.)=K, exp| —(K,L2 
+K,L,°+K,L,7)}, (10) 


where 


K,= (3/292)>> > (1/N,7) (x, —x,)?, (11) 
jy 
and K, and K, have similar forms. If the 
unstretched material is isotropic on the average 
one must have K,=K,=K,. Then Eq. (1) 
follows, with 


K=(1/E)DX(1/N,)[ (x — x)? 


Ty 
+ (vy, —y,)?+ (2, —2,)?]. (12) 


This, together with Eq. (9), gives the desired 
explicit form for K. 


II. CALCULATION OF K FOR A STATISTICALLY 
SPECIFIED NETWORK 


In practice, one can only hope to know certain 
statistical properties of physically occurring 
networks. Any one of several statistical specifi- 
cations of the network structure may serve to 
determine K; one of these is suggested by the 
results quoted in the preceding section. 

We define the ‘‘ vector-mean fractional exten- 
sion” of a network segment as follows. The 
instantaneous vector extension of the network 
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segment connecting the rth and rth junctions is 
(with arbitrary choice of sign) 


E, »y=I,—T, 
=i(x,—x,)+j(y,—y,)+K(2,—2,), (13) 


where i, j, k are the usual unit vectors in the 
coordinate directions. The mean value of this 
vector, for all possible configurations of the 
network, is 


(E, yay _ i(Xo, Saat Xov) +j(Vor — Yor) +k(zo, ime Zov) . (14) 


The vector-mean extension of the segment is the 
magnitude of this vector; its vector-mean frac- 
tional extension, A,,, is this magnitude divided 
by N,,/, the total length of the segment: 


Avy = [ (Xor —Xov)?+ (Yor — Yor)” 


+ (20, — Zor)? }}/N,,/. (15) 


In terms of vector-mean fractional extensions 
the number of configurations possible for the 
network becomes simply 


€=K, exp| ALIN? , (16) 


t>y 


Similarly, we have from Eq. (12) 
K= LDN [de FP, 
ty 


where 


he» = [ (x, —x,)?+ (y,© —y,)? 


+ (2, —2,)?}}/N,1 (18) 


is the vector-mean fractional extension of the 
ry-segment when the material is unstretched. 
The sum in Eq. (17) is over all segments in a 
unit volume of the network; each segment con- 
tributes to the rigidity of the material an amount 
proportional to the product of the number of its 
links and the square of its vector-mean fractional 
extension when the material is unstretched. This 
result is independent of the detailed structure 
of the network. 

It is thus evident that K can be determined 
without knowledge of the detailed structure of 
the network: it will suffice to know how many 
segments there are per unit volume of the 
network having any given number of links and 
any given value of A. Let F(N)dN be the 
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number of segments per unit volume for which 
the number of links lies between N and N+dN, 
and let G(N, A)dd be the fraction of these seg- 
ments which have a vector-mean fractional 
extension between \ and A+d\ when the ma- 
terial is unstretched. Then Eq. (17) becomes 


x 1 
k= f anf d\F(N)G(N, NX. (19) 
0 0 


This determines K in terms of two distribution 
functions, F(NV) and G(X, A); it will be shown 
later that under certain circumstances it is not 
even necessary to know the form of F(J). 

It is possible to consider as part of the network 
chain segments that are not actually linked into 
it. For such unattached segments all orientations 
of the vector E,, are equally probable, and 
(E,,)w and A,, are zero. It follows from Eqs. (17) 
and (19) that such unattached segments do not 
contribute to the rigidity of the material. The 
same is true of more complicated structures that 
are not attached to the network or to a fixed 
point; it is also true of all ‘‘loose ends’’—seg- 
ments and complex structures joined to the rest 
of the network at only one point. Only when a 
segment is subject to constraints that cause its 
ends to have different mean positions does it 
contribute to the rigidity of the material. It is 
easily seen that this requires that the segment 
lie on a non-retracing path through the network 
from one fixed point to another. The term 
‘active network”’ will be applied to the aggregate 
to these segments. , 

One can hope to determine K by establishing 
the forms of F(N) and G(N,A)—perhaps by 
reference to general principles, perhaps by con- 
sidering the process by which the molecular 
network is formed. One approach to the problem 
is the following. We may picture milled uncured 
rubber as a mass of unconnected Gaussian 
molecular chains which form a liquid-like mass, 
rather than a coherent network. In this liquid 
each chain segment can take on any orientation 
with equal probability; the vector-mean frac- 
tional extension of each segment is zero, and K 
is zero—the system may be very viscous, but it 
has no fixed form. Now let us consider the 
instantaneous distribution of extensions of the 
molecular segments in this liquid. The proba- 
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bility that a segment of N links has an extension 
r is proportional to the number of configurations 
of the segment consistent with that extension. 
The number of configurations of the ends of the 
segment consistent with separation r is propor- 
tional to r?; the number of configurations of the 
segment consistent with a given configuration of 
the ends is given by Eq. (2). Thus the probability 
that the segment will have an instantaneous 
extension 7, an instantaneous fractional extension 
Xi, is proportional to 


f(N, r) =P exp{ — (3/2PN)*} 
= N*Pd? exp{ —$NA?}. (20) 


It would be quite in accord with some current 
ideas to picture the process of cure as the 
instantaneous introduction of bonds which con- 
vert these originally unconstrained segments into 
a coherent network, and to suppose that the 
distribution of instantaneous fractional exten- 
sions of the segments at the moment of cure is 
the same as the distribution of vector-mean 
extensions in the completed network. In that 
case G(N,) for the network will differ from 
the expressions in Eq. (20) only by a normalizing 
factor : 


G(N, d) = (4/2) (3.N/2)) exp{ —3NX}. (21) 
Then 


k= aN-NF(N) J dd4/1? 
0 0 

X (3N/2)!M exp{—3NM}, (22) 
whence 


K= f dNF(N)=G. (23) 
0 


Thus K is found to be equal to G, the total 
number of segments per unit volume in the 
network, whatever may be the form of the 
distribution function F(V).4 

It should be emphasized at once that we do 
not consider this to be an acceptable solution of 
the problem. It is thus, however, that one must 


‘It should be noted that the integration over A 10 
Eq. (22) extends to infinity. The Gaussian distribution 
is only an approximation for any chain of finite length, 
and does not exclude values of \ greater than 1. Ine 
error resulting from the integration over these fictitious 
fractional extensions is completely negligible if N is large. 
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translate into terms of the complete network 
theory of rubber the theories of Wall® and of 
Kuhn which are not explicitly network theories, 
and a calculation of Flory and Rehner’ which 
employs a specialized and incompletely devel- 
oped network theory. Wall derives an equation 
of the form 


Z=GkT[L.—1/L2 ] (24) 


by a formalism that, as it stands, cannot be 
interpreted as a network theory. One can, how- 
ever, bring his calculation into correspondence 
with the network theory outlined in the pre- 
ceding paragraph by modifying his method of 
counting configurations in a way that does not 
affect the final result for Z, and then reinter- 
preting his formalism. In particular, it is neces- 
sary to interpret G, which in Wall’s paper is the 
number of molecules per unit volume, as the 
number of segments of the active network per 
unit volume. Influenced by Wall’s theory, Flory 
and Rehner have also obtained Eq. (24) by 
assuming mean segment extensions which would 
follow from Eq. (21). 

In relating his formalism to the network theory 
of rubber Kuhn employs the idea that the net- 
work junctions are essentially fixed in space. 
(We shall emphasize elsewhere how poor an 
approximation this is to the actual case.) If the 
molecules of Kuhn’s theory are identified with 
the segments of the network, and are considered 
to be constrained only at the ends, then one 
obtains® Eq. (23). As thus applied, the theory 
assumes that the fixed (and therefore the vector- 
mean) extensions of the segments have the dis- 
tribution given in Eq. (20). Kuhn’s original 
formalism yielded instead 


K=7/3Go, (25) 


where Go is the original number of molecules per 
unit volume. This result differs from Eq. (23) 
because it corresponds to a different picture 
of the process of cure. In effect, Kuhn treats 


°F. T. Wall, J. Chem. Phys. 10, 485 (1942). 

°W. Kuhn, Kolloid Zeits. 76, 285 (1936). Kuhn's 
original formalism is brought into closer correlation with 
the network theory of rubber in his recent papers; see 
particularly W. Kuhn, J. Polymer Sci. 1, 380 (1946). 
1943) Flory and J. Rehner, Jr., J. Chem. Phys. 11, 512 

* This was first shown by L. R. G. Treloar, Trans. 
Faraday. Soc. 39, 36 (1943). 


RUBBER 673 


each molecule as though its ends were converted 
into fixed points during the process of cure, 
the instantaneously existing extensions in the 
rubber liquid being converted into fixed ex- 
tensions of the solid. These extensions would 
have the distribution of Eq. (21). Next, the 
midpoints of these molecules are suddenly con- 
verted into fixed points which divide each of the 
molecules into two segments. These segments 
will not have had distributions of extensions 
described by Eq. (20), since they are not portions 
of a completely unconstrained molecule; thus the 
resulting distribution of fixed extensions of these 
segments will not be given by Eq. (21). Finally, 
the quarter-points are suddenly converted into 
fixed points, each original molecule thus being 
divided into four segments, with a distribution 
of extensions different from that of Eq. (21), but 
related to it. If Eq. (25) is rewritten in terms of 
the number of segments of the completed model 
it becomes 

K =7/12G. (26) 


Kuhn has interpreted Go as the number of seg- 
ments between points of chemical bonding, and 
G as the number of segments into which the 
network is divided by chemical bonds and ‘“‘steric 
junction points.”’ In any case, the G segments in 
question are treated as extending between fixed 
points, and the difference between the constant 
coefficients in Eqs. (23) and (26) is due to the 
different distributions of segment extensions im- 
plicit in the formalisms in the two cases. 

All derivations of Eq. (24), with its specific 
value of K, depend on the assumed identity of 
two distribution functions: (1) the distribution 
of extensions of molecular segments in solid 
unstretched rubber, and (2) the distribution of 
instantaneous extensions which these same seg- 
ments would have if they were completely 
unconstrained. This hypothesis appears to have 
been considered to be self-evident by its propo- 
nents, and no justification for it has been pro- 
posed.® It is, however, open to very serious 
criticism, and it appears to us that it cannot be 
relied upon. Our objections to it are connected 
with the fact that physically occurring networks 


®It is not always clearly stated that the extensions 
considered in the solid are the vector-mean extensions, 
but this interpretation is required if the result is to be 
understood in terms of the network picture. 
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are not formed all at once, but by a process of 
gradual growth and reinforcement which we 
shall now describe. 

The first stage in the cure of a liquid of rubber 
molecules is that of gel formation. Flory” has 
shown how this occurs as long flexible molecules 
are bonded together. First these molecules are 
linked together into more complicated molecules 
containing two, three, or more chains. With 
increasing number of bonds, the more compli- 
cated molecular species occur with increasing 
frequency. When the number of bonds becomes 
sufficiently large there suddenly appears an 
essentially infinite molecule—a gel, or network 
extending throughout the material. This is at 
first a very tenuous structure, with low tensile 
strength and very small rigidity. 

The second stage of cure, in which the gel is 
built up into a tough and relatively rigid ma- 
terial, is that with which we are here concerned. 
The amount of material in the network is in- 
creased when bonds are formed between the 
network and previously unattached material. 
Formation of these bonds does not increase the 
rigidity of the material, for a structure linked to 
the network by a single bond is only a “loose 
end’’; it is the formation of bonds between pairs 
of elements within the network that increases 
the rigidity of the system. These bonds may link 
two points of the active network, thus increasing 
its rigidity but not the amount of material 
involved in it; or they may link a loose end to 
the active network at a second point, or to 
another loose end, thus incorporating new ma- 
terial into the active network. All these processes 
increase both K and the number of segments in 
the network, though not evidently in the same 
ratio. Formation of a bond between two points 
in the same loose end increases neither the 
amount of material in the active network nor 
the value of K. 

Our objections to the assumption of Eq. (21) 
are then the following: 

(1) The rubber network is built up by pro- 
cesses which involve segments constrained by 
their inclusion in the network—segments to 
which the distribution of Eq. (20) does not apply. 
In fact, it is only the formation of intra-network 
bonds that causes an increase in K. 


10P. J. Flory, J. Am. Chem. Soc. 63, 3096 (1941). 





AND E. GUTH 









(2) The mean extension of a molecular seg- 
ment in the network is not likely to be equal to 
its instantaneous extension at the moment when 
it was incorporated into the network. The mean 
extension after union may be greater or less than 
that at the moment of union ; one can reasonably 
expect it to be less, on the average, since the 
chain is not being linked to a rigid framework but 
to points of a yielding structure, which it can 
pull toward each other. 

(3) Even if the segments had extensions de- 
scribed by Eq. (21) immediately after their 
incorporation into the network, they would often 
be subdivided into smaller segments as the 
process of cure continues. It is easily seen that 
after such subdivision has occurred the distribu- 
tion of extensions of the resulting segments will 
not be that of Eq. (21). For example, Kuhn’s 
theory, which can be interpreted in terms of such 
successive subdivisions of molecular chains, does 
not lead to Eq. (23). 

In view of these factors it seems evident that 
Eq. (21) predicts segment extensions that are of 
the right order of magnitude, but that it would 
be rash to assume any higher degree of validity 
for that equation. Correspondingly, Eq. (24) 
should be right as to order of magnitude, but 
subject to change in the numerical factor. We 
have therefore re-examined the problem from a 
point of view which avoids these difficulties. 


III. CHANGE OF K BY BOND FORMATION 


The essential features of our calculation of K 
as it depends on cure are (a) its use of the more 
realistic picture of the process of cure, as outlined 
in the preceding section, and (b) the application 
of the distribution of Eq. (20) only where this is 
logically justified—that is, in the case of uncon- 
strained molecular segments. Our fundamental 
assumptions are as follows: At any time during 
cure the material will consist of a molecular 
network and unattached material. We shall 
assume that: 

(1) The network can be treated as a Gaussian 
network, as defined in reference 2. In particular, 
it will be assumed that every unbranching chain 
segment of the network—and of the unattached 
material—has a configuration function described 
by Eq. (2). This assumption limits our discussion 
to the ideal soft rubbers for which Eq. (1) can be 
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deduced. It then follows from the preceding 
analysis (and will appear again in our further 
development of the theory) that the rigidity of 
the material is increased during cure only when 
a bond is formed between two elements within 
the network. 

(2) The probability that a bond will be formed 
between two bondable elements of the system is 
proportional to the probability that these two 
elements, during the thermal agitation of the 
material, will lie in the same small volume in 
space. This latter probability can be derived, 
for arbitrary pairs of elements of the general 
Gaussian network, using results of reference 2. 

It appears that these two assumptions, to- 
gether with an assumption as to the lengths of 
the unbranching molecular chains in the material 
before cure, should completely determine the 
solution of the problem; in particular it should 
determine K as a function of the total number 
of bonds formed in the material during cure. 
The complete analysis of this problem offers 
formidable mathematical difficulties connected 
with the division of the complete network into 
the ‘active network” and “loose ends.’’ Our 
more limited approach yields a relation between 
K and the number of bonds formed between 
elements of the network—not the completed 
network, but the network as it exists at the time 
of bond formation. This is a fraction of the total 
number of bonds that we can only estimate. 

Our analysis proceeds in three steps, as follows: 

(a) We .consider the nétwork as it exists at 
any stage in the process of cure, and compute the 
relative probability that the next bond to be 
formed within the network will be formed be- 
tween any two given points of the network. This 
is found to be a simple function of two parameters 
that give a sufficient description of the relation 
of the two elements: p, the separation of their 
mean positions in the material, and 9;;, which 
determines the variability of their separation. 

(b) We next compute the change in entropy 
of the system which would result from the 
formation of such a bond. This again depends 
essentially on p and N;;. 

(c) Finally, we compute the average change 
in K per bond formed between elements of the 
network, at this stage of the cure. The averaging 
process requires the assumption of a distribution 
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function g(p, Jt) such that 
dn = g(p, N)dpdMN (27) 


in the relative number of bondable element pairs 
in the whole network characterized by parameter 
values p to within dp, and ® to within dM. 
Fortunately, because of the special character- 
istics of Gaussian networks, it turns out that the 
result of the averaging process is the same for all 
distribution functions g(p, Jt) of a large class; it 
is the same in the opposing limits of very small 
and very large network constraint, and it is 
essentially the same for intermediate cases when 
the network is of regular structure. 

One is thus led with considerable assurance, 
though not with mathematical certainty, to the 
conclusion that the average change in K, per 
bond formed within the network, is equal to 1, 
and this independently of the detailed structure 
of the network or the state of cure of the ma- 
terial. It follows that 


K=B,, (28) 


where B, is the total number of bonds which, 
during the process of cure, were established 
between pairs of elements in the network as it 
then existed, but exclusive of those established 
between elements of the same loose end. Dis- 
cussion of the differences between this result and 
Eq. (23) will be deferred to Section IV. We now 
indicate the calculation in more detail. 


A. Probability of Bond Formation Between 
Elements of a Gaussian Network 


Let &, n, ¢ be the relative Cartesian coordinates 
of network elements i and ;: 


E=xXi—X,, N=Yi-Yi, F=8:—2; (29) 


In reference 2 it has been shown that the proba- 


bility that these relative coordinates have values 


§ to §+dé, n to n+dn, ¢ to (+d¢ is 
dR=R:;(é, n, ()dtdnde, (30) 


where the probability density function Rj; is 
given by 


‘ | 3 | 3 1 
ii(€, 9, 0) = 2aPN,; exp 2F N:; 


XL(E= fo)" (no? + 6-094] (31) 
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In this expression £, mo, fo are the mean and 
most probable values of & 7, ¢, and Mi; is a 
parameter which determines the breadth of the 
distribution of the actual fluctuating values 
£, n, ¢ about these mean values. Given a com- 
plete description of the network, the values of 
fo, no, £0, Itiz can be calculated by methods 
described in reference 2. 

The probability that the element 7 will lie in 
a small volume AV about the element j—in a 
volume AV about the point = 7={=0—is 


s 
2 





3 
dR=R;;(0, 0, O)A v-| 
2n?PN:;; 
| 3 1 ; 
er — alert t rela . (32) 


ij 


The values of &o, 70, 0 will depend upon the state 
of stretch of the material. Cure will of course be 
carried out when the material is unstretched, 
and the values of £0, No» £0 are E00, N00) S00; 
respectively. According to our second funda- 
mental assumption, the relative probability that 
a bond will be formed between these elements is 
then 


9 


feo[-e}. os 
i NR ;) 





A Mi; -| 
ts 2aPN; 5 


where 


eo Léo0?+ noo + Foo" _}? (34) 


is the distance between the mean positions of 
the elements when the material is unstretched. 


B. Change of Entropy by Bond Formation 


Before introduction of a bond between ele- 
ments z and j of the network, the relative 
coordinates &, n, ¢ of these junctions can take on 
all values. Introduction of a bond restricts the 
possible configurations of the network to those 
in which &, n, ¢ lie in a small range of values near 
zero, of total volume AV in (&, , ¢)-space. Now 
R:,(é, 7, £) is proportional to the number of 
configurations of the network consistent with the 
indicated values of the relative coordinates of 
the two elements. Introduction of a bond be- 
tween elements 7 and j will thus change the 
number of configurations possible for the network 






by a factor 
g=R,,(0, 0, oav / 


+@ +o +2 


The resultant change in the entropy of the sys- 





tem is 
AS=k Ing=k In —-)a v| 
2PM; ; 
3k 1 
— ee fo]. (36) 


The difference in the entropy of the network 
with the bond and without the bond falls into 
two parts. The first term on the right of Eq. (36) 
is independent of the form of the material; the 
second depends on the form of the material 
through the quantities £, no, fo. It follows from 
Eq. (8) that 


f= £ooL:, no= noolsy, fo= Cool :. (37) 





Thus 
3 \3 
AS=k in| (— ) av| 
2PM; ; 
3k 1 
—— —[fto Lot nor Ly + forL.”}. (38) 
22 N:; 


Only the second, shape-dependent term in the 
entropy will be of further interest to us, since it 
alone contributes to the rigidity of the material 





C. Average Change in K per Bond Formed 


We have now to find the average change in 
the shape-dependent entropy of the material, 
and the corresponding change in K, per bond 
formed in the material. 

First of all, we restrict our attention to element 
pairs with the same values of Jt and p*. These 
element pairs will have all orientations in space 
with equal probability when the material is 
unstretched. Averaging over all these element 
pairs, we have then 


(£007)av = (noo?)av = (f007)av = 3 P*- (39) 


It follows that the average change in the shape- 
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dependent entropy of the material, per bond 
formed between element pairs with the same 
values of Jt and p, is 


(AS)u = —(k/2P)(e?/M)[LP+L/+L]. (40) 


The corresponding change in the stress Z re- 
quired to produce unilateral stretch of a unit 
cube at constant volume is 


(AZ )w = (0°/PM) -RT LL, —1/L."). (41) 


Comparison with Eq. (1) shows that the average 
change in K, per bond formed between elements 
with given Jt and p, is 


AK = p?/PI. (42) 


Equation (42) shows that certain types of 
bonds will give AK =0 and wil! make no contri- 
bution to K. Bonding of one piece of unattached 
material to another, or to the network, makes 
no contribution to K because the bonded ele- 
ments were previously subject to no constraints 
on their relative motion: Ji= 2%. Bonds within 
the same piece of unattached material, or be- 
tween elements in the same loose end of the 
network, will not contribute to K because the 
bonded elements had previously the same mean 
position: p=0. In the calculation of K one is 
thus primarily interested in bonds formed be- 
tween elements of the network, exclusive of those 
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formed between elements of the same loose end; 
only these bonds contribute to K. The total 
number of such bonds introduced during cure, 
per unit volume of the material, we denote by 
B,; this can be expected to be nearly equal to 
the total number of bonds formed between ele- 
ments-of the network. B will denote the total 
number of bonds formed in unit volume of the 
material, without regard to whether or not they 
contribute to K. In a similar way we shall 
distinguish two distribution functions, g(p, It) 
and ga(p, Jt). The function g(p, Jt) describes, as 
shown in Eq. (27), the distribution of values of 
p and Yt for all element pairs in unit volume of 
the material; ga(p, Jt) describes the distribution 
of p and Jt among element pairs within the 
network only, exclusive of those in the same 
loose end—in short, pairs of elements such that 
a bond between them would increase K and B,. 

We can now compute the average change in K 
per bond formed in a unit cube of the material: 
(dK /dB),. The relative probability that the 
bond will be formed between elements with 
given p and Yt is equal to the product of the 
number of such element pairs in the material, 
as described by the distribution function g(p, It), 
and the probability per element pair that a bond 
will be formed, as given by Eq. (33). Then the 
required average is 


f tei f ANe(p, NA(p, RAK) 


(dK /dB)y =-— 





f ty J dNe(p, NA(p, N) 


’ 


f ai J dNe(p, N) PN? exp | — (3/2E)(p2/N)} 
0 0 





(dK /dB) », a 


J aof areco, 2)M-" exp! —(3/28)(6/N} 


Similarly, one has 


1 i dof ANga(p, N) PN? exp { — (3/21) (p?/M)} 
0 0 





(dK /dBa)wy _ = 


f dp f ANgalp, NN-2? exp{ —(3/2P)(p?/N)} 
0 0 
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Equation (45) will be more useful to us than 
Eq. (44), because of the simpler form of ga(p, It). 
The function g(p, It) differs from ga(p, Jt) in that 
it applies not only to the element pairs within 
the network, but also to a large number of 
element pairs with p=0 or Jt= %, because of the 
presence of unattached material and loose ends. 
One can. reasonably expect ga(p, Jt) to be a 
smoothly-varying function of its arguments; 
g(p, It) differs from ga(p, Jt) by the addition of 
integrable infinities at p=0 and J{t=0. The 
numerators in Eqs. (44) and (45) are identical, 
since g(p, Jt) and ga(p, It) differ only at points 
where other factors in the integrand vanish; 
the denominators differ by a factor (dB,/dB)w. 
In the early stages of cure, when a relatively 
small fraction of the material is included in the 
network, this factor will be small; in the later 
stages, as more and more unattached material 
and loose ends are incorporated into the active 
network, it will approach the value 1. To avoid 
a discussion at this point of the distribution of 
material between the active network, unattached 
material, and loose ends, we shall concentrate 
our attention on Eq. (45). 

If the form of ga(p, Jt) is regarded as subject 
to arbitrary choice one can make the right-hand 
side of Eq. (45) take on any chosen value. On 
the other hand, complete statistical analysis of 
the process-of cure should lead to a definite 
prediction of the form of g, at each stage of the 
cure, and a corresponding value of (dK /dB,)w. 
In default of this complete analysis, we have 
been lead to consideration of a number of special 
models. It is a surprising fact that these models, 
though they represent quite different limiting 
cases, lead to the same value of (dK /dBa)x—a 
value which can therefore be regarded with some 
assurance as representing the result of a more 
complete analysis, at least approximately. We 
turn now toaconsideration of these special models. 


D. Special Forms of the Distribution Function 


In physically occurring networks the mean 
positions of the network elements will be dis- 
tributed through space with uniform average 
density. Correspondingly, one would expect that 
the relative coordinates £0, noo, foo of pairs of 
these mean positions will be uniformly dis- 
tributed. This is not strictly true if the system 
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contains loose ends and unattached material 
which contribute a disproportionate number of 
element pairs with £00=o0=fo00=0. It is easily 
shown, however, that it is true for the active 
part of any random Gaussian network, and it 
remains true if loose ends are randomly attached 
to such a network, so long as one does not count 
pairs of elements in the same loose end. Thus 
we must expect a uniform distribution of values 
of £o0, noo, foo among the element pairs described 
by the distribution function g,(p, Jt). The corre- 
sponding distribution of values of p is propor- 
tional to p*. Physically reasonable forms of g, 
must then all have the property that 


P(p) -{ AN ga(p, M) =cp*, (46) 
0 


where ¢ is a constant. 

Independent Elements. The first limiting case 
to be considered is that in which all bondable 
elements take on their positions essentially inde- 
pendently of each other. This limit is approached 
by a network formed of very long molecular 
chains bonded together at a relatively few points. 
Each element of such a network will be sur- 
rounded for the most part by elements of other 
chains—elements which can be reached only by 
a relatively long path through the network. The 
relation of these chain elements is then almost 
the same as if they were separate molecules in an 
ordinary liquid: the probability that element 1 
will be found in a given position is independent 
of the position of element j, except as this may 
be affected by the direct steric hindrances that 
come into play only when the elements are 
immediately adjacent. In neglecting the rela- 
tively small probability of bond formation be- 
tween elements of the same network segment, 
and the very small correlations between the 
motions of elements of different segments of the 
network, one obtains the limiting case of inde- 
pendent elements. 

It is shown in reference 2 that the probability 
distribution for the coordinates of the 7th element 
of a Gaussian network is of the form 

3 | | s 3 
P (xi, Vi, 3;) = 4——f e 
; 2aPN, | 


exp > 
Pl ore, 
X C(x; — x0)? + (¥i— Yor)? +(8:—30:)?]}, (47) 
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where Jt; is a parameter which determines the 
breadth of the distribution of the coordinates 
Xi, Yi, 2; in the same way that Jt;; does that of 
the relative coordinates £, , ¢ of the ith and jth 
junctions. A similar parameter Yt; describes the 
breadth of the distribution of coordinates of the 
jth junction. It can be shown (reference 2, 
Section VII) that when the ith and jth junctions 
are independent (the limiting case of weak 
coupling), then 


NijG=NAMN,, (48) 


whatever may be the value of p. The value of 
N%; will vary from element to element, and 9;; 
will vary from pair to pair, but there will be no 
correlation between the values of Jt;, and p: the 
distribution of values of 3t;; will be the same for 
all p. This, together with Eq. (46), shows that 
for this limiting case 


galp, NR) =p? f(N), (49) 


where f(Jt) is a function of still undetermined 
form. On substituting Eq. (49) into Eq. (45) 
and carrying out the integration over p, one finds 


(dK /dByw = SANF(N)/SANFM) =1, (50) 


whatever the form of the function f(J). 
Strongly Coupled Elements. We next consider 
the opposite case to the preceding, that of pairs 
of elements as strongly coupled as is consistent 
with the separation of their mean positions. The 
strongest possible correlation in positions is 
shown by elements of the same chain segment. 
We therefore consider first a single Gaussian 
chain of N links, with ends fixed at the points 
(0,0,0) and (0,0, Z). Let elements 7 and j of 
this chain be separated by a segment of n links. 
It is easily shown (see, for instance, Section IX, 
reference 2) that for this pair of elements 


p=(n/N)L, Nij=n(N—n)/N. (51) 
lf the chain is very long, N>>n, then 
MN ;/p=N/L. (52) 


To this approximation the value of t/p for the 
strongly coupled pairs of elements in the chain 
is independent of p, just as the value of 9 is 
independent of p in the case of independent 
elements. 

In a network consisting of many long segments 
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‘the ratio 3t/p will vary from segment to segment. 


There will then be a distribution of values of 
N/p for pairs of elements within the segments, 
but it will be the same distribution for all values 
of p. (Strictly this will be true only so long as p 
is less than the extensions of the long segments. 
Since only pairs of elements with small p are 
likely to be bonded together, this limitation is of 
little importance.) If we pass to the limit of 
strong coupling by assuming this to be true for 
all element pairs in the material, we have 


ga(p, Nt) =pf(N/p) ; (53) 


the form of the multiplying factor in p is deter- 
mined by Eq. (46). On introducing this expres- 
sion into Eq. (45), changing the variables of 
integration to 


x=N/p, y=p?/N, (54) 


and integrating over y, one finds 
(dK /dBa)w = Sdxf(x)/fdxf(x)=1, (55) 


whatever the form of f(x). 

A generalized form of Eq. (53) leads to the 
same result. If the distribution of values of 
N/p” is independent of p, then 


£u(p, N) =p?-"F(N/p”). (56) 
If n#2, introduction of the variables 


x=J/p", y=p?/N, (57) 


and integration over y leads again to Eq. (55). 
For n=2 the integrals do not converge. 

Regular Cubic Network. The relation between 
the limiting cases treated above becomes especi- 
ally clear on consideration of a special network 
—the regular cubic network discussed in detail 
in reference 2, Section IX. This consists of 
segments, all of M links, connected into a net- 
work with the connectivity of a cubic lattice. 
This network extends between the faces of a 
cube with edges of length LZ, on which the fixed 
points of the network are regularly arranged; 
the faces of the cube will be taken to be the 
planes x=0, x=L, y=0, y=L, z=0, z=L. The 
number N of segments in each chain extending 
between opposite faces of the cube will be 
assumed to be large. Each junction of the net- 
work will be distinguished by three indices, r, s, ¢, 
which give the number of network segments 
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lying between that junction and the planes’ 


x=0, y=0, z=0, respectively. Each of these 
indices can range from 1 to NV—1. 

It is shown in reference 2 that the mean 
coordinates of junction (r, s, t) are (rL/N, sL/N, 
tL/N). The separation of the mean positions of 
junctions (7, s, /) and (u, v, w) is thus 


p(r, s, t; u, v, w) =(L/N)[(r—u)? 
+(s—v)?+(t—w)? }}. 


The parameter Yt;-of Eq. (47) is essentially the 
same for all junctions except those very near the 
fixed points on the bounding surfaces: 


Nir, s, NLM /4. 


(58) 


(59) 


The parameter Jt;; for junctions (r,s, ¢) and 
(u,v, W), not too near the bounding surfaces, is 


Nir, s, t; u, v, w) 
M 1 1 


~ 21 e[(r—w)?-+(s—v)?+ (¢—w)?} 





oz 


(60) 


For the junctions of this network one has then 
the relation 


N= M/2{1—(1/4N)(L/p)}. (61) 


If the junctions of this network were inde- 
pendent one would have, by Eq. (48), Jt=//2. 
Thus Eq. (60) shows ‘that the coupling of 
immediately adjacent junctions of the network 
has an appreciable effect on the distribution of 
their relative coordinates, reducing Jt by some 
30 percent. This coupling, however, decreases 
rapidly with increasing difference in the ordinal 
numbers of the junctions; junctions more than 
a few segments removed from each other take 
on their positions almost independently. This is 
true whatever the values of L and N. The 
condition for effective independence of two ele- 
ments is thus that p>L/N. If L is relatively 
small—that is, if the network is not greatly 
extended—junctions with small p will neverthe- 
less be essentially independent. On the other 
hand, it is clear from Eqs. (33) and (60) that 
bonding is likely to occur even between elements 
with p> { MP}}, which is the r.m.s. separation of 
the ends of an unconstrained chain segment of 
M links. Thus mos’ of the bonds introduced at 
random into this network will link essentially 
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independent elements if 


L/N«{MP}}; (62) 


the network elements can be treated as inde- 
pendent if the network constraints are so small 
that the segments have vector-mean fractional 
extensions small compared with the ordinary 
r.m.s. average extensions they would have if 
they were unconstrained. 

It is easily shown that Jt for two elements in 
the same segment of a regular cubic network 
approaches the number of intervening links as 
this number becomes small. In this limit one finds 


N=n, p=(n/M)-(L/N). (63) 


Such elements are strongly coupled in the sense 
already discussed. If 


L/N>{MP}}, (64) 


bonds introduced at random will for the most 

part link such strongly coupled element pairs; 

the limiting case of strongly coupled elements is 

approached by a regular cubic network in which 

the network constraints are strong and result in 

large vector-mean extensions of the segments. 
The intermediate case, 


L/N={ MP}} (65) 


probably corresponds more closely to physically 
occurring networks than either of the limiting 
cases. It can be discussed as follows. Equation 
(61) strictly applies only to pairs of network 
junctions, for which p is not a continuous vari- 
able, and is never less than L/N. If we consider 
all element pairs, Jt must be a continuous func- 
tion proportional to p when this becomes small. 
The relation 


ye— 


(66) 


2 


p+L/Nr 


M | p 


approximates Eq. (61) very well wherever that 
equation is applicable. For small p it gives 

=7/2(MN/L)p, whereas Eqs. (63) show that 
for element pairs in the same segment one should 
have 3t&(MN/L)p; nevertheless it can be ex- 
pected to be a good representation of the situa- 
tion in a regular cubic network when all element 
pairs, whether in the same or in different network 
segments, are taken into account. 
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The form of g. corresponding to Eq. (66) is 


M p 
aie 2) = #4(R—— —__}), (67) 
2'p+L/Nr 





where 6 is the usual singular function which gives 
unit integral over any range including the origin 
but vanishes except when the argument is zero. 
Equation (45) then becomes 


1 J dpp‘|(p+a) /cp}** exp{ — (3/21) p*(p+a) /cp} 
0 





(dK /dBa) wy “; 


«© 
0 


where 
c=M/2, a=L/Nrz. (69) 
On changing the variable of integration to y, with 


? = (3/2/c)p(p+a), (70) 





(68) 


i) dpe'L(p-+a)/cp}*? exp{ — (3/2) p*(p-+a)/cp} 





and introducing the parameter 
b= {3a?/2cl?}4 = {3L?/rN?MP}}, (71) 


one obtains 


f dyy* expl—y?]{1+[b/(62+492)!]} 


(dK /dBa)w = 3 





In the limit of independent elements, 6-0, or 
strongly coupled elements, b>, the last factor 
in the two integrals of Eq. (72) becomes constant 
over the range of integration, and one obtains 
the results already noted in Eqs. (50) and (55). 
For other positive values of } this factor decreases 
monotonely from 2 to 1 as y increases from 0 to 
*. The presence of this factor decreases the 
numerator by a larger factor than it does the 
denominator, since the larger values of y con- 
tribute relatively more to the numerator. Thus, 
one can easily prove, 


(dK /dBa)w <1 (73) 


for all positive values of 6. The intermediate case 
described by Eq. (65) corresponds to b=3, for 
which numerical evaluation of Eq. (72) yields 
(dK /dBa)w=z0.96. The minimum -value of this 
quantity, ~0.92, is obtained with 6=1.5. 

It thus appears to be a good approximation 
in any case, and a safe upper limit, to take 
(dK /dBa)w=1. Integration of this equation with 
respect to B, yields the already stated result of 
this approach to the problem: 


K=B,. (28) 


f dyy? expl—y"]{1+[b/(b-+4y*)*J} 
0 


(72) 





IV. COMPARISON OF THE THEORETICAL RESULTS 


We have now to discuss the relation between 
the result of Wall’s theory, K =G, and the result 
obtained by us in the preceding section. 

Let the system to be considered consist origi- 
nally of N unbranching molecular chains. Each 
bond introduced into the system will increase 
the total number of chain segments by two: a 
bond formed between two different segments (or 
chains) will divide each segment into two, and a 
bond formed between two elements of a single 
segment will divide this into three. If B bonds 
are introduced, the total number of segments in 
the system will become 


G=2B+N. (74) 


If Wall’s postulate concerning the distribution 
of extensions of molecular segments [Eq. (21) ] 
is to be introduced into the network theory, it 
must certainly apply only to the vector-mean 
extensions of the segments of the active network; 
all other segments, whether of loose ends or 
unattached material, have \=0 and contribute 
nothing to K. Let the active network, stripped 
of all loose ends, contain G, segments, and let 
there be J junctions in this network. Each of 
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these junctions will form the end point of either 
three or four segments. It is then clear that 


$I <Ga<2J <2B. (75) 


Each chain attached to the network will include 
at least two segments, the terminal segments, 
which are not part of the:active network. If all 
chains are attached to the network one has also 


Ga<G—2N=2B-N. (76) 


The degree of inequality will depend on the 
nature of the network; the difference will be 
large if loose ends make up a considerable part 
of the structure. 

The quantity B, of Eq. (28) is unlike G, in 
that it cannot be determined by observation of 
the completed network. B, is the number of 
bonds introduced as intra-network bonds when 
they were formed; examination of simple cases 
will show that a given network can be built up 
by alternative processes for which B, has differ- 
ent values. One may term G, a “‘structural”’ 
quantity, since it can be determined from the 
network structure, and B, a “historical’’ quan- 
tity, since to determine it one must know the 
history of the network—that is, the process by 
which it is formed. The relation between B, and 
G,, like that between B, and K, or G, and K, is 
thus only a statistical relation. 

Concerning the relation between B, and Gg, 
the following remarks can be made. A_ bond 
formed between different loose ends will con- 
tribute 1 to B,, and 1 or more to G,; the increase 
in Gq may be large if the loose-end structures are 
complicated. A bond formed within a single 
segment of the active network will increase B, 
by 1, but will not change G,; the loop thus 
formed in the segment is a loose end, and when 
it is stripped off the added bond will not sub- 
divide the original segment. A bond formed 
between two different segments of the active 
network, or between a loose end and a segment 
of the active network, will increase B, by 1, and 
G, by two or more. Processes of the last types 
will certainly dominate in the latter stages of 
network formation, and it is reasonable to expect 


Ga=2B,. (77) 


It thus appears that our theory predicts a value 
of K which is less by a factor of about 2 than 
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that given by Wall’s 
placing G by G,. 

A comparison of theory with experiment has 
been carried out by Flory." Flory treats the 
inequality of Eq. (76) as an equality, but finds 
for a series of Butyl rubber vulcanizates that the 


theory, amended by re- 


1 Pp, J. Flory, Chem. Rev. 35, 51 (1944). In this review 
article Flory makes a number of incorrect or misleading 
statements about our work which we take this opportunity 
to correct. , 

(1). Flory states that Eq. (3.3) of reference 1 is in error 
by a factor of 3. Actually it is correct, and in agreement 
with the result which he gives in his own paper. As was 
noted in our explanation of Eq. (2.5), (L?)a, is, in the three- 
dimensional case, the mean-square value of a component of 
the extension of the chain; in Flory’s notation it is 1/28°. 

(2). As concerns a very simple network model introduced 
in reference 1 for purposes of visualization, Flory seriously 
misstates the justificatioa given for its use and gives a 
false idea of the importance to be attached to its intro- 
duction. Instead of presenting as a justification the quite 
trivial observation indicated by Flory, we gave an extensive 
analysis of the behavior of a general network of Gaussian 
flexible chains. Along with other results, this led to the 
conclusion that the entropy forces exerted by an arbitrary 
network, and hence the entropy itself, was the same as 
that for a very simple “network”’ of independent chains 
from which cross linkages were entirely absent. For 
illustrative purposes in certain derivations we used this 
simple model, but all results hold for the general model 
with Gaussian chains. 

(3) In speaking of our use of ‘‘springs’’ Flory apparently 
expresses a distaste for a discussion of the problem in 
terms of average forces rather than of thermodynamic 
functions. In statistical mechanics these methods are of 
course equivalent. 

(4) Flory states that certain equations relating to this 
model are dimensionally incorrect. Actually they are 
correct, though the notation may disguise this fact. J/ is 
the effective number of chains crossing unit area in the 
unstretched material; its dimensions are reciprocal length 
squared. N is the number of links per unit length of 
molecular chain in the unstretched state of the material; 
it has the dimensions of reciprocal length. In equations 
involving this model (L*),4,= 3N/? has therefore the dimen- 
sions of length. The physical significance of M and N is 
discussed in the Appendix to reference 1, p. 481. 

(5) After noting that almost all the chains of a well-cured 
piece of rubber will be linked to the network, Flory 
indicates that this result is in contradiction to our estimate 
that 25 percent of the total material is ‘ ‘actively involved” 
in the network. We do not wish to attach undue signific: ance 
to our estimate, which was admittedly rough. It is, 
however, necessary to note that Flory’s estimate and ours 
do not apply to the same quantity. Flory’s statistical 
theory lumps together all material of all chains connected 
into the network, whether by one bond or by many. In 
the ‘active part of the network,” with which we were 
dealing, one has to do only with portions of chain which 
contribute to the forces. A molecule linked to the network 
by only one bond does not contribute at all to the active 
part; molecules linked at two points contribute, on the 
average, about half their length. 

(6) Our failure to arrive at a definite value of K in 
reference 1 was not due to such use as we made of the 
independent chain network, but to our desire to avoid 
unjustifiable postulates. One can introduce special postu- 
lates which lead to a value of K into that model quite as 
readily as into any other model. We prefer the approach 
to the problem given in this paper. 
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theoretically predicted rigidity, with K=G,, is 
less than the observed rigidity by a factor of 
about 2.1. If Flory’s results are correct there is 
then a discrepancy of a factor of about 4 between 
observation and the theory of Section III. 

Flory has suggested several possible origins for 
this discrepancy between theory and his observa- 
tions. He suggests that entanglements—that is, 
steric hindrances of the sort neglected in the 
theory of Gaussian networks—may have the 
same effect as bond formation. In a paper now 
in preparation we shall discuss the effect of such 
entanglements. It appears to us that their effect 
should be small when the network segments are 
long and the material is but little stretched ; 
their effect should be greater in more highly 
cured materials, and at greater extensions. The 
fact that Flory’s observations are made at 300 
percent elongation would favor increased im- 
portance of these steric hindrances. On the other 
hand, Flory’s stress-strain curves have essentially 
the form of Eq. (1), whereas steric hindrances 
would modify this form; furthermore, the dis- 
crepancy between experiment and theory is less 
for the more highly cured materials, rather than 
greater, as we would expect if the effect is one 
of steric hindrances. 

We attach greater importance to another 
point made by Flory, concerning a peculiarity 
of Butyl rubber—the limited number of points 
at which cross linking can occur. Our theory has 
assumed that there are so many points of 
bonding that all relative positions of bondable 
elements occur with equal probability. If the 
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bondable elements in the chains are widely 
scattered the small values of p will be relatively 
few. In such a case, the average increase in K 
per bond formed may be much greater than that 
given by our theory. For instance, a treatment 
of the regular cubic network of Section III, with 
the assumption that only the junction elements 
are bondable, leads to calculation of a value of 
(dK /dB,)« which may be much larger than 1 if 
the junction points are not close together. A 
repetition of Flory’s experiments, with other 
types of rubber, is much to be desired. 

It seems also to be possible that Flory’s 
estimate of the number of bonds introduced 
during cure is in error by a considerable factor. 
Flory has determined this quantity by use of the 
theory of gelation. Walling” has attempted to 
check the theoretical predictions of Stockmayer'® 
as to the point of gelation, in the case of the 
systems methyl methacrylate-ethylene dimetha- 
crylate and vinylacetate-divinyladipate. Walling 
succeeds in checking the qualitative features of 
Stockmayer’s formulas but notes that ‘‘ observed 
extents of reaction at the gel-point are in general 
found to be times those calculated, 
particularly in experiments in which early gel- 
points are to be anticipated.’ If such a discrep- 
ancy is present in Flory’s observations, his esti- 
mate of the number of bonds formed in cure of 
his materials may be several-fold too low. This 
factor alone would go far toward accounting for 
the discrepancy between theory and experiment. 


several 





2 C, Walling, J. Am. Chem. Soc. 67, 441 (1945). 
1 W.H. Stockmayer, J. Chem. Phys. 12, 124 (1944). 
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will be expressed as —of¢s,, where s, is the entropy of fusion 
per unit when W;—1. Actually, the configurational entropy 
contribution, which may account for the major part of the 
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THIS section will accept reports of new work, provided 
these are terse and contain few figures, and especially 





entropy of fusion, should be taken proportional to ¢—1 
(see Eq. (8) of reference 3). 


3P, J. Flory, J. Chem. Phys. 10, 51 (1942). 
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few halftone cuts. The Editorial Board will not hold itself The total entropy change involved in the formation of Ro 
responsible for opinions expressed by the correspondents. — the crystallite will be 
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High Polymers The equilibrium condition at which T=T,,, the tempera- 6 
II. Simplified Derivation of Melting-Point ture (“melting point”) for complete disappearance of crys- pose 
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PauL J. FLORY either o or ¢ and equating to zero, giving V: 
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shown by x-ray studies, the structural units of polymers If the polymer chains are composed exclusively of A tran: 
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substances. While a comprehensive treatment covering x js the number of repeating units per molecule. Substitut- ina 
various complicating circumstances which may be antici- ing in Eq. (2) and proceeding as above, there is obtained in mani 
pated in crystaliine polymers is somewhat cumbersome, the limit of vanishing ¢ corre 
the essentials devoid of these complications can be pre- 1/Tn—1/Tn9 = R/ xh. TI 
sented briefly as follows. , ester 
PRE apt EE a a ee ee However, the more detailed treatment leads to an addi- oe 
polymer c g ‘ a : was 
. . : tional term arising from the increased concentration of ends 
considered, these units are capable of forming a crystal j aay? prese 
lattice. Other units B, unable to enter the crystal lattice of molecules in the liquid phase. Thus, of th 
characteristic for A, also may be present. It will be as- 1/Tm—1/Tn?=(R/xhy)(1+1/2), (9) the s 
sumed that the structure of the-liquid can be described in where z is the number of configurational segments per in tl 
terms of a pseudo-lattice, each cell of which is occupied by _ structural unit. If the molecules are of random length such press 
a polymer segment. (The configurational segment is as- that the probability of continuation of the chain from one ever, 
sumed to be identical with the structural unit in the present —_unit to the next is 9, may 
treatment.) The formation of a crystallite composed of 1/Tm—1/Tm®=(R/h.)[—In p+ (1—p)/s], (6) W 
o-chains in its cross section and f-units in length, requires which is approximately equivalent to Eq. (5) if x is replaced plot 
first of all the location of a planar array of o-cells each by its number average 2p. deca 
occupied by type-A was. Each of these A units inp be A direct measurement of the heat of fusion of a high theo 
succeeded by an uninterrupted sequence of ¢—1-additional polymer invariably yields a value dependent on the degree with 
A units if the array is to be acceptabie asa nucleus for the of crystallinity, the latter being difficult to determine. By 0s. It 
formation of a crystallite of the prescribed length. The means of the above relationships the latent heat (and the mett 
probability of such a constellation will be W;’, where Ws entropy) for the transformation from totally crystalline prece 
a the probability of a yeas of fA segments ina given + amorphous polymer may be deduced readily from T: 
chain, starting with an arbitrarily selected unit. The corre- mere melting-point measurements on polymers of known = 
sponding entropy term will be ko In Wy. coastivation. Con 
It is further required that the segments (units) involved . : . ’ 
‘ ‘ ‘ ! The work on which this paper is based comprises a part of a program x 
become aligned in perfect array, packed in the more dense _ of fundamental research on rubber and plastics being carried out under I 
arrangement of the crystal lattice. The entropy associated a. of Naval Research and the Goodyear Tire ‘ 
with this step wiil be assumed, as an approximation, to be 2 Contribution No. 142 from the Goodyear Tire and Rubber Company S 
proportional to the number of units, i.e., the entropy change Re aporatny. = 
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Thermodynamics of Crystallization in 
High Polymers 


III. Dependence of Melting Points of Polyesters 
on Molecular Weight and Composition’? 
ROBERT D. Evans, HAROLD R. MIGHTON, AND PAUL J. FLoRY 


The Goodyear Research Laboratory, Akron, Ohio 
July 24, 1947 


LTHOUGH fairly definite melting points have been 

assigned to a number of crystalline polymers and to 
a few copolymers, no adequate systematic study has been 
carried out either on the effects of molecular weight or on 
the influence of copolymer composition on the melting 
point. Recently we have undertaken in this laboratory a 
comprehensive investigation of the dependence of melting 
point and crystallinity on polymer constitution. Linear con- 
densation polymers are especially well suited for this pur- 
pose and we have chosen the polyester, decamethylene 
adipate, for the preliminary experiments reported here. 

Various copolymers have been prepared in which a por- 
tion of the decamethylene or the adipate residue is replaced 
by another glycol or dibasic acid radical. Decamethylene 
adipate polymers, varying in molecular weight and in the 
nature of the end groups have been prepared also. Melting 
points of these polymers and copolymers can be repro- 
ducibly measured using a micro-melting point apparatus, 
consisting of a low power microscope equipped with an 
electrically heated hot stage. The sample, located in the 
hot stage, is observed between crossed Polaroid sheets using 
transmitted light while the temperature is raised slowly. 
So long as crystallinity remains, the sample appears bright 
in a dark field. The fairly sharp disappearance of this 
manifestation of depolarization marks the “melting point,” 
corresponding to complete disappearance of crystallinity. 

The melting points of decamethylene adipate copoly- 
esters in which 20-mole percent of the glycol or acid residues 
was replaced by the indicated comonomer residue, are 
presented in Table I. The same molar percentage of each 
of the first three comonomers lowers the melting point by 
the same amount in accordance with the theory presented 
in the preceding communication.* The melting-point de- 
pression for the copolymer containing sebacate units, how- 
ever, is appreciably less, indicating that the sebacate unit 
may occur within the crystallites. 

When the proportion of the co-ingredient is varied, a 
plot of 1/7 vs. In X4, where X4 is the mole fraction of 
decamethylene adipate units, is linear in agreement with 
theory. Eutectic compositions can be located in analogy 
with simple binary mixtures. From the slope of the 1/7 
vs. In X4 plot a heat of fusion of 3800 calories per deca- 
methylene adipate unit is computed from Eq. (4) of the 
preceding paper.* 


TABLE I. Melting points of decamethylene adipate copolyesters. 








Comonomer group 7 in poises at 218°C m m. pt. (°C) 
None 6800 81 -81.5 
Isophthalate 600 64.5-65.5 
cis-quinitol 310 5 66.5-67.5 
Resorcinol 64 -65 
Sebacate 72 -75 
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TABLE II. Effect of molecular weight on the melting point 
of decamethylene adipate. 








° 


m. pt. (°C) 
mm. pt. (°C) cale from Eq. (7) 


7 in poises at 
109°C 





0.171 42.5-45.5 45.3 
0.616 62 -64 62.0 
1.60 67 -68 67.8 
18.1 72.5-74 71.2 
31.3 74 -75 74.0 


85 est. 6800* 81 -81.5 80.0 








* Viscosity determined at 218°C. 


The dependence of the melting point on the average 
number of structural units (the unit consisting of one 
decamethylene and one adipate radical) is shown in Table 
II for polymers terminated with benzoate groups. In agree- 
ment with Eqs. (5) and (6) of the preceding communica- 
tion,’ 1/T,, varies linearly with 1/2,; i.e., 


1/Tn— 1/353 =7.69 X 1074/Zp, (7) 


from which the melting points given in the last column of 
Table II have been computed. Comparing Eq. (7) with (6)* 


h,/(1+1/zs) = 2600 cal., 


which, in conjunction with the above value for hy, leads to 
a not unreasonable value between 2 and 3 for z. 

Decamethylene adipates having hydroxyl or carboxyl 
(or both) end groups, melt higher than benzoate terminated 
polymers of the same degree of polymerization. Occurrence 
of these smaller end groups within the crystal lattice affords 
a likely explanation for this observation. 

We have concluded from these preliminary experiments 
that melting points of polymers depend in a systematic 
fashion on polymer constitution and that quantitative rela- 
tionships can be traced. The results appear to be in accord 
with the theory outlined above. 

1 The work on which this paper is based comprises a part of a program 
of fundamental research on rubber and plastics being carried out under 
contract between the Office of Naval Research and the Goodyear Tire 
and Rubber Company. 

2 Contribution No. 143 from the Goodyear Tire and Rubber Com- 


pany, Research Laboratory. 
3 See Paul J. Flory, J. Chem. Phys. 15, 684 (1947). 





Resonance Energies from Thermal Data 


T. L. CoTTRELL AND L. E. SUTTON 
Physical Chemistry Laboratory, Oxford, England 
July 3, 1947 


N a recent paper Dewar! discussed the discrepancy 

between the values for resonance energies of benzene 
derivatives found by the Pauling-Sherman method from 
heats of combustion and by the Kistiakowsky method 
from heats of hydrogenation, and suggested that it is 
due to differences in internal kinetic energy between cyclic 
and open-chain compounds. In particular, he predicted 
that the heat of formation of cyclohexane should be 
greater than that calculated from the bond-energy sum 
by an amount equal to the difference between the reso- 
nance energy of benzene obtained by the two methods. 

Recent data show that in fact the heat of formation of 
cyclohexane is the same as that calculated from bond 
energies obtained from the normal paraffins. The dis- 
crepancy between the two ways of obtaining resonance 
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energies of benzene derivatives proves to be due to the 
use of incorrect bond-energy values. No special explana- 
tion in terms of differences in stability between cyclic 
and open-chain structures is necessary, since the required 
bond energies may be derived from open-chain compounds. 

The bond energies in the normal paraffins reach con- 
stant values with increasing chain length. From Rossini’s 
data? for the heats of formation of the paraffins, and values 
of 126.3 and 51.71 kcal. as the heats of formation per 
gram atom of‘carbon and hydrogen, respectively, in their 
standard states from monatomic gases at 298.1°K%, the 
the C—C and C—H bond energies are found to be 60.3 
and 87.1; kcal. The sum of these bond energies gives the 
heat of formation of n-hexane exactly, and using them we 
obtain 1407.6 kcal. for the heat of formation of cyclo- 
hexane from atoms, which is in very good agreement with 
the value of 1407.8 obtained from Rossini’s data.‘ 

With the above C—C and C—H values, that for the 
C=C energy may be obtained from Rossini’s data*® for 
the heats of formation of the monodlefin hydrocarbons, 
which were obtained by combining heats of hydrogenation 
with specific heat and heat of combustion data. For 
olefins with the double bond in the 1 position, this tends 
to a constant value of 101.2 kcal. With these bond energies, 
the apparent resonance energy of benzene is 41 kcal., in 
disagreemeni with hydrogenation result of 36. If, however, 
the heats of formation of the isomeric olefins are considered, 
and the other bond energies taken as constant, the C=C 
energy is found to vary from, for example, 101.3 in 1- 
henene to 107.2 in 2,3-dimethyl-2-butene. This is, of 
course, the well-known substitution effect noted by 
Kistiakowsky,® which persists to a considerable extent at 
0°K.7 The hexene in which the double bond most nearly 
resembles that found in an unsubstituted ring is clearly 
cis 3-hexene, for which the C=C energy is found to be 
102.9 kcal. Using this C=C value, we obtain for the 
resonance energy of benzene_36 kcal. in agreement with 
the hydrogenation value. Similarly for 0-xylene, two of 
the C=C bonds may be taken as equivalent to that in 
2-methyl-2-butene, and the third as equivalent to that in 
cis 3-hexene, which reduces the derived resonance energy 
by 12 kcal. to 34. Again, with Rossini’s value for the heat 
of combustion of styrene,* and appropriate C=C energies 
obtained from non-cyclic compounds as before, the reso- 
nance energy of styrene obtained by the two methods 
agrees. That it is legitimate to consider the double bond 
in a ring structure as similar to a corresponding double 
bond in an olefin is shown by the fact that the heat of 
hydrogenation of non-cyclic compounds of the type 
CHR=CHR (cis) is the same as that of cyclohexene. 

We should point out that while for convenience we have 
assumed that bond-energy differences are localized in the 
C=C bond, this is a purely formal treatment; some change 
may occur in other bond energies. 

1 Trans. Faraday Soc. 42, 767 (1946). 

2 Prosen and Rossini, J. Research Nat. Bur. Stand. 34, 263 (1945). 

8’ Coates and Sutton, to be published. 

— Johnson, and Rossini, J. Research Nat. Bur. Stand. 37, 51 
" 4 Preaeis and Rossini, J. Research Nat. Bur. Stand. 36, 269 (1946). 

® Kistiakowsky et al., J. Am. Chem. Soc. 58, 137 (1936). 


? Rossini, Chem. Rev. 27, 1 (1940). 
8’ Prosen and Rossini, J. Research Nat. Bur. Stand. 34, 59 (1945). 
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The Detection of Free Radicals in 
Hydrogen Atom Reactions 


G. M. HARRIS AND A. W. TICKNER 


Department of Chemistry, University of Saskatchewan, 
Saskatoon, Saskatchewan, Canada 


July 2, 1947 


LARGE number of investigations of reactions of 

hydrogen atoms with various substances have been 
made making use of the discharge tube method of Wood 
and Bonhoeffer.! With organic vapors, the mechanisms of 
such reactions have invariably been explained in terms of 
free organic radical chains.2 However, no direct evidence 
for the existence of free radicals in such reacting mix- 
tures has heretofore been reported. The Paneth-Rice 
technique for the detection of free radicals,’ successfully 
applied in the past to photo-chemical and thermal de- 
composition reactions, has been recently adapted by the 
present authors to the hydrogen atom reaction problem. 

A mirror of metal known to be insensitive to atomic 
hydrogen but sensitive to free organic radicals was de- 
posited by condensation in a vacuum on a water-cooled 
glass finger, and introduced into the reaction chamber of 
the apparatus. Bismuth was adopted as the mirror metal 
in our work because of the ready availability (from uranium 
extraction residues) of its radioisotope of atomic weight 
210, the activity of which is eminently suitable for tracer 
use. Preliminary experiments confirmed the non-removal 
of a radioactive bismuth mirror in the presence of atomic 
hydrogen alone, thus eliminating the possibility of the 
formation of a stable volatile hydride under the conditions 
of our experiments. 

A volatile compound known to be easily decomposed at 
room temperature by atomic hydrogen, dimethyl mercury, 
was selected for trial. Previous work‘ had suggested a free 
radical mechanism to explain the products, which are free 
mercury and a methane-ethane mixture. It was readily 
demonstrated that no direct reaction occurred between the 
dimethyl mercury vapor and the bismuth mirror in the 


* absence of atomic hydrogen within the temperature range 


20° to 45°C. With the mirror maintained at 20°C, which 
was about 5° cooler than the prevailing room temperature, 
the presence of the hydrogen atom-dimethyl mercury re- 
action mixture resulted in a rapid clouding up of the mirror 
with condensed mercury vapor, and within a few minutes 
it was completely concealed by mercury droplets. Radio- 
activity measurements indicated that part of the bismuth 
was transported to a trap cooled in liquid air, some dis- 
tance removed from the reactor, and recoverable from this 
trap as a volatile compound. However, after 20 minutes 
treatment most of the bismuth still remained on the cooled 
finger, apparently amalgamated with the mercury. 

The experiment was repeated, this time maintaining the 
bismuth mirror at a temperature of 45°C. In this case, the 
mirror disappeared almost completely within the 20-minute 
reaction period, and the resulting volatile bismuth com- 
pound was shown to have been transported to the liquid 
air trap by means of radioactivity measurements, as before. 


The mercury formed in the reaction condensed as a grayish 


fog on the surrounding walls of the reactor, which were at 
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room temperature, and consequently about 20° cooler 
than the bismuth mirror. 

Our experiments thus prove that a volatile bismuth 
compound, presumably the trimethyl, is formed when di- 
methyl mercury is decomposed by atomic hydrogen in the 
vicinity of thin deposits of metallic bismuth. The fact that 
the mercury produced in the reaction was prevented from 
amalgamating with the mirror metal when the latter was 
slightly warmer than the reactor walls, provides strong 
evidence that the dimethyl mercury-hydrogen atom re- 
action was a gas phase process, and not merely a hetero- 
geneous reaction on the surface of the bismuth deposit. 
The volatile bismuth compound must result, therefore, 
from direct reaction between alkyl radicals and the mirror. 
Full details of our technique.and of correlative experiments 
now under way in this laboratory will be published else- 
where at a later date. 

1G. Glocker and S. C. Lind, The Electrochemistry of Gases and Other 
Dielectrics (John Wiley & Sons, Inc., New York, 1939), p. 115. 

2E. W. R. Steacie, Atomic and Free Radical Rections (Reinhold 
Publishing Corporation, New York, 1946). 

3F. Paneth and W. Hofeditz, Ber. 62B, 1335 (1929); F. O. Rice and 


K. K. Rice, The Aliphatic Free Radicals (Baltimore, 1935). 
+G. M. Harris and E. W. R. Steacie, J. Chem. Phys. 13, 559 (1945). 





Absorption by Hydrogen Sulfide in the 
3.5u-4.5u Region* 
EarL A. WILSON, JR. AND PAuL C. Cross 
Brown University, Providence, Rhode Island 
July 24, 1947 
HE percent absorption by 90 cm of H2S at 1-atmos- 
phere pressure is plotted against wave number in 
Fig. 1. This spectrum was obtained with a fluorite prism 











WAVE MUMBER © (nce-Secer sect) 


Fic. 1. Absorption by 90-cm atmos. of hydrogen sulfide 
in the 3.5u—4.5u region. 


spectrometer which had a resolution of about 6 cm™ in 
this region. 

Theoretically, it is possible for three vibrational transi- 
tions (vg, vz, and 25) to contribute to this absorption. It 
is especially noteworthy that considerable absorption is 
observed in the region 2620 cm=!—2450 cm™, where the 
P-branch of vz may be expected to appear according to the 
calculations of Hainer and King,' but contrary to the 
reports of Nielsen and co-workers.’ 

A detailed tabulation of the absorption peaks at this 
resolution would be of little use unless accompanied by an 
analysis of the band envelope. A complete analysis will 
probably require new measurements at higher resolution. 
_* This work was performed under a contract between Brown Univer- 
sity and the Office of Naval Research. 

'R. M. Hainer and G. W. King, J. Chem. Phys. 15, 89 (1947). 

*H. H. Nielsen and E. F. Barker, Phys. Rev. 37, 727 (1931); A. D. 


Sprague and H. H. Nielsen, J. Chem. Phys. 5, 85 (1937); R. H. Nobel 
and H. H. Nielsen, Phys. Rev. 71, 484 (1947). 
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The Magnitude of the Experimental Error in 
the Adsorption of Nitrogen at 78°K 
on Various Surfaces 


W. V. LOEBENSTEIN AND V. R. Dettz 
National Bureau of Standards, Washington, D. C. 
July 7, 1947 


HE adsorption isotherms of nitrogen at — 195°C have 

been widely used in recent years to determine the sur- 
face areas of solids. The published values for the surface 
areas have covered the range from about 2000 square meters 
per gram to about 0.1 square meter per gram. If the surface 
area is large, a small sample of about one gram is sufficient 
in the experimental procedure for the adsorption measure- 
ments. However, for materials of small surface areas, it is 
necessary to take large samples and, hence, use adsorption 
tubes with comparatively large dead spaces. It is the pur- 
pose of this note to ascertain that part of the experimental 
error which may result from the use of the various sized 
adsorption tubes which are necessary to accomodate ad- 
sorbents of widely varying surface areas. 

Assuming the perfect gas laws to apply to the gas phase 
measurements, the following expression gives the number 
of moles, n, of gas adsorbed for a given initial point in an 
adsorption measurement. 


n= pi Vze/300R— p. Ve/300R— pe Vn/78R, (1) 


where p;=initial pressure of nitrogen introduced into the 
buret system. 
p-=equilibrium pressure of nitrogen after exposure 
to sample at 78°K. 
Vp=volume of nitrogen (S.T.P.) in the buret system. 
Vp=volume of nitrogen (S.T.P.) in the dead space of 
the adsorption tube. 
R=gas constant (in appropriate units). 
Room temperature is taken as 300°K and the temperature 
of the liquid nitrogen as 78°K. 
If 62 be the variation in the number of moles adsorbed, 
then the general expression for the percentage error is: 


{ (Ve)pit(Va—3.85 Vo) dp. 
1006 i— pe )bV 3.85p.)5V 
even 100 + (pi-— pe)EVa+( >P-) D 2) 
" - Va(pi— pe) —3.85peV 





where 6p;, 5p, 6Vze, and 6Vp are the corresponding varia- 
tions in pi, pe, Ve, and Vp, respectively. It is to be noted, 








Tare I. 
Pressures 
B.E.T. (mm Hg) Percentage 
Surface equi- orer 
Sample area Vu 6Vp_ initial librium Vg  1006n 
(g) (m2/g) (ml) (ml) Pi Pe (ml) n 

Cement* 21.5 1.5 15.3 0.15 104.0 39.2 123. 1.1 
Cotton* 31.1 0.7 121. 1.21 184.7 34.6 125. 10.4 
Bone char (1) 1.92 121. 7.0 0.07 562. 156. 150. 0.20 
Bone char (2) 1.92 121. 7.0 0.07 428. 0.4 47. 0.27 
Glass spheres 12.29 0.93 16.3 0.16 232.8 67.6 50.1 1.9 
Activated 
cocoanut shell 
charcoal 0.71 = 1720. 5.7 0.06 392.7 0.3 149.8 0.12 
Titanium oxide 1.08 9.9 4.1 0.04 1314 789 499 21 





* Data kindly furnished through the courtesy of R. L. Blane of this Bureau. 
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since variations may be positive or negative, that all the 
signs in the numerator of the above fraction are taken to 
be positive in order to obtain the maximum percentage 
error. The expression for the pressure drop, pi—f., will 
always be positive. 

An accuracy of +0.1-mm Hg may be readily attained 
in the range of pressure measurements that are dealt with. 
Hence, 6p; and 6p. may each be set equal to 0.1. An ac- 
curacy of +0.1 milliliter in Vg may also be easily realized, 
since a calibration of the buret space with mercury is 
possible before the apparatus is assembled. Therefore, 
5Vzg=0.1. The dead space volume is generally determined 
by a calibration with helium. A reasonable value in the 
accuracy of such a measurement is 1 percent. The value 
for 5Vp may be determined, consequently, by multiplying 
the value for Vp by 0.01. 

The data contained in Table I were observed in a num- 
ber of different adsorption tubes. These calculations were 
made for the first point of each isotherm. The calculation 
for subsequent points becomes more complicated, and 
the equation above is not directly applicable. It is obvious 
that the error in obtaining subsequent points is cumulative. 
Nevertheless, it is felt that the treatment above is in- 
dicative of the relative accuracy of this type of measure- 
ment. 

It should be noted that the absolute accuracy of the 
surface determination includes other factors than those 
treated here. These other quantities are dependent on the 
validity of the assumptions in the B.E.T. theory and on 
the density of packing of the adsorbed nitrogen mole- 
cules. Since the dead space is obtained by a helium cali- 
bration, the above consideration shows that it is necessary 
to make this calibration with an accuracy adequate to the 
particular application. 





The Formation of Per-Compounds in 
Relation to Electronegativities 


A. D. WALSH 
Laboratory of Physical Chemistry, Cambridge, England 
July 3, 1947 


I is pertinent to add to the letter of M. Haissinsky,' 

under the above title, the following remarks. They are 
largely taken from recent communications? in which we 
have surveyed some of the many factors that determine 
bond strengths. 

It is clear that the electronegativity of the groups at 
either end of a bond is a very important factor in de- 
termining the strength of the bond. The bond strength 
commonly increases with the electronegativity of the 
bonded groups. Examples are numerous** and the point 
has recently been stressed by Gordy* who has given an 
empirical relation between bond strength and the product 
of the electronegativities of the groups concerned. It is 
equally clear, however, that the increase of bond strength 
with electronegativity of bonded groups can only be 
true up to a certain point: if the electronegativities are 
too great there will be insufficient electron density in the 
‘shared’ region between the nuclei that is all important 
for bonding purposes. In other words, high bond strength 
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also requires adequate atomic orbital overlap, and too 
high an electronegativity of the bonded groups may result 
in reduced overlap. We therefore expect the rule of in- 
creasing bond strength with increasing electronegativity 
product to break down when we come to consider bonds 
between highly electronegative groups. Thus it is not sur- 
prising that Gordy’s relation fails for the OO bond in 
H,2O02 and probably for the FF bond.? As regards the OO 
bond the indications® * are that the strength increases with 
substitution of less electronegative groups for H. Simi- 
larly, the FF bond is weaker than the FCI bond, FF is 
weaker than FH, OO weaker than OH, and NN weaker 
than NH. If we accept the structure (J) advocated by 
Ingold and Ingold® for N2O,, the dissociation energy of 
the NN bond is 13 kcal./mole. The bond energy (Z) must 
be distinguished from the dissociation energy”® but is 
likely also to be very low since the stretching force constant 
(k) is only 1.5105 dyne/cm. In hydrazine,? however, 
E(NN) appears -to be ~64 kcal./mole, with k(NN) 
= 3.6105 dyne/cm. The great reduction of NN strength 
in N2O, relative to N2O, can be explained as caused by 
reduced overlap in the NN bond of N2O, consequent 


O O 
N-N (I) 
O O 


upon the high dipole ca ) in each NO, group. Glyoxal 
affords a similar example,? the CC bond strength appearing 
to be considerably weaker than considerations of the con- 
jugation between its two C =O groups alone would indicate. 

The conclusion of M. Haissinsky that elements whose 
electronegativities are smaller than a certain amount 
form per-compounds while those whose electronegativities 
are greater than this amount do not form per-compounds 
finds a natural explanation in terms of this factor of re- 
duced atomic orbital overlap consequent upon tco great 
an electronegativity. It is an important point of peroxide 
chemistry that the OO bond, to be strong, needs groups of 
low electronegativity attached to it. 

1M. Haissinsky, J. Chem. Phys. 15, 152 (1947). 

2? Walsh, J. Chem. Soc., in press (2 papers). 

*W. Gordy, J. Chem. Phys. 14, 305 (1946). 

‘Walsh, Trans. Faraday Soc. 42, 264 (1946). 


5 Ingold and Ingold, Nature, 159, 743 (1947). 
* Walsh, Trans. Faraday Soc. 43, 60 (1947). 





Addendum: A Relation between Bond Order 
and Covalent Bond Distance 


[J. Chem. Phys..15, 284 (1947)] 


H. J. BERNSTEIN 
Division of Chemistry, National Research Council, 
Ottawa, Canada 


HE paper, by G. N. Copley,' presenting the same 

model for obtaining multiple carbon bond distances 
as that which appeared in the appendix of this paper, has 
just come to my attention. Copley treats CC, CN, CO 
bonds and requires an empirical correction factor to his 
relation between bond order and bond distance. 


1G. N. Copley, Chem. Ind..663 (1941), 
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Structure-Optical Studies 


I. X-Ray Diffraction by Addition Compounds of 
Halogens with Hydrophilic Organic Polymers 
C. D. West 


Polaroid Corporation, Cambridge, Massachusetts 
July 14, 1947 


ODINE, bromine, and iodine bromide can be incorpor- 
ated in high concentrations in films or fibers of many 
hydrophilic organic high polymers, crystalline and amor- 
phous, natural and artificial. Examples are starch, cellulose, 
polyvinyl alcohol (PVA), the polyamides called Nylon (not 
previously reported), and many of their derivatives. These 
polymers are all transparent in the visible, while the 
halogens by themselves absorb light in this region. It is 
easy to show that an increase in the absorption coefficient 
of the halogen generally accompanies its incorporation into 
the polymer, indeed a well-defined new, strong band often 
appears in the visible (blue starch-iodine). By a high con- 
centration of, e.g., iodine on PVA, is meant about 30 per- 
cent by weight of free iodine, accompanied by an amount 
of hydrogen or alkali iodine which may be varied between 
wide limits. The presence of halogen in these preparations 
in the known orthorhombic, or any other three-dimensional 
crystal form, is readily excluded by x-ray diffraction study. 
The polymers named (save starch), including such sub- 
stantially amorphous PVA derivatives as the borate and 
the commercial butyral can be put into well-oriented uni- 
axial states by methods involving unidirectional stretch- 
ing, if they do not naturally occur in this state. When the 
halogen is incorporated into such well-oriented prepara- 
tions, then strong uniaxial positive dichroism accompanies 
the enhanced light absorption.! 

We have photographed, by monochromatic x-ray diffrac- 
tion under widely varying conditions, many such dichroic 
preparations and find, nearly without exception, a single 
constant characteristic feature on the pictures: a pair of 
diffuse continuous layer lines, which are also observable in 
the second and third orders under favorable conditions. 
This can only indicate diffraction by independent line 
gratings, which are here extended straight chains of halogen 
atoms all lying parallel to the orientation or absorption 
direction. The periodicities of the halogen chains are sub- 
stantially independent of the original periodicity of the 
oriented organic polymer (if any) and are for iodine 3.10, 
for bromine 2.70, for iodine bromide 2.90A, with individual 
values possibly varying up to 0.15A from these figures. 
Besides these characteristic diffractions, accessory diffrac- 
tions in other directions and varying to some extent with 
the mode of preparation are not uncommon, e.g., with 
PVA-iodine; these, however, need not affect the present 
arguments. Only rarely are the characteristic diffractions 
absent when they would be expected, such exceptions could 
be accounted for by picturing the halogen chains as short, 
or crooked, or both. With unoriented preparations, includ- 
ing those of starch-iodine, a diffuse halo corresponding in 
Position to the first characteristic continuous layer line can 
usually be seen; by themselves such haloes can prove little 
in the present problem, but they do have a certain con- 
firmatory interest. 


TO 
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The foregoing periodicities are much shorter than the 
unbonded separations of halogen atoms (4.3A for I—I, 
etc.) and are somewhat greater than the single bond 
separations (2.7A for I—I, etc.). Shortened distances be- 
tween known atom pairs, AB, usually mean either a hydro- 
gen bond between A and B, or the bonding of both A and 
B to a central atom C, or a direct bond of some kind be- 
tween A and B. The foregoing data suggest that only the 
third possibility comes into question here; the halogens in 
these preparations form bonds with each other at 180°, 
weaker than single covalent bonds but much stronger than 
van der Waals attractions, they are therefore properly 
called linear polymeric halogens or polyhalogens for short. 
Such a picture is consistent with the observed small fu- 
gacity of halogen in these preparations, and with the ob- 
served optical properties. Concerning the probable mode 
of interaction (if any) between halogens and oxygen or 
nitrogen atoms in the organic polymer, which confer hydro- 
philic properties on the latter, one can for the present only 
speculate within two limitations: first, such interaction 
cannot depend on any crystalline or long-range order of 
the ligher atoms; second, the incorporation of the halogen 
does not affect the infra-red absorption spectrum of the 
organic polymer over the range 3-15 mu, so far as we are 
able to determine. 

The several previous explanations for these addition com- 
pounds, of which the latest was made by Rundle and co- 
workers,” require revision to account for the characteristic 
x-ray diffractions which we consider to be the single most 
important clue to their structure. 

1E. H. Land and C. D. West in J. Alexander, Colloid Chemistry 
(Reinhold Publishing Corporation, New York, 1946), Volume 6, 
Chapter 6, pp. 160-90. 


2 R. Rundle, J. Foster, and R. Baldwin, J. Am. Chem. Soc. 66, 2116 
(1944). 





Activation in Unimolecular Reactions* ** 


O. K. Rice! 
Clinton Laboratories, Oak Ridge, Tennessee 
July 9, 1947 


HE decomposition of nitrogen pentoxide remains 

first order to extremely low pressures, of the order 
of 0.01 mm. Kassel? has shown that to account for this it is 
necessary, using the most probable model for the system 
of oscillators comprising the molecule, to suppose that 
the diameter for collisional activation and deactivation is 
about 20 x 10-8 cm. Recently Ogg** has made an ingenious 
suggestion which appears to account for this characteristic 
of the reaction, as well as a number of other features, 
without any excessive collision diameters. It is supposed 
that the equilibrium 


N.O;2NO;+NO, 
is established, and that the dissociated products of this 
equilibrium can occasionally react 
NO;+NO.—NO;+0:+NO 


to give the decomposition. Working out this mechanism 
shows that the reaction should be first order to arbitrarily 
low pressures; the apparent observed falling off in the 
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rate constant at very low pressures presumably indicates 
that an appreciable fraction of the N.O; is dissociated. 

A brief historical note, which will indicate where this 
explanation stands in relation to the theory of unimo- 
lecular reactions, may be in order. In the early days of 
the discussion of these reactions it was supposed that the 
equilibrium fraction of activated molecules would be 
given by the simple exponential e~@/"7, where Q is the 
activation energy. The rate of activation could not be 
faster than the rate of deactivation which would occur if 
there were no reaction and consequently no draining 
away of activated molecules. This in turn would be equal 
to the number of collisions made by the activated fraction 
of molecules, e~@/*", with other molecules. It was in 
general impossible, on this basis, to account for a suffi- 
ciently rapid rate of activation to maintain the reaction. 
It was soon pointed out,* however, that in the case of a 
molecule with many internal degrees of freedom, the 
fraction of molecules with energy greater than Q is much 
greater than e~@/87, All this energy may be considered as 
activation energy if the molecule is capable of trans- 
ferring energy between the internal degrees of freedom, 
and in particular to the bond which is to break in the 
reaction. Thus it is possible to account for much larger 
rates of activation. This comes about because there are 
many ways in which the excess energy necessary for ac- 
tivation can be distributed among the various oscillators 
of the molecule, thus increasing the probability, or the 
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entropy, of the activated molecules. Anything which in- 
creases this entropy may increase the rate of activation. 
For example, if the energy levels of an activated molecule 
are much closer together on the average than those in the 
normal state of the molecule, the rate of activation is 
increased. This can be important when the number of 
degrees of freedom is small.4 

It is now apparent how Ogg’s suggested mechanism 
fits into the scheme. The activated state of the molecule 
is actually the separated pair NO2+NO; which has an 
extremely high entropy. The rate of activation can, there- 
fore, be large. In fact, the difference in entropy between 
the activated and normal states increases with decreasing 
pressure in such a way as to render the rate of activation 
independent of the pressure. But the assumption that 
activation occurs while the parts of the molecule are 
separated is essentially a device for increasing the rate 
of activation. 

* This document is based on work performed under Contract No. 
W-35-058-eng-71 for the Atomic Energy. Project at the Clinton 
Laboratories. 

** A remark on the note, ‘‘The Mechanism of Nitrogen Pentoxide 
Decomposition,’”’ by R. A. Ogg, Jr., J. Chem. Phys. 15, 337 (1947). 

1On leave 1946-47 from the University of North Carolina. 

2L. S. Kassel, The Kinetics of Homogeneous Gas Reactions (Reinhold 
Publishing Corporation, New York, 1932) p. 182. 

3G. N. Lewis and D. F. Smith, J. Am. Chem. Soc. 47, 1514 (1925); 
J. A. Christiansen, Proc. Camb. Phil. Soc. 23, 438 (1926); C. N. Hinshel- 
wood, Proc. Roy. Soc. London 113A, 230 (1926); O. K. Rice and H. C. 
Ramsperger, J. Am. Chem. Soc. 49, 1617 (1927); see also Kassel, 


reference 3, p. 94. 
40. K Rice, J. Chem. Phys. 9, 258 (1941). 





